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A  fuel  mass  flowmeter  wa^developed  with  the  objective  of  accurately  meal- 
suring  the  mass  flow  rate'M)of  various  fuels  and  their  mixtures,  for  flow  ! 
rates  from  91  to  11, 364  kg/hr  (200  to  25,  000  lb/ hr ).  The  delivered  sys¬ 
tem  included  (1  )  an  area  averaging  zigzag  flow  velocimeter  cell,  in  which 
the  fuel  was  interrogated  in  opposite  directions  by  an  alternating  sgguence 
of  ultrasonic  pulses,  to  determine  the  area  averaged  flow  velocity  Y\ 

(2)  a  densitometer  based  on  the  Clausius-Mosotti  law,  in 
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cwhich  th/e  fuel  density  p  is  determined;  (3)  electronics  which  computed 
p'  andivl  separately,  aftd  then  generated  a  square  wave  output  having  a 
frequency  proportional  to  the  product  of  and  V,  namely,  M  =  pVA, 
where  A  =  velocimeter  cell  cross  Sectional  area.  *^The  transducers  can 
operate  beyond  the  temperature  range  -38  to  +150  C,  and  the  flow  cell, 
with  transducers,  at  pressures  up  to  1500  psig  and  in  the  3-axis,  20  g, 

5  to  2000  Hz  vibration  environment  specified  by  MIL  Std.  810-514.  1, 

except  for  noise  problems  at  a  few  discrete  fr equencies. ^Ultrasonic _ 

and  turbine  flowmeter  data  were  compared  using  a  fuel  substitute 
(Stoddard  solvent  _)_for  M  between  ~  43  and  14,  300  kg/hr.  In  water  cali¬ 
bration  tests  ,  at  V  from  0.  1  to  38  ft/s  (water  M  from  — -150  to  60,  000 
lb/ hr )  [30  mm  /  s  to  11.6  m/s,  or  ~70  to  27,  000  kg /hr]  flow  velocity, 
linearity,  and  accuracy  of  better  than  0.  5%  full  scale  was  observed 
using  the  equipment  delivered  at  the  end  of  the  program.  To  determine 
repeatability,  a  10-day  test  was  conducted  with  water  flowing  through 
an  offset  cell  of  6  mm  inside  diameter,  at  a  velocity  between  0.  5  and 
4.  4  ft/s  (~-  1  50  and  ~  1  500  mm/ s  ).  The  velocity  was  controlled  by  the 
pressure  of  a  variable  height  water  column.  Repeatability  (one  stan¬ 
dard  deviation)  was  0.  2%  of  reading  at  ~  1500  mm/s  and  2%  of  reading 
at  -  150  mm/s.  In  the  final  Stoddard  solvent  (7024  B  II)  tests,  the  rms 
deviation  from  linearity  was£l.7%  ofthe  reading  for  flow  rates  above 
5  %  of  the  maximum  M  encountered  in  these  tests,  14,300  kg/hr  (31,450 
lb/hr).  Based  on  a  full  scale  rating  of  11,  364  kg/hr  (25,  000  lb /hr)  rms 
deviation  from  linearity  was  demonstrated  to  be  about  1%  of  full  scale, 
with  a  2  s  response  time.  Up  to  0.  5%  of  this  nonlinearity  may  be  due  to 
the  turbine  sensors'  own  error  band.  Above  9000  lb/hr,  for  3  runs,  the 
average  nonlinearity  is  less  than  0.  77%  of  the  actual  flow  rate.  The 
stability  of  the  waveforms  observed  for  the  zigzag  flow  cell  at  the  maxi¬ 
mum  flow  rate  encountered  in  the  fuel  substitute  test  (14,  300  kg/hr,  or 
31, 450  lb/hr)  and  in  the  water  tests  (V  =  1 1.  6  m/ s,  or  38  ft/ s)  suggests 
that  a  longer  zigzag  path  could  be  used  in  the  future.  This  would  im¬ 
prove  the  resolution  at  low  rates  and,  together  with  some  other  modifi¬ 
cations  to  the  cell  and  transducer,  ought  to  improve  linearity  in  the 
lower  part  of  the  system's  range.  Alternatively,  linearization  may  be 
accomplished  in  the  software,  at  least  for  nonlinearities  mainly  depend¬ 
ent  on  V. 
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I.  INTRODUCTION 


A.  Need  for  Mass  Flow  Rate  Measurements 

The  purpose  of  the  proposed  development  effort  was  to  provide 
NAVAIRDEVCEN  with  a  highly  accurate,  rugged  and  reliable  fuel  mass 
flowmeter.  The  ideal  flowmeter  would  be  sufficiently  versatile  so  that 
it  could  be  employed  for  engine  diagnostics  in  future  engines  of  the  type  an¬ 
ticipated  for  advanced  type  A  aircraft  (subsonic)  and  later,  type  B  aircraft 
(supersonic).  Accordingly,  a  flowmeter  was  sought  which  had  the  potential 
for  being  optimized  with  respect  to  many  capabilities  and  characteristics, 
including  the  following: 

o  Accurate  measurement  of  mass  flow  rate  of  JP-4,  JP-5 
and  their  mixtures 

o  Repeatability 

o  Hysteresis 

o  Vibration  (cell,  electronics) 

o  Rangeability  >  250: 1 ,  possibly  ~  800: 1 

o  Flow  cell  materials,  shape,  dimensions  and  weight;  quick 
disconnect  type  fluid  connectors;  and  cell  mounting 
means 

o  Response  time  =  2s 

o  Square  wave  (digital)  signal  output,  capable  of  multiplex 
operation 

o  Contaminated  fuel:  operation  without  clogging  or  degrading 

o  Temperature  =  -38  C  to  +150°C  (fuel),  and  even  higher  soak- 
back  temperature  (est.  200°C) 

o  Pressures  to  1600  psig 

o  Laminar,  transitional  and  turbulent  flow  profiles 

o  Minimum  number  of  conductors  between  flow  cell  ar.d  cockpit- 
installed  electronics  (<  4  desired) 

o  Power  source:  28V  dc;  design  shall  include  115V,  400  Hz  option 

o  Totalized  volumetric  (Q),  or  mass  (M)  for  convenience  in  trouble¬ 
shooting  and  calibration 

These  needs  were  expressed  by  NAVAIRDEVCEN  in  1977,  and  were 
addressed  during  the  course  of  the  contract.  The  relevance  of  these  needs 
is  apparent  if  one  compares  them  (Section  I-D)  with  needs  expressed  three 
years  later  by  other  organizations  (e.  g.  ,  NASA)  concerned  with  gas  turbine 
engine  diagnostics  and/or  control. 
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Many  types  of  flowmeters  are  presently  available,  such  as,  magnetic, 
fluidic,  nuclear  magnetic  resonance,  turbine,  vortex  shedder,  pitot  tube, 
orifice  plate,  coriolis,  venturi,  target.,  ultrasonic,  etc.  Some  of  these 
would  have  to  be  combined  with  a  densitometer  to  yield  mass  flow  rate. 
(Dowdell,  1974;  Sturdivant  etal.  ,  1978).  However,  if  one  seeks  a  mass  flow¬ 
meter  which  has  no  moving  parts,  minimum  obstruction  to  the  flow,  small 
size  and  light  weight,  wide  rangeability,  etc.  (see  Appendix  A)  the  choice 
narrows  considerably,  and  in  fact  at  the  program' s  outset  in  Dec.  1977  it  was 
obvious  that  no  equipment  that  was  commercially  available  could  satisfy  all 
the  requirements  and  objectives.  Some  of  the  objectives  are  still  unattained, 
partly  due  to  funding  and  schedule  constraints,  and  partly  due  to  the  state  of 
the  art. 


Besides  at  least  a  twenty-five  year  diagnostic  need  for  measuring 
mass  flow  rate  in  military  and  commercial  aircraft  engines  (Kritz,  1955), 
there  remains  a  control  need  which  differs  somewhat  from  the  diagnostic 
need.  While  diagnostic  needs  exist  both  on  the  ground  (engine  test  stands) 
and  in  the  air,  control  needs  are  obviously  most  important  in-flight.  For 
control  applications,  accuracy  may  sometimes  be  relaxed  slightly,  perhaps 
to  +_  2%  of  range,  but  response  time  is  much  more  stringent,  20  ms  sometimes 
being  required,  i.  e.  ,  100  times  faster  than  the  present  contract's  diagnostic 
response  time  objective. 

Aside  from  the  need  for  accurate  M  data  per  se,  a  totalized  or  time- 
integrated  measurement  of  j£  Mdt  yields  the  mass  of  fuel  burned,  from  which 
fuel  remaining  can  be  calculated.  In  other  words,  a  fuel  mass  flowmeter 
can  serve  as  a  fuel  quantity  instrument. 

B.  Background  on  Related  Work 

In  April  1971  Panametrics  responded  to  the  Army's  Eustis  Directorate 
advertised  RFP  for  improved  diagnostic  instrumentation  for  aircraft  engines. 

On  28  June  1971  Eustis  awarded  Panametrics  our  first  contract  to  develop  an 
ultrasonic  fuel  mass  flowmeter. 

Our  original  proposal  offered  a  unique  combination  of  a  doppler  flow 
velocimeter  (the  doppler  shift  is  proportional  to  V/c,  where  V  =  flow  velocity 
and  c  =  sound  speed  in  the  fuel)  and  an  acoustic  impedance  (pc)  densitometer. 

The  product  of  (V/c)  (pc)  =  M/A,  the  mass  flow  rate  per  unit  area.  As  our 
tests  and  scattering  analysis  proceeded,  it  became  evident  from  the  data  that  clean  • 
fuel  could  not  scatter  enough  sound  to  enable  a  fast  response  ultrasonic  doppler 
system  to  provide  accuracy  better  than  a  few  percent,  say  ~  5%.  N.  E.  Pedersen 
soon  found  an  alternate  way  to  measure  flow  velocity.  This  alternate  way 
involved  a  pseudo  random  noise  (PN)  code  which  phase-modulated  5  MHz  cw 
waves.  The  area-averaging  flow  cell  used  in  this  first  Eustis  program  is 
shown  in  Fig.  1,  exploded  view,  and  Fig.  2  shows  this  cell  assembled,  and  also 
shows  the  PN  flowmeter  block  diagram. 
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We  have  continued  to  use  modulated  cw  in  some  of  our  flowmeters, 
but  for  several  technical  reasons  we  now  use  amplitude  modulation  rather 
than  phase  modulation.  The  amplitude  modulation  method,  now  patented,  is  des¬ 
cribed  by  its  inventors,  Pedersen  and  Bradshaw  in  their  patent  (1977)  and  in  the 
Proc.  of  the  1977  NBS/AIP/ISA  Symp.  on  Flow.  Our  second  Eustis  contract,  cov¬ 
ering  the  period  January  1973  to  June  1975,  supported  the  fabrication  and  testing 
of  an  amplitude -modulated  flow  velocimeter,  i.  e.  ,  the  use  of  coherent  rf 
bursts.  During  this  second  contract,  we  learned  of  the  capabilities  of 
Simmonds  Precision's  (SP)  capacitance -type  densitometer.  Test  data  made 
available  to  us,  plus  published  accuracy  data  for  JP-4  and  JP-5  (Stuart,  1974) 
showed  that  the  capacitance  densitometer  was  likely  to  be  more  accurate  than 
our  pc  densitometer  available  at  that  time.  We  therefore  obtained  on  loan 
one  of  SP' s  densitometers  and  used  it  in  calibration  tests  at  Avco  Lycoming 
in  December  1974.  Precision  of  the  SP  densitometer  was  found  to  be  0.  1% 
or  better.  Precision  of  the  entire  M  system  (Panametrics  ultrasonic  flow 
velocimeter  plus  SP  dielectric  densitometer)  was  0.  25%  FS,  and  linearity 
~  0.  1%  FS. 

In  view  of  these  results,  in  June  1976  we  were  awarded  a  third  contract 
from  Eustis  to  particularize  the  design  for  GE's  T-700  engine,  for  the  M  range 
1  00  to  1 000  lb/hr. 

Fabrication  and  testing  of  the  5  MHz  Eustis  electronic  design  shown 
in  Figs.  3  and  4a,  and  the  cell  of  Fig.  4b  was  completed  shortly  before  the  present 
Navy  program  began  in  Dec.  1977. 

The  status  of  the  Eustis  work  at  the  end  of  the  third  contract  was 
reported  by  Lynnworth,  Pedersen,  Seger  and  Bradshaw  (1978).  This  report  des¬ 
cribed  in  detail  how  the  fuel  mass  flowmeter  system  (previously  demonstrated  in 
the  second  Eustis  program,  and  consisting  of  an  ultrasonic  flow  velocimeter  and 
a  capacitance-type  densitometer)  was  adapted  to  T-700  gas  turbine  engine  re¬ 
quirements.  Mechanical  adaptations  included  the  reduction  in  size  and  weight 
of  the  flow  cell  to  fit  the  configuration  and  mounting  constraints  of  a  T-700 
engine.  Electronic  developments  included  repackaging  into  a  portable  field- 
type  case,  configuration  switching  and  digital  delay  equalizers,  and  the  sub¬ 
stitution  of  a  microprocessor  for  a  previously -used calculator  chip,  resulting 
in  a  choice  of  response  times,  or  integration  times,  selectable  from  ~  50s 
down  to  ~  0.  5s.  Test  results  were  presented  for  7024  Type  IIB  calibrating 
fluid  over  the  mass  flow  rate  range  M  =  40  to  1200  PPH,  for  temperature  near 
70°F  (2l°C),  and  for  pressures  up  to  1500  psig.  Performance  of  the  system 
was  discussed  with  respect  to  laminar,  transitional,  and  turbulent  flow  profiles 
associated  with  Reynolds  numbers  from  ~  50  to  ~  50,  000  and  with  respect  to 
geometrical  constraints  imposed  by  the  T-700  configuration. 
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C.  Distinctions  Over  Prior  Work 

At  the  time  the  present  work  was  proposed,  July  1977,  it  was  planned 
to  modify  the  Eustis  system  to  meet  the  Navy  requirements.  The  differences 
in  requirements  and  objectives  between  Eustis  and  NAVAIRDEVCEN  are 
summarized  in  Table  1.  Shortly  after  the  time  the  Navy  program  got  under¬ 
way,  in  early  1978,  it  became  apparent  to  the  principal  investigators  that 
a  new  phase  locked  loop  electronic  approach  quite  distinct  from  the  Eustis 
rf  burst  system  merited  consideration.  Accordingly,  both  of  these  approaches 
are  reviewed  in  section  II B.  See  also.  Appendix  C. 

Three  years  later,  in  July  1980,  the  new  broadband  phase  locked  loop 
(PLL)  technique  had  been  demonstrated  over  a  wider  range  of  flow  rates, 
with  numerous  fluids  and  with  a  wider  variety  of  transducer  and  flow  cell 
combinations,  than  had  been  achieved  with  the  rf  burst  or  any  other  earlier 
Panametrics  flowmeter  approach.  The  PLL  technique,  however,  being  broad¬ 
band,  inherently  has  resolution  and  acoustic  noise  limitations  which  are  not 
as  good  as  the  rf  burst  method.  For  example,  the  present  PLL  resolution 
limit  is  1  ns,  as  compared  to  the  rf  burst  resolution  of  0.  2ns  achieved  pre¬ 
viously  at  a  carrier  frequency  of  5  MHz.  The  PLL  bandwidth  typically  ex¬ 
ceeds  100%  of  the  center  frequency;  the  rf  burst  bandwidth  is  typically  10%. 

In  retrospect,  considering  the  various  tradeoffs  between  the  PLL,  rf 
burst  and  other  methods,  we  conclude  that  the  PLL  method  was  indeed  the  most 
appropriate  method  to  use  in  this  contract.  The  performance  achieved  in  this 
contract  represents  results  on  two  area  averaging  flow  cells  (zigzag  path  in 
square  holed  sleeve,  and  offset)  with  the  first  PLL  velocimeter  of  the 
"Model  6000"  type.  It  is  reasonable  to  expect  that  improvements  in  the  flow 
cell,  the  transducers,  and  the  electronics  (hardware,  firmware  and  software) 
will  lead  to  improved  performance  in  the  future. 
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Table  i.  Comparison  of  Eustis  and  Navy  Requirements  and  Objectives 

REQT  OR 

OBJECTIVE 

EUSTIS 

NAVAIRDEVCEN 

F  uels 

JP-4,  -5  mixtures; 
Contaminated 

Same  types  of  fuel(s) 

M,  lb/hr 

100- 1000 

200-25, 000  (SO, 000) 

Accuracy 

+  0.  5%  FS 

+  0.  5%  R,  3000-25000  PPH; 

+  5%  R  at  100  PPH 

Hysteresis 

- 

0.  3%  of  reading  (R) 

Repeatability 

- 

+  0.  1%  R,  3000-25000  PPH 
+  2%  R,  200  PPH 

Cal.  Temp.  Range,  °C 

+  10  to  +  46 

-  38  to  +  150 

Response  Time,  sec. 

5 

2 

Signal  Output  F ormat 

- 

Sq.  wave,  f  or  M;  muxable 

Vibration 

810-C,  Proc.  I** 

810-514. 1 

Power 

115V  60  Hz 

28V  dc;  optional  design, 
115V,  400  Hz 

Wiring 

- 

4  conductors,  2  pwr,  2  sig; 
optical  option 

Flow  connectors 

1500  psi,  AN-3,  -4 

1600  psi,  quick  disconnect 

Configuration 

T -700 

V/STOL  Type  A(B) 

Electronic  housing 

Field  type  case  for 
helicopter  cabin 

Flow  cell  weight,  lb 

<  5  lb  if  SS,  <  1  lb 
if  aluminum 

Delivery  date  for 

October  1977 

October  1979 

complete  system 

*4  hr  exposure;  calibrate  before  and  after. 

Navy,  300  hr;  calibrate  during 

contaminated;  fuel  contaminated  according 

to  Table  X,  MIL-E-5007D. 

**No  flow,  room  temperature. 

-  ..  .  . 
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D.  Comparison  with  Present  Objectives  of  Other  Organizations 

The  continuing  need  for  improved  M  instrumentation  is  reflected 
in  needs  expressed  at  present  by  government  and  industry.  Because  of 
the  growing  scarcity  of  fuel,  and  its  rising  cost,  even  minor  improvements 
in  engine  efficiency  are  significant.  Hence  manufacturers  of  jet  engines 
seek  M  accuracy  of  0. 1%  (i-  e.  ,  error  <  0. 1%  of  reading). 

NASA  RFP  3-188129,  issued  in  early  1980,  reiterated  the  need  for 
M  measurements  in  fuels  such  as  JP-3,  -4,  -5,  -8  (Type  A-l)  and  Type  A 
with  rangeability  on  the  order  of  100:1,  an  accuracy  requirement  of  0.  25% 
of  reading  with  0.  1%  of  reading  as  a  goal,  and  a  response  time  requirement 
of  25  ms.  These  and  other  present  design  guidelines  and  specifications  are 
contained  in  Appendix  F  for  reference  and  comparison  purposes. 

E.  Summary  of  Theory 

The  theory  of  the  present  M  system  may  be  explained  in  five  steps: 

1.  The  mass  flow  rate  M  is  proportional  to  the  product  of  the  average  flow 
velocity  V  and  the  fuel  density  p. 

2.  V  is  measured  ultrasonically  by  determining  the  small  difference  in 
transit  times  in  the  upstream  and  downstream  directions.  Ideally, 

100%  of  the  flowing  cross  section  is  interrogated,  so  that  the  transit 
times  provide  area-averaged  values  independent  of  the  flow  profile. 

3.  p  is  measured  dieiectrically  (capacitively)  based  on  the  Clausius-Mosotti 
law  that  relates  p  to  e,  the  fuel's  dielectric  constant.  A  Simmonds 
Precision  densitometer  provides  an  analog  voltage  which  is  empirically 
related  to  p. 

4.  A  microprocessor  -  controlled  computer  performs  the  necessary  cal¬ 
culations  to  yield  M. 

5.  A  square  wave  is  generated  at  a  frequency  proportional  to  M. 

These  steps  are  covered  in  detail  elsewhere  in  this  report,  and  in 
the  literature.  It  may  be  useful,  however,  to  expand  on  the  V  and  p  ideas 
in  the  following  paragraphs. 
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For  the  contrapropagating  transmission  method  of  flow  measurement 
utilized  in  this  program,  transit  time  decreases  in  the  downstream  direction 
and  increases  upstream.  Flow  rate  is  proportional  to  (t^  -  w^erc 

the  values  of  t  are  measured  transit  times.  A  valid  approximation  when 
(V/c)^<<  1  is  given  by  an  equation  of  the  form 

V  = 

where  P-  length  of  liquid  path  and  L-  axial  projection  of  path  in  flowing  liquid. 

To  retain  high  accuracy  in  V  despite  variable  llow  profiles  it  is  neces¬ 
sary  to  calibrate  and/or  use  equipment  which  properly  weights  the  flow  profile. 
The  delivered  equipment  utilizes  an  area  averaging  cell.  In  contrast  to  this 
method,  meter  factors  for  conventional  flow  cells  are  usually  based  on  the 
Reynolds  number  Re.  Theoretically,  the  meter  factor  K  =  1  for  cells  with 
mid  radius  paths  or  cells  suitably  interrogated  over  100%  of  their  area,  i  or 
large  diameter  pipes  interrogated  conventionally  (tilted  diameter  path)  K—  0.750 
for  laminar  flow  (Re<  2000),  K  »  0.85  for  transitional  flow  (2000<  Re<  4000) 
and  K  I  '(1.119-  0.01  1  log  Re )  for  turbulent  flow  (Re  >  4000).  For  example, 

K  =  0.930  at  Re  =  104,  0.940  at  105,  0.950  at  106,  0.960  at  107. 


The  electronic  approach  to  measuring  V  utilizes  a  modified  Panametrics 
Model  6000  instrument.  This  instrument  transmits  pulses  sequentially 
upstream  and  downstream  at  two  phase  locked  repetition  frequencies  whose 
periods  are  proportional  to  the  times  of  flight  tj  and  betw'een  transducers. 

By  timing  a  large  number  of  periods,  t  and  ty  are  measurable  to  high  pre¬ 
cision  (resolution  -  Ins).  The  flow'  velocity  V  is  computed  by  an  equation 
similar  to  that  given  above.  (See  Appendix  C.)  The  adjustable  interrogation 
rate  is  factory- wired  to  avoid  errors  due  to  reverberations  in  transducers,  pipe 
wall  or  fluid.  The  number  of  interrogations  to  be  averaged  is  chosen  to  maxi¬ 
mize  resolution  within  the  allowed  response  time. 

The  dielectric  constant/density  relation  used  in  SP1  s  densitometer 
is  given  by  the  Clausius  -  Mosotti  law  (c  -  1  )/(e  +  2)  =  C p  .  This  is  rearranged 
to  yield  density  as  p  =  (e  -  1 )  /[  a(e  -  1 )  +  b]  where  a  and  b  are  empirical  constants. 
By  measuring  capacitance  between  concentric  tubes,  p  is  determined. 

Having  determined  p  and  V,  the  mass  flow  rate  M  is  determined  from 
their  product,  taking  into  account  the  cell  cross  sectional  area  A,  and  using 
consistent  units ,  of  course. 
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II.  DESIGN  CONSIDERATIONS  AND  ALTERNATIVES 

A.  Hydrodynamic  Considerations /Flow  Profile 

A  full  discussion  of  flow  profile  is  beyond  the  scope  of  this  report. 
The  interested  reader  is  referred  to  books  by  Schlichting  (1968)  or  Streeter 
(I96l)  for  tutorial  coverage  of  this  topic. 

For  our  present  purposes  it  may  suffice  to  simplify  or  abbreviate 
the  discussion  by  listing  the  following  relationships: 


Reynolds  Number 
<  2000 

2000  to  4000  (or  10000) 
>  4000  (or  > 10000) 


Flow  Regime 


Laminar 

Transitional 

Turbulent 


Flow  Profile 

Parabolic 
Unstable 
Power  law 


Reynolds  number  Re  =  pVD/q  where  p  =  density,  V  =  average  flow 
velocity,  D  =  conduit  diameter  and  n  =  viscosity.  The  ratio  rj/p  is  often 
denoted  v,  the  kinematic  viscosity. 

To  obtain  a  general  perspective  on  the  variation  of  v  with  temperature 
T  for  numerous  fluids,  and  therefore  a  perspective  on  how  T  influences  Re, 
refer  to  Fig.  5.  The  right-hand  portion  of  Fig.  5  yields  Re  in  terms  of  v  and 
VD.  Thus,  at  a  given  v,  say  v  =_0.  01  stokes  (water  at  room  temperature) 
if  D  =  2.  54  cm,  then  increasing  V  from  4  to  40  cm/s  will_increase  Re  from 
1000  to  10000,  i.  e.  ,  laminar  to  turbulent.  Similarly,  at  VD  =  100  cm^/s, 
if  v  decreases  from  0.  1  to  0.  01  stokes,  Re  again  ranges  from  1000  to  10000. 

According  to  Fig.  6c,  the  value  of  viscosities  of  interest  in  this 
program  ranges  from  about  0.  3  to  30  cs,  or  0.  003  to  0.  3  stokes.  According 
to  Fig.  7,  the  V  range,  for  a  1  in.  x  1  in.  duct  (25.  4  mm  x_25.  4  mm),  is  from 
1.  7  to  300  in.  / s  (4.  32  to  76_2  cm/ s).  Therefore  the  min  VD  =  4.  32  x  2.  54  = 

10.97  cm  / s,  and  the  max  VD  =  1935  cm^/s,  or  approximately  10  to  2000cm^/s. 
Refering  again  to  Fig.  5,  we  see  that  Re  ranges  from  ~  40  to  ~  7  x  10^.  Thus, 
flow  cell  and/or  analytic  (e.  g.  ,  p processor)  compensation  for  profile  is 
required  in  order  to  obtain  accurate  results. 

We  may  remark  that  in  the  course  of  contrapropagating  ultrasonic 
flowmeter  developments  since  the  1920's,  the  importance  of  flow  profile  has 
been  recognized  since  ~  195 0.  A  partial  list  of  investigators  who  published 
their  ideas  on  how  to  deal  with  flow  profile  is  given  in  Table  2.  These  con¬ 
tributions  are  fully  referenced,  and  several  are  illustrated,  in  Mason  and 
Thurston  (1979),  pp.  407-525. 


8 


NADC-80254-60 


The  present  contract  mainly  utilized  the  eighth  area  averaging 
method  of  Table  2.  The  principal  investigators  have  used  this  area  averaging 
method  in  three  previous  M  contracts  with  satisfactory  results.  Depending 
on  requirements,  however,  other  cell  designs  may  be  preferred  for  reasons 
of  size,  weight,  resolution,  ease  of  installation  and  maintenance,  absence 
of  cavities  in  which  gas  could  be  trapped,  or  cost.  See,  for  example.  Fig.  9, 
especially  offset  (IB)  and  clamp-on  (8A)  designs. 

B.  Electronic  Alternatives  Considered 

In  this  section  we  briefly  discuss  the  rf  burst  system  as  originally 
proposed,  e.  g.  ,  based  on  the  Eustis  system  reported  in  1978,  and  then  we 
briefly  discuss  the  newer  phase  locked  loop  (PLL)  approach.  By  early  1978 
the  PLL  flowmeter  approach  had  been  developed  as  a  Panametrics  product 
to  the  point  of  demonstrating  stability  and  reasonable  linearity.  By  May  1978 
it  became  apparent  that  the  PLL  was  likely  to  emerge  as  the  superior  approach 
to  meet  the  velocimetry  objectives,  and  in  fact  this  approach  was  eventually 
selected  and  utilized  in  the  delivered  mass  flowmeter  system. 


1  •  RF  Burst 

The  rf  burst  flow  velocimeter  due  to  Pedersen  and  Bradshaw  transmits 
coherent  rf  bursts  (  -VM  *  •  •  )  of  50  ps  duration  upstream,  and  then  down¬ 
stream.  Each  received  burst  is  narrowband  filtered.  The  phase  difference 
between  upstream  and  downstream  directions,  A0  ,  is  proportional  to  the 
flow  velocity  divided  by  c^.  In  the  Eustis  system,  to  eliminate  undesired 
effects  of  sound  speed  c,  we  also  filtered  at  the  am  frequency,  2  kHz,  and 
the  phase  of  the  2  kHz  component  provided  the  c  compensation.  (In  another 
version  of  this  system,  developed  for  measuring  the  flow  velocity  of  natural 
gas,  the  rf  burst  center  frequency  is  100  kHz.  See  Pedersen  etal.  ,  1977). 

A  pprocessor  combines  p,  V,  and  the  flow  cell  geometry  to  yield 
the  desired  mass  flow  rate  M,  in  the  Eustis  system. 

The  flow  velocity  V  may  be  shown  to  be  proportional  to  AT /(Tj  +  T ^ )^ 
where  AT  =  T^  -  T^  and  the  T1  s  are  transit  times  measured  upstream  and 
downstream. 

A  block  diagram  of  the  Eustis  fully  coherent  electronic  system  which 
measured  T  -  T^  and  T^  and  T^  is  shown  in  Fig.  3.  It  operated  as  follows. 
Flow  related  parameters  were  measured  by  using  a  fully  coherent  electronic 
system  which  operated  synchronously  from  a  5  MHz  crystal-controlled  oscillator. 
The  transmitted  waveform  was  generated  by  means  of  synchronously  gating  the 
output  of  the  5  MHz  oscillator  for  a  50  ps  interval  at  a  2  kHz  repetition  rate. 
Upstream  and  downstream  common-path  transmissions  were  consecutive  and 
spaced  250  ps  apart.  Independent  measurements  were  made  of  the  phase  dif¬ 
ference  between  the  upstream  and  downstream  5  MHz  spectral  components  of 
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Table  2.  Contributors  to  ultrasonic  area  averaging  flowmeters  using  contra- 

propagating  transmission  (upstream- downstream  transit  time)  methods. 


Inventor(s)  or 


Author  (s) 

Year(s) 

Approach 

Remarks 

Swengel 

1950- 

1956 

100% 

Waveguide  antennas  communi¬ 
cate  over  sheet-like  unfolded 
path;  fas  25  kHe;  rectangular 
duct  dimensions  up  to  ~5  x  8  m. 

Kritz 

1955  a, b 

K 

Tilted  diameter  path;  f.=  10MHz; 
pipe  diameter  =  100  mm;  graph  of 

K  vs  Re;  questioned  whether  a 
path  other  than  the  tilted  dia¬ 
meter  could  provide  better  im¬ 
munity  to  flow  profile. 

Katzenstein  & 
Katzenstein 

1961 

1 00%  axial 

Axial  interrogation  illustrated 
for  small  tube. 

Knapp 

1964 

4-  chord 

Proposed  (but  did  not  demonstrate) 
use  of  4- chord  Gaussian  quadra¬ 
ture  method. 

Noble 

1968 

100% 

Demonstrated  axial  interrogation 
method;  f  6  MHz;  acrylic  tube 
diameters  =  2,38  and  6.35  mm. 

Malone  & 

Whirl  ow 

1971 

4-chord 

Quadrature  methods  of  Gauss, 
Chebycheff  and  Lobatto. 

Lynnworth 

etal. 

1972- 

1979 

100%  zigzag 

Zigzag  interrogation  over  square - 
or  rectangular- enveloped  path;  f  = 
0.5  to  5  MHz;  duct  dimensions 
from  12.7  to  25,4  mm. 

Pfau 

1973 

~  Midradius  path 

Computed  optimal  distance  to  be 
0.493  of  the  radial  distance  for 
both  laminar  and  turbulent  flows. 

Lynch  and 

Brown 

1974 

K 

Suggested  K  be  modified  elec¬ 
tronically  as  a  function  of  flow 
velocity  V. 

Baker  & 
Thompson 

1975 

~  Midradius  path 

Errors  calculated  for  chordal 
paths  at  and  near  midradius; 
tested  in  air  at  f »  40  kHz  in 

1 00  mm  duct. 

Lynnworth 

1975, 

1977a 

1 00%  axial 

Transmitter  and  receiver  trans¬ 
ducers  in  acoustically  isolated 

parts  of  flow  cell. 
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Table  2.  (Cont1  d) 


Inventor(s)  or 


Author(s) 

Year(s) 

Approach 

Remarks 

Lynnworth 

1977b, 

1978 

~  Midradius 

Meter  factors  calculated  for 
narrow-beam  chordal  paths 
at  and  near  midradius;  folded 
path  along  inscribed  equilat¬ 
eral  triangle  (r  %  0.  5a). 

L\  nnworth 

1977 

100%  axial 

Reflectors  installed  in 
enlarged  bore. 

Johnson  et  al. 

1977a 

Tomography 

Analysis  of  multiple  chord 
data  allows  algebraic  recon¬ 
struction  of  three  dimensional 
fluid  flow. 

Pedersen  etal. 

1977 

K 

V  computed  from  V  and  Re 
based  on  power  law. 

Lynnworth 

1980 

Axial  sample 

Reflected  path  at  ~  3/4  of 
radial  distance  (r  =  0.  76a). 

Lynnworth  etal. 

1980 

K 

K  switch -selectable;  K  varied 
electronically  as  a  function 
of  V  or  Re. 

Nleisser 

1980, 

1981 

100%  axial 

Heat  meter  uses  fixed  offset 
geometry  with  cylindrical  tee 
manifold,  1  MHz  transducers. 

Pedersen 

1981 

Tomography 

Volumetric  flow  rate  deter¬ 
mined  in  conduit  of  arbitrary 
shape  by  integration  over 
conduit  area  using  linear 
equation  combinations. 

N  ADC -8  02  54- 60 


the  received  waveforms  and  the  phase  sum  of  the  2  kHz  modulation  components 
of  these  waveforms.  These  two  independent  measurements  were  combined 
to  yield  the  Mach  number  V /c  and  V,  where  V  =  flow  velocity  and  c  =  sound 
speed.  High  signal-to-noise  ratio  was  achieved  by  means  of  narrowband 
receivers  and  fully  coherent  detection,  thereby  providing  distint  advantages 
in  precision  and  stability  over  conventional  flowmetering  techniques  which 
utilize  broadband  pulses. 

One  of  the  limitations  on  the  single  frequency  (e.  g.  ,  5  MHz)  rf  burst 
system  as  just  described,  is  that  at  high  flows,  the  phase  shift  may  exceed 
180°  or  even  360°,  causing  ambiguity  unless  remedies  are  introduced.  If 
a  low  carrier  frequency  is  used,  to  prevent  A 0  from  exceeding  180°,  then 
at  low  flows,  adequate  resolution  of  A0  becomes  exceedingly  difficult.  The 
phase  ambiguity  problem  is  aggravated  by  the  present  contract's  wide  temperature 
range  which  leads  to  over  a  2:1  sound  speed  range  from  ^750  to  ~  1600  m/s. 

We  considered  using  two  different  carriers,  say  0.  5  and  5  MHz,  with 
some  overlap,  to  avoid  phase  ambiguity  at  high  flow  yet  retain  adequate  phase 
resolution  at  low  flows.  However,  implementation  of  this  dual-frequency 
rf  burst  velocimeter  appeared  more  difficult  than  the  PLL  method,  which 
avoids  these  problems  of  ambiguity  and  resolution. 

2.  Phase  Locked  Loop 

Tb-  phase  locked  loop  (PPL)  method  as  utilized  in  our  flow  velocimeter 
originated  with  N.  E.  Pedersen.  In  some  respects  it  resembles  singaround 
circuits  that  have  been  used  for  some  thirty  years  in  sound  velocimeters  and 
later,  flow  velocimeters.  The  basic  idea  in  these  systems  is  to  launch  up¬ 
stream  and  downstream  waves  at  repetition  frequencies  whose  difference  Af 
is  proportional  to  flow  velocity  V.  Pedersen' s  concept,  which  was  demon¬ 
strated  prior  to  its  use  in  this  program,  included  the  use  of  a  VCO  which,  at 
steady  flow,  would  normally  control  the  launching  of  pulses  at  frequencies  f^ 
and  f  having  periods  equal  to  twice  the  time  of  flight  in  the  two  directions. 

The  phase  lock  technique  assured  that  lock,  once  achieved,  would  maintain 
the  desired  relationship.  In  operation  two  quasi-dc  control  voltages  and 
corresponded  to  f  and  f  .  Their  difference  AV  was  an  analog  measure  of  the 
flow  velocity  V.  •the  initial  PLL  sequentially  interrogated  upstream  and 
downstream  for  equal  intervals,  e.  g.  ,  0.  1  s  each  way. 

Two  problems  with  the  early  model  PLL  flowmeter  were  lack  of 
interchangeability  of  nominally  identical  VCO' s,  and  repetition  frequencies 
that  in  some  cases  were  so  high  that  cell  reverberations  and/or  transducer 
ringdown  interfered  with  obtaining  high  precision.  These  problems  were 
solved  in  the  delivered  equipment  as  described  in  Section  UIB  and  Appendix  C. 

C.  Flow  Cell  Alternatives  Considered 

Figure  9  provides  a  general  perspective  on  flow  cell  alternatives  that 
may  be  considered  for  contrapropagating  transmission  ultrasonic  flowmeters. 

In  any  particular  application,  the  choice  is  usually  governed  by  pipe  size 
(flow  rate),  accuracy,  geometrical,  installation  and  maintenance  constraints, 
cost  and  delivery  schedule. 
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If  we  separate  applications  according  to  small,  medium  and  large 
size  pipe  (defined  dimensionally  as  1/8"  to  1/2",  1/2"  to  2",  and  2"  to  2m, 
respectively)  then  the  flow  cell  designs  that  will  commonly  be  selected 
are  those  depicted  in  Fig.  10.  It  will  be  recognized  that  the  flow  cell 
finally  selected  for  this  program  basically  resembles  the  "medium" 
style  of  Fig.  10,  with  the  "small"  style  recommended  for  higher  accuracy 
at  the  lowest  flow  rates  of  interest.  See  Table  3. 


Table  3.  Calculated  At  in  nanoseconds  in  small  and  medium 

size  conduits  at  T  =-50°C  (p  =  0.9  g/cm3,  c=  l600m/s) 
and  at  T  =  150°C  (p  =  0.  65  g/cm3,  c  =  750  m/s,  for 
L  -  4"  (10  cm). 


M,  PPH 

Small  Design 

1/8"  Schedule  80  Pipe.  ID=  0.  3" 

Medium  Design 

1"  x  1"  Square  Hole 

-  50°C 

+150°C 

-50°C  +150°C 

200 

48 

302 

3  19 

3000 

722 

4549 

51  321 

25000 

- 

- 

425  2678 

The  At  values  in  Table  3  are  calculated  for  the  "worst  case"  lowest 
temperature  condition  from  the  equation  At  =  2L V/c2  where  V  is  replaced  by 
M/pA.  That  is  to  say,  At  =  2LM/pAc^,  where  A  =  conduit  cross  sectional 
area.  At  the  high  temperature,  the  At  values  increase  over  the  low  tem¬ 
perature  values  by  a  factor  of  (pmax/pmin)(c  /Cmin)2  =  (0.  9/0.  65) 
(1600/750)2  =  6.  3. 

1.  Axial  Path,  Offset  and/or  Manifolded  Configurations 

This  category  includes  "medium  diameter"  designs  having  an  internal 
dimension  on  the  order  of  1"  (25.  4mm),  and  also,  "small  diameter"  designs, 
with  a  conduit  size  in  the  range  of  standard  1/8"  to  1/2"  pipe,  for  low  flow 
rates.  As  early  as  1954,  Kalmus  reported  on  a  method  for  interrogating 
effectively  along  a  substantially  axial  path  in  a  plastic  tube.  Practical  axial 
interrogations,  however,  awaited  later  methods,  i.  e.  ,  offset  designs. 

a.  Medium  Diameter,  Optionally  Including  Densitometer 
Tubes 

The  concept  of  interrogating  along  the  axis  of  an  offset  section  of 
conduit  has  been  considered  by  several  investigators  since  Katzenstein  and 
Katzenstein  ( 1 96 1 ).  Some  of  these  investigators  have  pointed  out  the  area- 
averaging  advantage  of  a  conduit  whose  full  cross  section  is  interrogated. 

In  May  1978,  N.  E.  Pedersen  suggested  using  an  axial  interrogation 
cell  in  which  transducers  would  be  installed  at  each  end  of  a  pipe  section, 
and  communicate  right  through  the  Simmonds  Precision  (SP)  densitometer. 
Pedersen  also  suggested  using  symmetrical  streamline  inlet  and  outlet 
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transitions,  and  possibly  additional  flow  control  surfaces.  Because  the 
flow  cell  design  had  to  be  frozen  before  this  concept  could  be  fully  investi¬ 
gated,  the  delivered  cell  design  resembles  the  Eustis  velocimeter  cell  of 
Fig.  4.  However,  we  tested  the  axial  path  right  through  the  SP  tubes,  shortly 
after  they  were  received  in  July,  1979.  Results  at  no  flow  were  encouraging: 
good  signal,  little  reverberation  (see  Fig.  11).  Accordingly,  more  work  with 
this  concept  may  prove  worthwhile,  whether  the  subsequent  densitometer 
is  based  on  dielectric  constant  or  torsional  wave  principles  (Lynnworth,  1977b, 
Fig.  13).  See  also.  Appendix  G. 

b.  Small  Diameter  for  Low  Flow  Rates 

For  a  constant  volumetric  flow  rate  Q,  At  is  inversely  proportional 
to  A.  Therefore,  the  use  of  small  diameter  tubes  to  increase  the  flow  velocity, 
and  hence  increase  At,  has  been  considered  by  several  flowmeter  manufacturers. 
Some  designs  were  reported  by  Lynnworth  (1977a  and  1979),  including  contri¬ 
butions  from  Bradshaw.  Examples  of  such  cells  are  given  in  Fig.  12.  (See 
also.  Fig.  10,  top,  and  Appendix  G.  ) 

The  cells  of  Fig.  10  (top)  and  some  of  those  in  Fig.  12  have  several 
features  apparently  unavailable  from  other  manufacturers  of  ultrasonic 
flowmeters: 

o  45°  inlet  and  outlet  fittings,  for  less  pressure  drop  than  90° 
tees  provide  (Bradshaw,  ca.  1977,  unpublished). 

o  Acoustically  isolated  transmitter  and  receiver  parts  in  which 
the  transducers  are  installed  (Lynnworth,  1975). 

o  Metallurgically- sealed  transducer  housings,  optionally  con¬ 
taining  an  acoustic  impedance  matching  member  between  the 
piezoelectric  transducer  and  the  stainless  steel  front  layer 
(Lynnworth,  Fowler  and  Patch,  1979,  unpublished). 

Limits  due  to  the  pressure  drop  at  high  V,  nonlinearity  due  to  end 
effects,  asymmetrical  flow  profile,  and  potential  cavitation  are  drawbacks 
to  this  design. 

2.  45°  Zigzag  Path,  Minimally  Obstructed 

This  design  is  based  on  the  area-averaging  principles  used  in  all  three 
Eustis  contracts.  See,  for  example,  Figs.  1,2,4  and  10  (middle).  Collimated 
beams  interrogate  the  entire  cross  section  of  fluid  passing  through  the  square 
holed  sleeve  within  the  cell.  The  square  holed  sleeve  also  makes  it  easy  to 
interrogate  along  an  axially-extended  in-line  low-pressure-drop  path  by  bouncing 
the  collimated  beams  in  zigzag  fashion,  achieving  L  =  4D  in  the  delivered  cell. 
Thus,  the  square  holed  sleeve  provides  100%  in-line  area  averaging  to  minimize 
profile  effects,  and,  in  contrast  to  a  round  pipe,  also  provides  for  distortion- 
free  multiple  reflections  to  increase  sensitivity,  since  At  is  proportional  to  L. 
The  deliverable  flow  cell  design  finally  selected  for  this  program  is  shown 
in  the  photograph  of  Fig,  13,  and  in  outline  form,  in  Fig.  14.  This  cell 
includes  optional  ports  for  a  torsional  wave  densitometer,  as  will  be 
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described  in  Sections  II D2  and  IUC.  In  simpler  form,  without  these 
optional  ports,  the  selected  flow  velocimeter  cell  for  this  program  is 
shown  in  outline  form  in  Fig.  15.  This  simplified  illustration  also  shows 
the  positions  of  the  square  holed  sleeve  and  transducers. 

The  delivered  cell  uses  the  longest  (6"  or  15  cm)  square-holed  one-piece 
stainless  steel  sleeve  that  was  commercially  available  without  involving 
special  tooling.  This  square-holed  sleeve  has  an  axial  hole  1"  x  1" 

(25.  4  mmx  25.  4  mm)  through  which  all  the  flowing  fuel  passes.  The  conically 
chamfered  inlet  and  outlet  desirably  reduce  turbulence  and  pressure  drop, 
but  undesirably  reduce  the  maximum  axial  distance  over  which  zigzag 
interaction  can  occur  to  L  =  4"  (10  cm).  Ports  for  transducers  are  screened 
to  reduce  eddy  formation  at  or  near  the  ports. 

The  only  practical  problem  encountered  so  far  with  this  design  is 
that  air  bubbles  are  occasionally  trapped  in  the  transducer  ports.  To  avoid 
this  potential  problem,  the  cell  should  be  installed  with  the  transducer  ports 
below  the  horizontal  plane  passing  through  the  sleeve  axis,  e.  g.  ,  transducer 
ports  on  bottom.  If  the  cell  axis  can  not  be  maintained  in  a  horizontal  or  nearly- 
horizontal  plane,  then  the  alternative  geometry  of  Fig.  16  may  be  necessary. 

If  a  longer  cell  were  desired,  so  as  to  increase  L,  one  could  elongate 
the  housing  and  install  two  (or  more)  relatively  standard  square  holed  sleeves, 
each  having  a  maximum  length  of  6"  (15  cm). 

3.  Bypass  and  Valved  Combinations  -  Concepts  for  Maximum  Rangeability 

The  flow  velocity  range  is  given  by  (M  /M  .  )  (p  /p  .  )  = 

7  max  min  max  man 

(25000/200)  (0.  9/0.  65)  =  173,  or  nearly  200: 1.  Initially  it  had  been  desired 
to  cover  M  down  to  100  PPH,  which  would  have  doubled  the  rangeability 
(turndown  ratio)  to  nearly  400:1. 

Two  concepts  have  arisen  which  might  be  useful  in  some  circumstances. 
M.  McDonald  suggested  using  a  bypass,  so  that  at  low  flow,  all  the  fuel 
would  pass  through  the  single  sensing  region.  At  high  flow  only  a  small  and 
presumably  known  fraction  would  be  sensed. 

A  second  and  somewhat  related  concept  is  to  use  two  flow  cells,  one 
small,  one  medium  (Fig.  10),  and  an  electronically  actuated  diverter  valve. 
Here,  all  the  fuel  is  to  be  sensed  either  in  one  cell  or  the  other.  For  example, 
the  small  cell  would  cover  M  up  to  perhaps  500  PPH,  and  the  medium  cell, 

M  from  500  to  25000  PPH.  Provision  would  have  to  be  made  to  divert  auto¬ 
matically  but  not  toggle  when  M  was  right  near  the  switch  point,  500  PPH  in 
this  example. 
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D.  Densitometer  Alternatives  Considered 

Seven  types  of  densitometers  were  considered.  The  advantages 
and  limitations  are  summarized  below,  including  reasons  for  rejecting 
all  but  die  last  type,  the  SP  capacitance  -  type. 

1.  Nuclear  Densitometer.  Measures  fuel  density  substantially 
independent  of  composition  or  temperature.  Rejected  because: 
radioactive  source  required;  high  activity  required  for  fast 
response;  development  cost  considered  prohibitive. 

2.  Temperature- Dependent  Densitometer.  Measures  fuel  tem¬ 
perature,  computes  "nominal"  or  average  density  for  a 
specified  fuel  type.  Rejected  because:  errors  could  exceed 
5%  due  to  density  variations  in  a  specified  fuel  even  at  room 
temperature. 

3.  Float- Type  Densitometer.  Might  be  usable  in  fuel  tank  or  in 
ground  testing.  Rejected  because:  does  not  measure  p  in  flow 
cell;  cannot  measure  p  transients;  sensitive  to  vibration  and 
and  acceleration,  especially  in-flight  conditions. 

4.  Resonant  Vane.  Response  derived  from  Newton's  law,  F  *  ma. 
Provides  a  resonant  frequency  output  signal  which  is  nonlinearly 
related  to  the  fuel  density  surrounding  the  probe.  Linearized 
electronically.  Manufactured  by  Agar,  ITT -Barton,  Solartron, 
others.  Rejected  because:  size  was  too  large;  it  required  a  2  inch 
pipe,  the  vane's  support  ring  alone  being  1.5  inch  diameter. 

5.  Reflectance-Type  Acoustic  Impedance  Probe.  Responds  to 
characteristic  acoustic  impedance  (pc)  of  fuel,  where  c  =  sound 
speed  in  fuel.  In  principle,  pc  could  be  combined  with  a 

" velocimeter"  measurement  of  V/c  (the  Mach  number).  This 
combination  could  provide  a  flush -mounted,  noninrastoe 
measurement  of  fuel  M.  Investigated  by  Panametrics  in  pre¬ 
vious  contracts  but  abandoned  in  favor  of  capacitance  type. 

[More  recent  designs  of  shear  wave  impedometer  probes,  such 
as  those  illustrated  in  Fig.  16  which  provide  a  number  of  bounces 
between  transmission  and  reception,  may  prove  useful  in  the 
future.  However,  1%  or  better  accuracy  has  not  yet  been  demon¬ 
strated.  (Lynnworth,  1979;  Lynnworth,  Seger  and  Bradshaw,  1980.] 
Rejected  because:  initial  designs  were  not  sufficiently  pc-aensitive 
or  stable, but  were  sensitive  to  temperature  effects  ih  probe;  ac¬ 
curacy  may  degrade  if  contaminant  builds  up  on  surface. 
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6.  Slow  Torsional  Wave  Densitometer .  New  device  based  on 
discovery  in  1976  by  one  of  the  principal  investigators  that 
the  speed  of  a  torsional  wave  at  a  frequency  of  ~  0.  1  MHz 
that  is  guided  in  a  noncircular  waveguide,  is  inversely  pro¬ 
portional  to  the  density  of  the  surrounding  liquid.  Requires 
only  2  electrical  conductors,  which  could  be  the  same  ones 
used  for  the  V  measurements.  Under  development  for  nuclear 
application  in  two-phase  water/ steam  mixture  (Lynnworth,  1977b; 
Arave,  1978,  1979.)  [See  Recommendation  IV  A.  j  Rejected 
because:  accuracy  not  yet  demonstrated  to  be  good  enough;  probe 
requires  temperature  compensation;  mounting  may  be  difficult. 

7.  Capacitance-  Type  Densitometer.  Based  on  Clausius*  Mob otti 
relationship  between  fuel  density  and  dielectric  constant  (Stuart, 
1974).  Hardware  manufactured  by  several  vendors.  Simmonds 
Precision  (SP)  provides  the  best,  to  our  knowledge,  and  their 
equipment  has  been  calibrated  and  used  extensively  in  aircraft 
and  in  JP-4,  JP-5,  etc.  The  SP  equipment  was  used  in  our 
previous  Eustis  program  and  performed  satisfactorily.  It  was 
therefore  selected  by  us  for  use  in  our  present  program,  despite 
the  size  and  intrusiveness  of  the  coaxial  tubes  and  supports 
required  for  responding  to  the  fuel's  dielectric  constant,  and 
potential  problems  with  water  in  the  fuel,  and  some  shift  in  cali¬ 
bration  with  different  fuel  types,  and  high  temperature  nonlinearity. 
This  densitometer  was  selected  because:  it  appeared  to  meet  all 
accuracy,  rangeability,  fuel  types  and  other  requirements,  and 
was  compatible  with  the  rest  of  our  M  system. 


I 

i 


NADC-80254-60 


III.  MASS  FLOWMETER  FOR  PRESENT  PROGRAM 

A.  Selected  System 

Figure  18  shows  in  block  diagram  form  the  selected  M  system.  It 
consists  of  a  "spoolpiece"  containing  the  SP  densitometer  at  the  inlet  or  outlet 
end,  the  Panametrics  ultrasonic  flow  velocimeter,  and  optionally,  quick-con¬ 
nect  valved  fittings.  The  electronics  consists  of  a  modified  model  6000  flow¬ 
meter  instrument  connected  to  a  SP  densitometer  electronic  module,  the 
system  providing  a  square  wave  output  proportional  to  M.  If  it  is  desired  to 
operate  the  V  part  with  clamp-on  transducers,  or  other  cells,  this  can  be 
accommodated  in  the  software  (Appendix  G). 

B.  Description  of  Electronics  for  Flow  Velocity  and  Mass  Flow  Rate 
Determination 

The  ultrasonic  flow  velocimeter  electronics  consists  essentially  of  two 
basic  parts:  the  phase  locked  loop  controlled  signal  section  (Section  1)  and 
the  digital  computer  section  (Section  2).  The  first  section  generates,  trans¬ 
mits  and  receives  the  ultrasonic  signals.  It  also  produces  outputs  that  are 
fed  to  the  second  section.  The  computer  section  accepts  the  output  of  the 
first  section  and  calculates  the  flow  velocity.  The  computer  also  runs  any 
displays  or  outputs,  and  can  be  optionally  formated  in  different  units  and 
scales.  See  Fig,  8.  (Appendix  C  contains  a  more  detailed  description.  ) 

The  first  section  consists  of  a  modified  phase  locked  loop  (PLL). 

The  output  of  the  voltage  controlled  oscillator  (VCO)  is  used  to  generate  15V 
transmit  pulses  that  are  sent  to  the  upstream  or  downstream  transducer, 
depending  upon  the  state  of  the  configuration  switch  (CS).  The  CS  determines 
which  transducer  acts  as  transmitter  and  which  acts  as  receiver.  The  CS  changes 
state  (from  upstream  to  downstream  and  vice  versa),  when  a  fixed  number 
of  transmit  pulses  have  occurred.  This  way  both  upstream  and  downstream 
configurations  always  involve  the  same  number  of  transmitted  pulses.  After 
propagating  over  the  liquid  fuel  in  the  flow  cell,  the  received  pulses  are 
detected  and  amplified  by  an  automatic  gain  controlled  (AGC)  amplifier.  The 
amplified  received  pulse  is  detected  by  a  zero  crossing  detector  (ZCD).  The 
output  of  the  ZCD  is  fed  to  the  phase  detector  along  with  the  original  VCO 
output.  The  50%  duty  cycle  point  of  the  VCO' s  output  and  the  output  of  the 
ZCD  (receive  signal)  are  compared  (phase  detected)  and  any  difference  in 
time  between  them  is  used  to  correct  the  VCO' s  frequency  so  that  the  output 
of  the  phase  detector  is  maintained  as  close  as  practical  to  zero.  This 
ensures  that  the  time  between  transmit  pulses  is  essentially  twice  the  transit 
time  of  the  ultrasonic  pulse  from  transmitter  to  receiver  transducer.  (Typical 
received  signals  obtained  using  JP-5  at  various  flow  rates  are  shown  in  Fig.  22.  ) 
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Since  the  time  of  transit  of  the  ultrasonic  pulse  from  transmitter 
to  receiver  transducer  is  increased  when  the  ultrasonic  wave  propagates 
against  the  flow,  but  decreased  when  propagating  with  the  flow,  the  output 
frequency  of  the  VCO  will  be  changed  accordingly  as  the  PLL  tracks  the 
changing  transit  times  t^  and  t^.  (Various  pulse  widths,  blanks,  filters 
and  thresholds  are  selected  to  optimize  system  performance  and  minimize 
interference  from  spurious  pulses  or  noise.  )  This  means  that  when  the 
CS  is  in  one  configuration,  say  upstream,  the  frequency  will  be  less  than 
the  no  flow  condition.  When  this  CS  is  in  the  downstream  configuration, 
the  VCO  frequency  will  be  greater  than  the  no  flow  frequency.  Since  both 
upstream  and  downstream  configurations  use  exactly  the  same  number  of 
pulses,  N,  the  total  time  that  the  CS  is  in  either  configuration  will  proportionately 
increase  or  decrease  in  response  to  the  flow  velocity.  As  the  flow  increases 
the  total  upstream  time  2Nt^  =  T  increases  and  the  total  downstream  time 
2Nt^  =  T?  decreases.  The  output  of  the  CS  is  fed  to  the  second  section 
where  each  configuration  time  T^  and  T^  (number  of  pulses  times  period 
of  VCO  output)  is  measured. 

The  computer  section  consists  of  two  24  bit  counters,  a  10  MHz 
oscillator  and  some  arbitration  logic.  The  arbitration  logic  receives  the 
input  times  from  the  first  section  (the  output  of  the  CS)  and  switches  the 
10  MHz  clock  back  and  forth  between  the  two  24  bit  counters  for  the  duration 
of  each  configuration  time.  The  arbitration  logic  also  tells  the  computer 
which  time  (upstream  or  downstream)  it  is  measuring  and  when  the  mea¬ 
surement  is  done.  The  computer  then  takes  the  two  times  Tj  and  (total 
upstream,  total  downstream)  and  performs  calculations  on  them  based 
upon  the  number  of  pulses  in  each  direction,  the  physical  dimensions  of 
the  flow  cell  and  any  scale  factors  necessary  (see  Appendix  2).  The  cal¬ 
culation  yields  the  flow  velocity  V.  The  computer  can  then  use  the  flow 
velocity  to  derive  the  volumetric  or  mass  flow  rate  Q  or  M  respectively 
using  additional  information  about  the  cell  or  the  fluid.  The  computer  then 
outputs  the  result  in  a  given  format  and  resets  the  system  for  the  next 
measurement.  In  the  present  program  the  desired  format  is  a  square  wave 
whose  frequency  is  proportional  to  the  mass  flow  rate. 

The  response  time  of  the  digital  portion  of  the  electronic  system, 
if  defined  as  the  maximum  time  required  to  fully  respond  to  a  step  change 
in  flow  rate,  is  essentially  2(T  +T^).  For  example,  if  the  average  transit 
time  were  250ps,  and  N  =  1024,  the  response  time  would  be  4x250xl0^x 
1024  -  1.024  second,  or  approximately  Is. 

The  mass  flowmeter  electronics  multiplies  p  and  V  as  follows.  The 
flow  velocity  V  is  available  in  the  computer  from  the  previous  calculation. 

The  computer  is  also  provided  with  A,  the  cross  sectional  area  of  the  flow 
cell  so  that  the  computer  multiplies  V  times  A  to  obtain  Q,  the  volumetric 
flow  rate,  in  gallons/sec.  The  density  pis  obtained  from  the  SP  densito¬ 
meter  electronics.  The  output  of  the  SP  electronics  is  a  voltage  that  is 
related  to  density  and  that  voltage  is  used  to  calculate  actual  density  in 
lb/gal  by  the  following  relation: 

p=  2.505  Ib/gal  +  (19.05  Ib/gal  •  volt)(R  volts) 
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where  R  =  V  /V  and  VQut  is  the  output  voltage  of  the  SP  electronics 
and  V  is  lSe*  supply  voltag'e  for  the  SP  electronics  (15V).  The  output 
voltage6^  is  fed  to  a  12  bit  analog  to  digital  converter  that  is  controlled 
by  the  computer.  The  program  in  the  computer  tells  the  A/D  to  make  a 
conversion  and  upon  completion  of  the  conversion  the  computer  uses  the 
resulting  digital  code  to  calculate  V  .  The  computer  solves  the  above 
equation  for  p  and  then  multiplies  the  calculated  p  times  the  calculated  Q 
to  yield  the  mass  flow  rate  Ki  in  lb  per  second.  This  number  is  then  scaled 
and  used  as  a  divisor  for  a  computer  using  a  1  MHz  clock  so  that  the  resultant 
square  wave  output  of  the  counter  is  the  frequency  output  proportional  to 
mass  flow  rate.  This  square  wave  output  is  thus  the  required  measure  of  M. 


C.  Description  of  Flow  Velocimeter  Cell  and  Transducers 


The  flow  velocimeter  (V)  cell 
introductory  manner  in  Section  II CZ. 
are  as  follows:  gp 

Densitometer 


has  been  partially  described  in  an 
Size,  weight  and  other  characteristics 

Panametrics 

Cell  Flow  Velocimeter  Cell 


Length 


5"  ( 1 27  mm) 


10-7/8"  (276  mm) 


Diameter  and/or  2"  diam  -f  elec, 

other  dimensions  connectors 


2.  5"  diam  x  5"  high 


Inlet 

Outlet 

Weight 

Wetted  metals 
Drwg.  No. 


1-1/2"  NPT  male 
1-1/2"  NPT  male 
2.  21b  (1  kg) 

SS  304 

SP  #  Rev 473279 


1-1/2"  NPT  female 
1-1/4"  NPT  female 
10  lb  (4.5kg) 

SS  304 

Panametrics  138-401 


The  V  cell  is  functionally  equivalent  in  principle  to  the  design  shown 
in  Fig.  10,  middle.  However,  to  reduce  the  weight  of  standard  laterals, 
the  cell  body  was  welded  using  SS  304  pipe  sections,  nominally  1.  5"  schedule 
80.  A  1.  5"  3000  #  half  coupling,  SS  304,  was  welded  to  its  inlet  to  adapt  to 
the  SP  p  cell.  The  outlet  is  tapped  1-1/4"  NPT,  into  which  a  nipple  and 
Snap-Tite  quick-connect  valved  fitting  may  be  installed. 
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The  axial  portion  of  pipe,  in  the  simplest  case  (Fig.  15)  would  be 
tapped  for  1/4-20  set  screws  which  hold  the  square  holed  sleeve  in  place. 
For  permanent  installation  the  set  screws  could  be  welded  in  place. 


The  square  holed  sleeve  has  two  transducer  ports  milled  on  one 
side,  collimating  the  interrogating  beam  to  provide  the  desired  square  or 
rectangular  beam  envelope  which  then  zigzags  along  the  square  holed  sleeve. 
The  collimating  ports  are  covered  with  40  mesh  SS  304  wire  cloth  (0.  005", 
or  ~  0.  1  mm  wire  diameter),  spot-welded  in  place. 

We  had  also  considered  using  a  membrane  of  metal  or  other  material 
instead  of  the  wire  cloth,  to  further  reduce  the  tendency  of  the  transducer 
ports  to  generate  eddies.  Figure  19  is  a  plot  of  the  calculated  attenuation 
vs  thickness  for  several  candidate  membranes.  However,  it  was  decided 
to  use  screens,  as  they  had  been  used  in  previous  Eustis  contracts  with¬ 
out  any  problem  other  than  bubble  retention.  It  was  thought  that  this 
problem  could  be  aggravated  if  a  membrane  tended  to  seal  the  collimating 
ports.  The  screen  has  the  further  advantage  of  being  transparent,  so  if  the 
transducers  are  removed,  visual  inspection  of  part  of  the  square  hole  channel 
is  possible,  without  removing  the  entire  V  cell  from  a  flow  loop. 

The  initial  transducer  design  used  1.25"  (31  mm)  diameter  2  MHz 
discs  made  of  a  lead  zirconate  titanate  composition  with  a  curie  point  above 
300°C.  They  were  sandwiched  between  a  thin  (<  0.1mm)  corrosion- resistant 
metal  shim  and  an  insulated  graphite  backing.  The  backing  was  rigid  enough 
to  withstand  the  1600  psig  design  pressure,  yet  attenuating  enough  to  not 
limit  the  pulse  repetition  frequency.  The  transducer /backing  portion  was 
contained  in  an  all-metal-sealed  container  or  well  which  was  electron  beam 
welded  to  the  end  of  a  conventional  1.25"  threaded  SS  304  hex  head  plug.  The 
plug  was  modified  to  accommodate  a  Microdot  electrical  connector  and  a  gas 
tight  Swagelok  fitting  and  3/8"  (~4  mm)  tube  which  seals  and  protects  the 
coaxial  cable. 

The  hex  head  pipe  plugs  screwed  into  scooped  couplings  which  had  been 
tungsten  inert  gas  (TIG)  welded  to  the  cell  body.  When  the  flowmeter  is  in 
use,  the  transducers  should  be  oriented  so  that  air  does  not  get  trapped  in  the 
transducer  ports,  i.  e.  ,  transducers  should  be  downward  at  least  until  all 
gas  is  removed  from  the  liquid  under  test. 

For  extended-well  thin-window  flow  cell  geometric  parameters  were 
as  follows: 


fuel  path  length  P  =  7.5  (~  525  mm)  nominal 
axial  path  length  L  -  4.020"  (~  100  mm) 

duct  area  A  =  1.010  in.  ^  (~  6.  3  cm^) 
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Transducer  and  cable  delays  were  approximately  t  =  40  ns  (shims)  + 
20ns  (2  cables)  =  60 ns.  If  a  thick  window  or  buffer  rod  were  placed  between 
the  piezoelectric  element  and  the  fuel,  t  would  increase  by  2X/c^  where 
X  =  window  or  buffer  dimension  and  c ^  =  speed  of  sound  therein. 

In  the  delivered  system,  2  MHz,  25  mm  diameter  lead  metaniobate 
transducers  were  installed  in  SS  3 1 6  hex  head  1.25"  pipe  plugs.  These  were 
bored  to  provide  windows  of  thickness  X  =  30  mm.  Transducers  were  in¬ 
stalled  in  the  delivered  flow  cell,  the  appropriate  constants  being: 

P  =  10.8  in.  (274.  3  mm) 

L  =  4.020  in.  (102.  1  mm) 

A  =  1.010  in.2  (6.  581  cm2) 

t  =  10  us 
w 

These  values  were  entered  in  the  delivered  "constants  prom." 

Compared  to  the  initial  thin  window  design,  the  thick  window  SS  plug 
design  sacrifices  a  few  dB  of  signal  (due  to  impedance  mismatch)  in  exchange 
for  simpler,  more  rugged  construction  and  ease  of  self- testing.  Regarding 
the  latter  advantage,  each  thick  window  transducer  can  be  self-tested  in  the 
pulse-echo  mode  as  a  function  of  temperature,  e.  g.  ,  in  liquid  nitrogen  at 
-196°C,  boiling  water  at  +  100°C,  oven  at  +150°C,  etc.  The  thick  window 
transducer  is  more  resistant  to  corrosion,  impact,  abrasion  and  pressure, 
and  its  wetted  face  can  be  contoured  to  reduce  the  amplitude  of  triple  transits 
within  the  fuel  path  to  shorten  the  response  time  or  to  allow  more  interrogations 
within  a  given  response  time. 
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Some  additional  remarks  on  the  square  holed  sleeve  follow.  In 
Fig.  l6,  if  the  recessed  reflecting  face  is  to  extend  the  full  width  of  the  square 
hole  (i.  e.  ,  1  inch  or  25.4  mm)  then  the  minimum  diameter  of  the  square  holed 
sleeve,  Dg,  may  be  calculated  as  follows.  Referring  to  the  accompanying 
sketch  wherein  the  axial  projection  of  the  reflecting  face  is  denoted  w,  its 
maximum  depth  d  is  given  by  d  =  w tan ?.  (For  a  45°  zigzag  path,  9  =  22.5°.  ) 


=  ?  J  (f)  +  (f  +  w  tan  °) 


Now  if  w  =  e/2  this  simplifies  to 


D  = 
s 


7 


For  a  square  holed  sleeve  having  2w  =  e  =  1  inch  and  0  -  22.  5°  =  tt  /8  radians, 

=  J 0.5+  (1  +0.4142)2  =  .JO. 5+  2  =  1.58  %  it  >2  inches 


40.  1  mm 


NADC -80254-60 


D.  Description  of  Densitometer 

This  discussion  is  adapted  from  Stuart's  1974  paper  in  which  the  SP 
densitometer  is  described.  The  well-established  relationship  between  density 
and  dielectric  constant  for  aviation  fuels  has  long  been  used  in  aircraft  fuel 
gaging  systems. 

The  hydrocarbons  in  aircraft  fuels  are  non-polar  liquids.  For  a  non¬ 
polar  liquid,  a  theoretical  relationship  between  dielectric  constant  and  density 
is  the  Clausius -Mosotti  law 

(e  -  l)/(e+2)  =  Cp 

where  e  =  dielectric  constant,  p  =  density,  and  C  =  polarizability.  This 
relationship  is  rearranged  and  called  the  "fuel  regression  line"  in  capacitance 
fuel  gaging  literature.  One  fuel  regression  line  is: 

(c  -  l)/p  =b  +  a(e-l) 

where  b  =  0.  12192,  a  =  0.  04112,  and  p  in  lb/gal.  The  polarizability  varies 
only  slightly  from  one  hydrocarbon  to  another  hydrocarbon  for  the  principal 
constituents  of  aircraft  fuels.  Thus,  the  dielectric  constant,  as  expressed 
by  the  Clausius -Mosotti  law,  is  a  good  measure  of  fuel  density.  Stuart's 
paper  contains  a  plot  of  density  vs  dielectric  constant  for  some  typical  samples 
of  JP-4.  Stuart  illustrates  what  is  ccmmonly  called  "fuel  error."  This 
fuel  error  is  the  deviation  of  a  given  sample  from  the  theoretical  best  aver¬ 
age  straight  line  which  is  called  the  "JP-4  line  of  regression."  Further 
statistical  analysis  of  this  and  also  JP-5  data  is  *lso  given  in  Stuart1  s  paper. 
The  interpretation  of  this  plot  is  that  95.  5%  of  the  samples  measured  (JP-4 
and  JP-5)  fall  within  1.  3%  of  the  JP-4  line  of  regression.  From  these  data 
it  is  thus  possible  to  predict  a  2a  error  of  1.  3%  or  less.  The  data  has  been 
further  analyzed  and  la  and  3a  errors  are  included. 

To  more  certainly  establish  the  constancy  of  this  relationship  for  nor¬ 
mal  commercial  fuels,  more  recent  data  was  collected  by  Air  Canada  and 
United  Air  Lines.  This  data  show  that  of  55  samples,  51  are  within  +_  1%  of 
variation  from  the  JP-4  line.  This  tends  to  confirm  the  previous  prediction. 

It  is  known  that  the  Clausius-Mosotti  law  is  nonlinear,  statistically 
corresponding  to  0.5%  for  fuels  at  +70°C  and  -45°C  (range,  115°C.).  For 
the  range  -38  C  to  +150°C  the  corresponding  nonlinearity  needs  to  be  de¬ 
termined. 

The  particular  p  cell  fabricated  for  this  program  was  tested  by  SP 
in  JP-5,  varying  the  temperature  to  change  fuel  density.  According  to  SP' s 
Owen  Clay  (priv.  comm.  ,  July  1979)  the  data  were  used  to  establish  an 
equation  of  the  form 

P  =  2-  505g-ft+ll9-05^Kv]  xtRl 

where  R  (volts)  is  obtained  by  a  ratio  meter  to  read  the  output  of  the  SP 
signal  conditioner  in  reference  to  a  +  15V  power  supply:  R  =  output 
voltage  f  15.  24 
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IV.  TEST  PROCEDURES  AND  TEST  RESULTS 


A.  Accuracy 

The  principal  investigators  judge  the  flow  velocity  testa  to  be 
the  most  important  from  among  the  several  different  tests  specified  in 
the  contract.  The  reason  for  this  judgement  is  that  the  densitometer  has 
been  proven  elsewhere;  hydrostatic  tests  are  easily  satisfied  if  adequate 
thicknesses  of  flaw-free  engineering  materials  are  used;  transducers  and 
coupling  methods  are  now  available  which  exceed  the  temperature  objectives 
of  the  contract;  noise  due  to  vibration  should  be  filterable,  given  the  2  s 
response  time.  Accordingly,  except  for  electronic  tests  conducted  throueh- 
out  the  program,  more  effort  was  placed  on  flow  velocimeter  tests  than 
on  any  other  test.  Nevertheless,  due  to  funding  limits  and  piggy-back 
schedule  constraints,  we  were  unable  to  conduct  as  many  flow  tests  as 
desired;  more  flow  testing  is  therefore  recommended,  especially  for  tem¬ 
perature  and/or  viscosity  ranges  not  experienced  in  the  present  contract. 

Because  of  the  inconvenience  and  high  cost  of  conducting  fuel  flow 
calibration  tests  throughout  the  development  of  this  contract's  flowmeter, 
we  used  water  for  most  of  the  flow  velocimeter  testB.  The  complete  M 
system,  however,  was  tested  at  GE-Lynn  on  several  occasions  using  JP-5 
and,  finally,  the  fuel  substitute  7024  Type  II  (Stoddard  solvent).  In  these 
tests,  the  reference  M  was  based  on  turbine  readings  for  V  and  p 
which  was  computed  from  the  temperature  of  the  known  fuel  or  fuel  substitute. 
Density  had  been  checked  occasionally  with  a  hydrometer,  so  the  reference 
M  determination  was  considered  to  be  in  error  by  less  than  1/2%  of  reading, 
from  ~  1000  to  ~  25000  PPH  (~500  to  ~  10,  000  kg/hr). 

The  flow  cell  and  flow  velocimeter  electronics  unfortunately  did 
not  function  properly  on  the  final  occasion  when  the  fuel  substitute  cali¬ 
bration  test  was  conducted  at  GE.  A  circuit  deficiency,  discovered  after 
that  test,  handicapped  the  receiver  gain  by  18  dB.  An  additional  transducer 
handicap  of  12  dB  occurred,  apparently  due  to  a  partial  delamination. 
Therefore,  the  complete  and  final  delivered  system  was  tested  on  the 
fuel  substitute  only  up  to  16000  PPH.  After  this  final  fuel  substitute 
test,  however,  the  circuit  was  repaired,  the  transducers  were  replaced, 
and  the  flow  velocimeter  was  demonstrated  at  Panametrics  in  water  to  over 
25000 PPH,  the  contract  M  objective.  (Fuel  substitute  data  are  given  in  Table  4.) 

Weigh  tank  water  calibration  tests  were  later  conducted  to  38 ft/s  (11  6m/s) 
or  60,000  PPH  (Table  6,  Fig.  30). 

To  determine  the  10-day  repeatability  of  the  equipment,  a  closed- 
loop  system  of  water  flow  using  gravity  feed  was  utilized.  Although  water 
has  different  density  and  viscosity  than  fuel,  its  fluid  dynamics  are  adequate 
with  respect  to  the  goal  of  this  test.  Because  of  the  low  flow  rates  achievable 
with  the  gravity  feed  tester,  a  different  flow  cell  (offset  type)  was  used  in 
this  repeatability  test. 
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Following  Bernoulli' s  theory  the  velocity  of  fluid  flow  in  a  piping 
system  is  proportional  to  the  square  root  of  the  gravitational  difference 
in  head,  and  inversely  proportional  to  the  square  root  of  the  sum  of  the  friction 
factors  due  to  friction  loss  in  pipe  and  pipe  fittings  between  any  two  points. 

For  the  system  described  below  the  velocity  is  represented  by  the  equation 


Vcell<ft/sec)  = 


2  g  -  AHmanometer  ) 

0.586  +  93. 17  f  ,  +  53~4f 

tube  ss 


1  1/2 


(4-1) 


where  AH  (feet)  =  difference  of  height  between  inlet  and  outlet  of  the 
tubing  connecting  the  flow  cell 

AH  (feet)  =  pressure  loss  in  flow  cell 

manometer 

f  a  factor  for  friction  loss  in  flexible  tubing  (for  ~  =  0.000001) 
tube  a 

f  =  factor  for  friction  loss  in  stainless  steel  tubing  (for  ~  =  0.003) 
ss  d 

g  =  32.158-^-y 
sec^ 

Figure  26  shows  details  of  the  system.  The  recirculating  pump 
kept  water  overflowing  from  the  top  container.  Constant  head  was  therefore 
maintained  at  the  inlet  of  the  tubing  leading  to  the  flow  cell.  By  siphoning 
action,  water  left  the  top  container  and  followed  the  inlet  tubing  to  the  flow 
cell.  After  traveling  through  the  system,  the  fluid  returned  to  the  big  con¬ 
tainer  in  which  its  temperature  was  monitored  and  kept  constant  by  a  tem¬ 
perature  controller.  The  importance  of  controlling  both  the  room  and  fluid 
temperature  is  to  maintain  the  viscosity  of  water  at  one  particular  value 
during  the  duration  of  the  test  (the  factors  for  friction  losses  in  the  flexible 
tubing  and  the  flow  cell  depend  on  the  viscosity  of  water)  so  that  the  repeat¬ 
ability  of  the  flow  rate  of  water  can  be  obtained  from  day  to  day  at  each 
different  height  of  head  (AH).  The  applicability  of  Eq.(4-1)  to  the  test 
arrangement  of  Fig.  26  was  verified  by  timing  the  interval  required  to  collect 
100  mi  samples  of  water.  Test  data  are  plotted  in  Fig.  28  and  reported  in 
Table  5. 

The  inside  diameter  of  the  offset  flow  cell  was  5/l6"  (79  mm)  and  had 
the  following  characteristics: 

L  =  5-3/16"  (132mm) 

P  =  6-5/8"  (l68mm) 
t  =  1  ps 
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The  measured  velocity  is  determined  by 


/Volume  of\ 

...  .  \  container  '  1  1  ft 

measured'  /Time  to  fill)  /Cross-sectional  area\  12  in. 

'  container  '  '  of  flow  cell  ' 


(4-2) 


6 1.02  in. 

where  volume  of  container  =  1000  mf  x  — : — ^ - ;  time  to  fill  container  = 

10  m  f 

elapsed  time  (seconds)  for  water  to  fill  the  1000  m|  container;  cross  sectional 

..  it  (.3125  in.  )2 

area  of  flow  cell  =  — 1 - ■ - 

4 

The  water  temperature  was  kept  at  a  temperature  within  1°C  of  room 
temperature.  Room  temperature  was  22°C  at  each  head  setting.  The  mea¬ 
sured  velocity  in  Table  5  is  the  average  velocity  measured  over  10  days;  and 
the  measured  velocity  for  each  day  is  the  average  value  as  6  measurements. 
Each  measured  velocity  has  an  overall  uncertainty  of  +_  0.050  ft/sec 
(15.239  mm/sec)  based  on  a  standard  error  of  0.008  ft/sec  (2.438  mm/sec) 
and  an  allowance  of  0.026  ft/sec  (7.925)  mm/sec  for  systematic  error 
(Eisenhart,  1969). 


From  the  data  in  Table  5  the  velocimeter  portion  of  the  mass  flow¬ 
meter  is  found  to  have  the  following  specifications: 


Repeatability4  +  0.010  ft/sec  (+  3.048  mm/sec)  Using 
b  “  “  Offset 

Linearity  +  0.008  ft/sec  (+  2.438  mm/sec)  Cell 


Later,  using  the  delivered  cell  in  the  Foxboro  water  calibration 
tests,  where  V  =  38  ft/ s  (11.6  m/s),  the  following  results  were  ob¬ 
tained  (Table  6^aX 

Zero  flow  repeatability  +0.  01  ft/s  (offset  =0.  06  ft/s) 

Repeatability  at  full  scale  0.  06%  of  reading 


a  Repeatability  = 
(Abernethy,  1980). 


1/ One  standard  deviation  \  /One  standard  deviation\ 

(  of  10  days  electronics  ]  +  (  of  10  days  measured  ) 

reading  /  '  velocity  ' 


b  Linearity  =  Largest  difference  between  electronic  reading  and  best  fit 
regression  curve  at  a  measured  velocity  reference. 
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Table 

4  .  Calibration  test  using  Stoddard  solvent. 

Data  listed 

in  order  observed: 

first  run  down,  first  run  up,  second 

run 

down. 

Most  M 

values  are  averages 

of  2  readings; 

values. 

average  of  2  to  4  readings. 

Mc 

KPPH 

K&H 

T 

°F 

P 

PSIG 

Mp 

Mc 

16.37 

17.20 

77.2 

170+15 

1.051 

12.01 

12.85 

78.8 

100  +  10 

1.070 

10.00 

10.65 

78.8 

75  +  10 

1.065 

7.06 

7.45 

79.0 

35  +  10 

1.055 

5.00 

5.30 

79.4 

25  +  5 

1.060 

3.04 

3.30 

81.0 

20+3 

1.086 

2.00 

2.23 

83.5 

10  +  2 

1.115 

1.00 

1.10 

91.1 

10  +  2 

1.100 

2.505 

2.93 

83.2 

11  +  1 

1.170 

3.988 

4.33 

80.9 

23  +  3 

1.086 

5.990 

6.43 

80.1 

34+10 

1.073 

7.956 

8.37 

80.1 

50+  15 

1.052 

10.95 

11.75 

80.2 

- 

1.073 

13.85 

14.87 

80.5 

120  +  25 

1.074 

15.95 

16.95 

81.3 

160  +  20 

1.063 

14.95 

16.10 

81.3 

140  +  20 

L.  077 

12.90 

13.85 

81.8 

115+15 

1.074 

11.50 

12.37 

81.8 

85+  15 

1.076 

9.00 

9.59 

81.7 

60+12 

1.066 

5.51 

5.98 

82.1 

30+10 

1.085 

3.50 

3.78 

83.3 

20  +  2 

1.080 

1.27 

1.51 

87.4 

13+3 

1.189 

1.  00 

.  . 

1.22 

90.1 

10+1 

1.220 

M^  =  Reference  value  for  M  measured  with  Cox  system. 
=  Measured  with  Panametrics  system. 

P  =  Pressure  of  test  fluid  adjacent  the  flow  cell. 


Date  of  test:  March  1980 


faSSufeuMb* 
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Water  calibration  tests  were  conducted  at  the  Foxboro  Com¬ 
pany,  Foxboro,  MA,  on  23  and  24  September  1980.  Here,  the  delivered 
cell  and  the  V  part  of  the  electronics  were  calibrated  between  0. 1  and 
38  ft/s  using  a  weigh  tank.  The  fluid  was  water  between  80  and  90  F. 

Data  were  obtained  starting  at  the  maximum  velocity  available  for 
this  test,  38  ft/s  (  11.6  m/s).  At  this  velocity,  the  K  switches  in  the  modi¬ 

fied  model  6000  were  set  to  0.  9999  and  the  ultrasonic  V  essentially 
equaled  the  actual  flo\^  velocity.  [Theoretically,  K  =  1.000.  The  maxi¬ 
mum  Re  was  ~  3  x  10  .  The  maximum  mas£  flow  rate  was  Mmax  =  pVA  = 

(62.  4  lb/ft3)  (38  ft/s)  (1.  005  in.  /144  in.  /ft  )  =  16.  49  lb/s  a  60,  000  lb/hr 
(7.  5  kg/s  or  27,000  kg/hr).  However,  since  the  delivered^SP  densitometer 
is  designed  for  fuels,  not  water,  it  and  the  corresponding  M  computation 
were  not  used;  instead,  only  V  was  computed  and  displayed  in  this  calibra¬ 
tion  test.]  Other  than  the  K  switch  setting,  the  only  further  electronic 
adjustment  made  during  the  test  was  the  response  time,  which  was  switch- 
selected  in  the  model  6000  by  Foxboro  personnel  to  average  the  flow  readings 
over  intervals  approximately  proportional  to  the  time  required  to  collect 
a  pre- determined  mass  of  water.  Thus,  at  lower  flow  rates,  longer  re¬ 
sponse  times  were  used  to  average  128  readings,  256  readings,  etc. 

To  conduct  the  test,  V  was  decreased  from  38  ft/s  down  to  0.  1  ft/s, 
a  range  of  380:  1,  thereby  encompassing  turbulent,  transitional  (~  1/2  to 
1/4  ft/s),  and  laminar  flow.  A  few  data  points  were  also  obtained  on  in¬ 
creasing  flow  back  up  to  20  ft/s.  The  data,  given  in  Table  6  and  Fig.  30,  show 
that  V  zero  repeatability,  accuracy  and  linearity  of  better  than  1/2%  full 
scale  were  achieved. 

These  water  flow  calibration  or  scale  error  data  may  be  considered 
as  an  extension  of  the  Stoddard  solvent  tests  at  GE  (Table  4)  where  the 
maximum  M  was  16,  370  lb/hr.  In  other  words,  the  water  flow  data  in 
Table  6  and  Fig.  30  show  that  the  V  part  of  the  system  operates  properly 
at  velocities  which  correspond  to  fuel  mass  flow  rates  at  least  up  to  50% 
beyond  the  contract-specified  value  of  25,000  lb/hr.  Later,  on  31  March 
and  1  April  1981,  it  was  demonstrated  that  the  zigzag  cell  and  M  electronics 
could  operate  in  a  fuel  substitute  (7024  BII)  at  flow  rates  up  to  at  least 
31,450  lb/hr.  Data  in  Appendix  K  show  that  operation  to  the  limit  of  the 
test  stand  under  the  test  conditions  was  achieved  without  signal  degradation. 
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Table  6.  Water  calibration  test  at  the  Foxboro  Company. 
Date  of  tests:  23  and  24  September  1980. 


Flow  Rate 

Temp. 

v 

21  ct\l&l 

Error 

ft 3  /  s 

°F 

ft/ 0 

ultra 

ft/ B 

%  Reading 

%  Full  Scale 

.  2668 

86.4 

38.  045 

37.  90 

-.  38 

-.  38 

..  2672 

87.  1 

38.  096 

37.  93 

-.  44 

-.  44 

.  2109 

87.4 

30. 075 

30.  00 

-.25 

-.20 

.1797 

87.  7 

25.  620 

25.  58 

-.  16 

-.  11 

| 

.  1412 

88.  2 

20.  128 

20.  20 

+.  36 

+.  19  ' 

.  1047 

88.  6 

14.  923 

15.  00 

-h  52 

+.  20  ! 

.07136 

89.  2 

10.  173 

10.  20 

+.  27 

+.  07 

.03517 

90.0 

5.  014 

5.  071 

+  1.  14 

+.  15 

. 007844 

90.  2 

1.  118 

1.  094 

-2.  15 

-.  06 

. 004008 

89.  6 

.  5715 

.  540 

-5.  51 

-.  08 

. 001916 

88.  8 

..  2732 

.  217 

-20.  6 

-.  15 

. 00000 

.  0000 

-0.  07 

-.  002 

. 001198 

80.  5 

.  17073 

.  126 

-26.  2 

-.  12 

.01416 

80.  7 

2.  018 

2.  033 

+.74 

+.  04 

.  03530 

86.  8 

5.  033 

5.  150 

+2.  32 

+.  31 

.07890 

85.  5 

11. 249 

11.  40 

+  1.  34 

+.  40 

.  1072 

85.  6 

15.  284 

15.  40 

+.  76 

+.  31 

1  .1406 

85.  6 

20.  039 

20.  169 

+.  65 

+.  34 

.  00000 

! 

0.  000 

-0.  06 

-.  002 
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Test  Mode  1:  FUEL  FLOW  CALIBRATION  TEST  PLAN 


a 


1.  Turbine  flow  meter 

2.  Flow  straightener  (optional) 

3.  Sight  port  (optional) 

4.  SP  densitometer 

5.  Panametrics  flow  velocimeter 


6.  Sight  port  (optional) 

7.  Kenics  static  mixer  (optional) 

8.  Anti-vibration  mount 

9.  Temperature  sensor 

10.  Pressure  Tap 


JP-5  DATA  TO  BE  RECORDED 


GE  Mass  Flow  Rate.PPH  Panametrics  Mass  Flow  Rate,  PPH 
Pressure,  psig  Turbine  Corrected  for  T  p  V  Calc.  M  Meas.lvl 


Test  Mode  2:  WATER  FLOW  CALIBRATION  TEST  PLAN 

Test  arrangement  as  above  except  include  only  items  5,  9,  10;  use 
water  as  fuel  substitute. 

Record  data  on  actual  flow  velocity  (determined  from  weigh  tank  or 
volumetric  measurement)  and  on  ultrasonically-determined  flow  velocity. 

Test  Mode  3:  10- DAY  REPEATABILITY  TEST  PLAN 

Test  arrangement  like  in  Mode  2  except  use  head  difference  AH  in¬ 
stead  of  fuel  pump,  and  use  axial  path  offset  cell  instead  of  zigzag  cell  as 
item  5.  Record  data  as  in  Mode  2. 
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B.  RANGEABILITY 

The  maximum  (water)  flow  rate  at  which  the  flow  velocimeter  was 
tested  in  this  program  was  38  ft/s  (11,6  m/s).  (Higher  flows  could  prob¬ 
ably  be  measured  since  the  received  signals  at  38  ft/s  show  no  large  jitter. 
See,  for  example,  oscillograms  in  Fig.  22  for  JP-5  at  20  kPPH  and  water 
at  60  kPPH.  )  The  minimum  flow  rate  was  zero.  However,  the  zero  offset 
observed  prior  to  testing,  during  testing,  and  again  after  the  equipment 
was  returned  from  Foxboro  to  Panametrics,  was  0.06  ft/s.  This  zero  off¬ 
set  error,  V  ,  corresponds  to  a  At  offset  error  of 
e  e 

At  =  2LV  /c2 
e  e 

=  (2)  (0.  33  ft)(0.  06  ft/s )/ (5  x  103  ft/s)2 
=  1 . 6  ns 


which  is  close  to  the  +1  ns  resolution  limit  of  the  flow  velocity  electronics. 


this  V 


The  fuel  mass  flow  rate  zero  offset  error, 
and  At  ,  is  60  PPH. 


M  ,  corresponding  to 


The  M  objective  in  this  contract  was  from  200  to  25000  PPH,  a 
rangeability  of  125:1.  However,  due  to  the  density  range,  ~  0.  65  to  0.  9  g/cm 
the  corresponding  velocity  range  objective  is  thus  (0.  9/0.  65 )  (1 25)  =  173:  J. 

Now  in  the  Foxboro  calibration  test  (Table  6;  Fig.  30),  the  nonzero  flow 
range  was  0.  1  to  38  ft/s,  demonstrating  a  velocity  rangeability  of  380:  1. 

Due  to  the  zero  offset,  accuracy  over  this  entire  range  is  not  adequate  for 
all  purposes;  hence,  the  useful  rangeability  depends  on  the  required  accu¬ 
racy.  It  appears  that,  within  the  limits  of  one  zigzag  cell,  useful  range- 
ability  in  excess  of  100:1  was  achieved.  For  better  accuracy  at  the  lower 
flow  rates,  a  longer  L  is  required,  e.  g.  ,  as  provided  by  an  offset  axial 
path  flow  cell  (Appendix  G). 


C.  REPEATABILITY 


Reference  to  Table  6  shows  that  repeatability  at  zero  flow  was 
0.  06  ft/s,  or  0.  2%  full  scale  based  on  the  38  ft/s  maximum  velocity,  or 
0.  24%  based  on  a  full  scale  flow  velocity  of  25  ft/s  (fuel  lli  =  25  kPPH). 

D.  HYSTERESIS 


Reference  to  Table  6  shows  that  hysteresis  was  less  than  1.2%  of 
reading  (5  -  20  ft/s)  and  less  than  0.  25%  full  scale.  (The  values  listed  be¬ 
low  are  extracted  from  Table  6,  for  the  convenience  of  the  reader.  ) 
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Nominal 

%  of  Reading 

Error 

Hystereais 

Flow  Velocity 

Flow 

Flow 

(Difference) 

y- . » 

Deer. 

Incr. 

%  of  Reading 

%  Full  Scale 

+  1.14 

+  2.  32 

1.  18 

.  24 

+  .  27 

+  1.  34 

1.  07 

.  11 

*  .  52 

+  .  76 

.  24 

.  02 

t  . 

+  .  65 

.  29 

.  02 

r..  l.INLAHITY 

The  linearity  is  1 %  full  scale.  Note  that  the  "meter  factor"  was  set  to 
0.  in  the  water  flow  calibration  tests  and  agreement  between  the  ultra¬ 

sonic  measurement  of  flow  velocity,  and  that  determined  independently  by  the 
weigh  tank  method,  was  within  0.  5%.  This  demonstrates  the  area-averaging 
property  of  the  square  duct  when  100%  interrogated  as  described  elsewhere 
in  this  report.  No  flow  straightener s  or  static  mixers  were  used  in  this 
water  test.  In  the  water  calibration  tests,  and  in  the  final  fuel  substitute  tests, 

the  deviation  from  linearity,  expressed  as  percent  of  actual  flow,  is  larger  at 
the  bottom  portion  of  the  system's  range.  See  for  example,  p.40  or  Appendix  K. 

F.  RESPONSE  TIME 

In  the  delivered  instrument,  the  M  response  time  is  fixed  at  2  s. 
However,  for  V  only,  a  board  is  available  for  the  instrument  to  provide 
switch- selectable  averages  of  1,  2,  4,  8  ....  128  readings.  In  the  vi¬ 

bration  tests  at  ATL,  the  response  time  of  2  seconds  was  demonstrated 
by  rapidly  changing  a  valve  setting*  See  Appendix  E.  Referring  to  the  chart 
recording  on  p.  E-l,  note  that  deflections  for  each  cell  are  measured  from  a 
"V  =  0"  no-flow  reference  axis.  The  deflection  sensitivity  is  0.3  inches  (30  mV) 
per  English  unit  of  flow  velocity.  That  is,  30  mV  corresponds  to  V  =  1  ft/s 
(0.3  m/s).  When  the  water  flow  control  valve  is  suddenly  opened  or  closed, 
the  response  time  of  the  velocimeter  electronics  can  be  determined  by  mea¬ 
suring  how  long  the  output  takes  to  steady  down  (equilibrate)  at  the  new  flow. 

At  a  chart  speed  of  1  inch  per  minute,  1  minor  division  =  6  seconds.  Note 
thatthequasi-step  changes  in  the  output  are  completed  in  T  <  2s,  as  when 
flow  increased  from  1  to  2.5  ft/s,  or  when  the  valve  was  rapidly  opened  by 
three  turns.  In  other  laboratory  tests  at  Panametrics,  quasi-step  changes 
in  water  flow  rate,  produced  by  closing  a  ball  valve  in  less  than  1  s,  lead  to 
the  same  value  of  T,  even  when  M  changes  over  the  50:1  range  from  25,  000 
down  to  500  lb/hr.  If  the  flow  is  turned  off  suddenly,  the  output  reading  nearly 
reads  zero  immediately,  but  does  not  fully  settle  until  the  eddies  subside. 

See  also,  pages  G-10  and  G-ll. 
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G.  Vibration 


ATL's  vibration  test  report  comprises  most  of  Appendix  J. 
This  is  followed  by  a  drawing  to  define  x,  y  and  z  axes  of  vibration 
and  recordings  of  the  flow  velocity  data  during  vibration.  It  will  be 
understood  that  for  another  flow  cell  configured  to  and  mounted  on  a 
particular  engine,  results  will  be  different. 


Vibration  testing  at  ATL,  Burlington,  MA,  utilized  two  cells  in 
series  as  shown  in  the  Test  Plan.  Cold  water,  at  the  maximum  flow  rate 
available  from  the  tap,  up  to  3  ft/s  (~1  m/s),  flowed  first  through  the 
cell  being  vibrated  and  then  through  a  second  cell  at  rest  (Mode  1). 

Flow  velocimeter  electronics  monitored  the  velocities  V^  and  V^  respec¬ 
tively  in  the  two  cells.  The  influence  of  vibration  on  the  measured  V^ 
was  found  to  be  small,  in  general.  Besides  resonances  and  noise,  the 
only  effect  of  vibration  to  be  expected  from  theoretical  considerations 
occurs  when  the  vibration  frequency  f  in  the  axial  direction  is  approxi¬ 
mately  equal  to  the  frequency  at  which  upstream  and  downstream  inter¬ 
rogation  directions  alternate  electronically.  Such  vibratory  synchronism 
produces  an  additional  component  of  relative  velocity  V  between  the  cell 
and  the  fluid,  directly  proportional  to  g  and  inversely  proportional  to  f  . 
That  is,  the  peak  relative  velocity  V  =  7T  f  S  =  61.44  g/f  where  the 
units  of  V  are  in.  /  s,  f  ,  Hz;  S,  /tin.  (double  amplitude  displacement); 
and  g  is  tl?e  gravitational  constant  (g  =  0.  0511  f  S  =  0.  0162  V  f  ).  For 
example,  at  f  =  2000  Hz,  and  20  g’s,  S  =  100  /Jin.  (2.  54  film)  §n$ 

V  =  0.  64  in.  / s  (~16  mm/s).  Further  details  on  the  vibration  test  per¬ 
fumed  on  the  delivered  flow  velocimeter  cell  are  given  in  the  Vibration 
Test  Plan  and  in  Appendix  J.  No  damage  was  observed  due  to  the  vibra¬ 
tion  test.  However,  at  a  few  discrete  frequencies,  e.  g.  f  ~  1250  Hz, 
noise  was  so  intense  that  flow  data  had  to  be  discarded.  Some  future  work 
could  be  devoted  to  eliminating  such  undesirable  resonant  behavior. 


H.  Temperature 


A  transducer  identical  in  construction  to  the  two  in  the  delivered 
flow  cell  was  tested  in  an  environmental  chamber  in  the  pulse-echo  mode 
from  -196  to  +175  C.  Testing  included  cycling  and  dwells  of  several 
hours  at  the  extremes.  Both  low  and  high  extremes  are  beyond  the  con- 
tract-specified  values  of  -38  to  +150  C.  The  oscillograms  in  Appendix  I 
show  that  the  tested  transducer  exceeded  the  temperature  objectives  of 
the  contract. 


I.  Pressure 


The  flow  velocimeter  cell  was  hydrostatically  tested  at  1500  psig 
for  one  hour  and  was  found  to  be  leak  tight  (see  Appendix  H). 
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VIBRATION  TEST  PLAN 


1.  Water  pressure  source  and  drain 

2.  SP  p  electronics 

3.  Panametrics  electronics  (V  and  M) 

4.  Optional  Panametrics  flow  velocimeter  electronics 

5.  SP  p  cell  filled  with  fuel  or  fuel  substitute 

6.  SP  p  cell  through  which  water  flows* ** 

7.  Panametrics  V  cell 

8.  Optional  V  cell  not  under  vibration 

*Test  mode  1:  compare  V  and  V1  as  function  of  vibration,  at  constant  water 
flow  conditions. 

*Test  mode  2:  compare  p  and  V  and/or  M  readings,  under  vibration  and  flow, 
and  then  when  vibration  and/or  flow  is  momentarily  turned  off. 

DATA,  MODE  1:  V,  V'  &  H^O  flow  conditions 

DATA,  MODE  2:  p,  V  and/or  M  for  vibration  and/or  flow  =  zero  and  4  zero. 


*Test  mode  1  was  selected  for  the  ATL  tests  conducted  Sept.  1980. 

**SP  cells  may  be  eliminated  from  tests  to  avoid  water  damage;  SP  density 
signal  may  be  simulated  using  dc  power  supply  to  represent  p  =  1  g/cml 
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TEMPERATURE  TEST  PLAN 


Test  Mode  Is 


Test  Mode  2: 


Test  Mode  1 


1.  Flow  cell  full  of  liquid  (water  or  hydrocarbon) 

2.  Pulser /receiver,  Panametrics  model  5055,  or  flow  velocimeter  electronics. 

3.  Oscilloscope 

Procedure:  measure  received  signal  before  and  after  exposure  to  contra  ctually- 
specified  temperature  extremes. 

DATA:  Received  signal  (mV),  before _ ;  after  -65°F _ ;  after  304  F _ 


Test  Mode  2 


1.  Thick  window  transducer  tested  individually,  not  in  flow  cell. 

2,3.  Electronics  same  as  for  the  items  for  Test  Mode  1. 

Procedure:  measure  relative  amplitude  of  echo  from  end  of  thick  window  at  room 
temperature  and  then  during  exposure  to  contractually- specified  temperature  ex¬ 
tremes,  or  at  temperatures  beyond  these  extremes. 

HYDROSTATIC  PRESSURE  TEST  PLAN 


1.  Plug  flow  cell  outlet. 

2.  Pressurize  from  inlet,  using  water,  at  1600  psig.  Hold  for  one 
hour.  Report  any  leaks.  If  necessary,  re-seal  and  test  again. 
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J.  Self-Testing  of  Flow  Velocimeter  Electronics:  Differential  Path 

There  are  a  number  of  approaches  to  self-testing  parts  of  the  flow 
velocimeter  electronics.  Among  the  more  obvious  approaches  may 
be  mentioned  the  electronic  simulation  of  a  particular  flow  condition  as 
represented  by  a  preset  At  or  Af  injected  into  the  system,  and/or  the  use 
of  numerous  diagnostic  trouble  lights  (led1  s)  to  indicate  malfunction  of 
particular  circuits.  The  problem  with  the  former  is  that  the  transmitter 
and  perhaps  parts  of  the  receiver  are  not  included  in  the  test.  An  objection 
to  the  latter  is  that  if  the  operator  in  the  field  is  unable  to  repair  the  equip¬ 
ment  based  on  the  led  indications,  then  the  detailed  troubleshooting  infor¬ 
mation  provided  by  them  may  be  superfluous. 

In  order  to  self-test  the  complete  flow  velocimeter  electronics 
independent  of  the  actual  flow  cell,  and,  as  desired,  independent  of  or 
including  the  actual  transducers,  consider  the  following  differential  path  or 
interferometric  type  self-test  cells.  These  cells  may  use  two  or  four  trans¬ 
ducers  which  are  generally  fixed  in  a  standard  pipe  fitting  housing.  The 
housing  itself,  or  a  reflector  therein,  is  arranged  to  provide  different  fixed 
or  variable  paths  simulating  upstream  vs  downstream  interrogations.  Pipe 
section  extensions  may  be  connected  to  the  housing  to  increase  the  fluid 
paths  to  closely  approximate  the  delays  encountered  in  any  specific  application. 
The  self-test  cell  may  be  filled  with  water  or  other  fluid  of  interest,  provided, 
of  course,  that  the  necessary  precautions  are  taken  if  fuel  or  other  hazardous 
fluid  is  used. 

We  briefly  explain  four  self-test  cells  as  follows.  Referring  first  to 
Fig.  23a  note  the  use  of  two  fixed  transducers  and  a  reflector  between  them. 

The  transducers  are  operated  each  in  pulse  echo  mode,  the  reflector  blocking 
transmission  between  the  transducers.  Displacement  of  the  reflector  by  a 
small  distance  Ax  increases  the  "upstream"  simulated  delay  by  2  A x/c  and 
decreases  the  "downstream"  simulated  delay  by  the  same  amount.  Therefore, 
the  simulated  upstream  minus  downstream  time  difference  changes  by  4Ax/c. 

Depending  on  the  nature  of  the  pulser  and  receiver,  it  is  not  always 
convenient  to  operate  the  electronics  in  pulse  echo  mode.  For  a  through 
transmission  mode  self-test,  one  can  consider  the  self-test  cells  of  Fig.  23b 
and  c.  In  Fig.  23b,  the  differential  path  is  again  denoted  Ax,  the  path  dif¬ 
ference  simply  being  the  difference  in  fixed  distances  between  transmitters 
and  receivers.  The  time  difference  is  simply  Ax/c. 

Reference  to  F  ig.  23c  shows  that  displacement  of  the  reflector  along 
the  x-axis  by  Ax  changes  each  45°  folded  path  by  Zsfl  Ax.  If  the  test  fluid  is 
water,  the  flow- simulated  At  is  about  2gs/mm  or  50ns/mil  (I  mil-  10’3  inches) 
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It  may  be  of  interest  to  note  that  if  the  pulse  echo  idea  of  Fig.  23a 
is  combined  with  the  "geometry"  of  Fig.  23c  in  the  sense  of  having  the 
reflector  face  at  45°  to  the  displacement  axis,  then  the  same  At' s  are 
generated  with  this  cell  in  pulse  echo  or  through  transmission,  namely, 
ZsfZ  Ax/c.  See  Fig.  23d. 

With  respect  to  At,  the  simulation  of  flow  may  be  understood  by 
equating  actual  and  simulated  time  differences: 


At 


2LV 


Z*JZ  Ax 
c 

s 


where  c  =  sound  speed  in  actual  fluid  (e.  g.  JP-5)  and  c  =  sound  speed 

cL  S 

in  simulation  fluid  (e.  g.  water). 

Solving  for  the  displacement  corresponding  to  a  given  V  and  L, 


Ax 


As  a  numerical  example  corresponding  to  the  delivered  flow  cell  and 
the  maximum  flow  velocity  of  interest  (at  25000  lb /hr),  L  =  4  inches  =  0.333  ft 
(10  cm),  V  =  25  ft / s  (~  8  m's),  cg  =  5000  ft/ 8  (~*  1500  m's)  and  c  =  2500  ft's 
(•-'750  m/s).  Then 


A  (50  J0)(Q,  3  3)( 25 ) 
AX  (2500)(2500)  sfZ 


U.  057  inches  =  1.4  mm 


A  rapid  change  in  reflector  position  could  be  used  to  check  the  transient 
response  of  the  electronics.  For  example,  a  rapid  change  of  Ax  =  1.4  mm  in 
0.1s  could  simulate  aV  =  8m/s  in  0.1  s,  equivalent  to  a  fuel  acceleration  of 
AV/At  =  8m/s^  or  about  0.8g.  Similarly  a  rapid  path  change  of  AX  =  17  mm 
(~0.  7”) in  0.1  8  would  simulate  a  fuel  acceleration  of  lOg. 
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V.  CONCLUSIONS 

The  complete  delivered  M  system  was  tested  at  GE  in  the  fuel 
substitute  7024  BII  to  14,  300  kg/hr  (31,450  lb/hr).  The  V  part  of  the 
system  was  also  tested  in  water  to  V  =  11.6  m/s  (38  ft/s)  ot  M  s  27000 
kg/hr  (60,000  Ib/hr).  Based  on  a  full  scale  rating  of  11,364  kg/hr 
(25,  000  lb/hr)  rms  deviations  from  linearity  in  both  tests  were  demon¬ 
strated  to  be  about  1%  of  full  scale,  with  a  2  s  response  time.  Over 
approximately  the  upper  95%  of  the  M  range,  in  the  fuel  substitute 
tests  at  GE,  the  rms  deviation  from  linearity  is  at  most  1.7%  of  the 
actual  flow  rate.  Up  to  0.5%  of  this  nonlinearity  may  be  due  to  the  turbine 
sensors'  own  error  band.  Above  9000  Ib/hr,  for  3  runs,  the  average 
nonlinearity  is  less  than  0.  77%  of  the  actual  flow  rate.  The  performance 
of  the  mass  flow  rate  system  is  believed  to  be  substantially  independent  of 
fuel  composition  (at  least  for  JP-4,  -5,  their  mixtures,  or  7024  BII), 
and  substantially  independent  of  temperature  from  -55  to  +300°F,  pressure 
to  1600  psig  and  vibration  (MIL  Std.  810-514.1),  except  at  a  few  discrete 
frequencies.  Ten-day  repeatability  testing  between  ~0.5  and  ~5  ft/s 
(~0.  15  to  1.5m/s)  demonstrated  system  repeatability  of  ;t  0.01  ft/s  (4^3mm/s). 
The  flow  velocimeter  part  of  the  electronics  can  be  operated  with  a  variety 
of  flow  cells  such  as  clamp-on  or  axial  offset  or  conventional  diagonal 
traverse,  in  addition  to  the  square  hole  zigzag  design  with  which  the  present 
contract  was  mainly  concerned.  The  stability  of  the  waveforms  observed 
for  the  zigzag  flow  cell  at  the  maximum  M  encountered  in  the  fuel  substitute 
test,  and  in  the  water  tests,  strongly  suggests  that  a  longer  zigzc.g  path 
could  be  used  in  the  future.  This  would  improve  the  resolution  at  low  M, 
and,  together  with  some  other  modifications  to  the  cell  and  transducer, 
ought  to  improve  linearity  in  the  lower  partof  the  system's  range.  Alterna¬ 
tively,  linearization  may  be  accomplished  in  the  software,  at  least  for 
nonlinearities  dependent  on  V  alone. 
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Fig.  1, 

Fig.  2 

Fig.  3 

Fig.  4. 

Fig.  5 

Fig.  6 

Fig.  7 
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VIIL  LIST  OF  ILLUSTRATIONS 


Ultrasonic  mass  flowmeter  cell  shown  in  exploded  view.  Cell 
consists  of  aluminum  housing  milled  to  provide  1/2"  (2.  54  mm) 
square  duct  and  interrogation  path  at  45°,  5  MHz  transducers 
for  sensing  flow  velocity,  another  ultrasonic  transducer 

for  sensing  the  speed  of  sound  by  a  pulse  echo  technique,  and  a 
reflection  coefficient  densitometer  probe  using  a  10  MHz  trans¬ 
ducer.  Sponsor:  Eustis  Directorate,  1971-1973. 

Assembly  view  of  ultrasonic  mass  flow  cell  and  block  diagram 
of  flow  velocimeter  electronics  utilizing  a  pseudo  random  noise 
code  generator  which  permitted  the  fuel  to  be  interrogated  simul¬ 
taneously  upstream  and  downstream  at  the  same  carrier  frequency, 
5  MHz,  using  only  one  pair  of  transducers. 

Block  diagram  of  5  MHz  rf  burst  system  used  for  measuring  flow 
velocity  in  third  Eustis  program. 

(a)  Flow  cell  configured  for  GE's  T-700  engine  including  passage¬ 
ways  for  the  ultrasonic  measurement  of  flow  velocity  and  the 
capacitance-based  measurement  of  fuel  density,  (b)  Field-type 
case  which  housed  the  electronics  for  mass  flow  rate  determination. 
Source:  USARTL-TR -78-45  (October  1978). 

Reynolds  number  nomogram  for  determining  Re  »s  a  function  of 
kinematic  viscosity  and  the  product  of  area  averaged  flow  velocity 
and  duct  diameter  D.  Curves  at  left  show  the  temperature 
dependence  of  kinematic  viscosity  for  a  number  of  liquids  and 
gases.  Adapted  from  Physical  Acoustics  -  Vol.  14,  p.  490  (1979) 
including  adaptations  from  Streeter  (1961)  and  Moody  (1944)  cited 
therein.  Dashed  box  indicates  scope  of  present  contract. 

Graphs  showing  (a)  sound  speed  vs  temperature,  (b)  density  vs 
temperature,  (c)  kinematic  viscosity  vs  temperature,  and  (d) 
kinematic  viscosity  vs  sound  speed,  for  JP-4  and  JP-5 

—  •  3 

V  vs  M  for  p  -  0.  65  and  0.  9  g/cm  . 

Block  diagram  of  phase  locked  loop  flow  velocimeter,  operating 
in  conjunction  with  SP  densitometer  to  measure  mass  flow  rate. 
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VIIL  LIST  OF  ILLUSTRATIONS  (cont'd) 


Fig.  9. 
Fig.  10. 
Fig.  11. 

Fig.  12. 
Fig.  13. 

Fig.  14. 

Fig.  15. 

Fig.  16. 

Fig.  17. 

Fig.  18. 
Fig.  19. 

Fig.  20. 


Comparison  of  various  types  of  ultrasonic  flow  cells. 

Small,  medium  and  large  flow  cells. 

(a)  Axial  path  interrogation  of  liquid  contained  within  SP  concentric 
tube  densitometer  cell,  (b)  Oscillographic  test  results  for  SP  cell 
and  (c)  Comparison  with  unobstructed  axial  path.  Liquid:  water, 
room  temperature.  Electronics:  Panametrics  Puiser /Receiver 
5055  PRM  +  Krohn-Hite  Filter  3202,  1.5  MHz  high  pass. 

Flow  cells  made  of  standard  tubing  pipe  and  compression  fittings. 

Navy  flow  velocimeter  cell  photograph  showing  optional  ports  for 
set  screws  holding  square  holed  sleeve  and  for  optional  slow  torsional 
wave  densitometer. 

Outline  drawing  of  equipment  shown  in  Fig.  13  in  which  the  key 
components  are  identified. 

Outline  drawing  of  basic  flow  cell  in  which  optional  ports  are  omitted. 
The  position  of  the  square  holed  sleeve  and  the  two  ultrasonic  trans¬ 
ducers  are  indicated  by  dashed  lines. 

Alternate  design  for  flow  velocimeter  cell  portion  of  the  system, 
using  transducer  ports  orthogonal  to  flow  axis  and  reflectors  to 
control  45°  zigzag  interrogation  path. 

Photograph  showing  examples  of  pc  reflectometer  probes  similar 
in  principle  to  flush-mounted  pc  densitometer  probe  investigated  in 
first  and  second  Eustis  contracts. 

M  system  block  diagram. 

Photograph  of  prototype  flow  velocimeter  electronics  similar  to  the 
instruments  which  were  used  to  test  the  reference  (at  rest)  flow  veloci¬ 
meter  cell  in  the  vibration  tests  and  the  final  offset  flow  cell  tests  at  GE. 

Graph  of  insertion  loss  vs  thickness  of  sheets  that  were  considered 
as  eddy  reducers  in  vicinity  of  transducer  ports.  Curves  are  based 
on  calculations  for  normal  incidence  longitudinal  waves  at  5  MHz 
and  assume  that  the  sheet  is  fully  immersed  and  wetted  by  water  at 
room  temperature. 
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VIII.  LIST  OF  ILLUSTRATIONS  (cont'd) 

Fig.  21.  Accuracy  nomogram. 

Fig.  22.  Test  oscillograms  of  received  signals  transmitted  through  JP-5 
at  various  flow  rates. 

Fig.  23.  Flow  simulators  comprised  of  differential  path  cells,  (a)  Pulse 
echo,  (b)  Through  transmisstion.  (c)  Folded  transmission,  or 
"pitch-and- catch, "  arrangement,  (d)  Pulse  echo  with  reflector 
motion  at  45°  to  interrogation  axis. 

Fig.  24.  General  view  of  flowmeter  under  preliminary  test  at  GE-Lynn. 

Photo  courtesy  GE. 

Fig.  25.  Close-up  of  mass  flowmeter  cell  in  preliminary  test.  In  subsequent 
tests,  distance  from  elbow  to  cell  was  increased,  allowing  room 
for  either  flow  straightener  or  static  mixer  to  be  installed  in-line 
upstream.  Photo  courtesy  GE. 

Fig.  26.  Schematic  of  test  at  GE-Lynn,  March  1980 

Fig.  27.  Mass  flow  rate  calibration  data  using  fuel  substitute.  Reference 
M  obtained  using  Cox  ANC-24,  S/N  23803,  calibrated  from  990  to 
59000  pounds  per  hour,  accuracy  +  1/4%  of  reading. 

Fig.  28.  AH  flow  loop  with  offset  flow  cell  used  in  repeatability  test. 

Fig.  29.  Data  for  10-day  repeatability  test,  using  water  in  offset  flow- 
cell,  July  1980. 

Fig.  30.  Water  calibration  test  data  obtained  by  the  Foxboro  Company, 
Foxboro,  MA,  on  23  and  24  September  1980. 
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Fig.  1.  Ultrasonic  mass  flowmeter  cell  shown  in  exploded  view  Cell  con¬ 
sists  of  aluminum  housing  milled  to  provide  1/2"  (2,54mm)  square  duct  and 
interrogation  path  at  45°.  5  MHz  transducers  for  sensing  flow  velocity,  another 

ultrasonic  transducer  for  sensing  the  speed  of  sound  by  a  pulse  echo  technique, 
and  a  reflection  coefficient  densitometer  probe  using  a  10  MHz  transducer. 
Sponsor:  Eustis  Directorate,  1971-1973. 
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Fig.  2.  Assembly  view  of  ultrasonic  mass  flow  cell  and  block  diagram  of 
flow  velocimeter  electronics  utilizing  a  pseudo  random  noise  code  generator 
which  permitted  the  fuel  to  be  interrogated  simultaneously  upstream  and  down 
stream  at  the  same  carrier  frequency,  5  MHz,  using  only  one  pair  of  trans¬ 
ducers. 
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Fig.  5.  Reynolds  number  nomogram  for  determining  Re  as  a  function  of  kinematic  viscosity  and  the  product  of 
area  averaged  flow  velocity  V  and  duct  diameter  D.  Curves  at  left  show  the  temperature  dependence  of  kinematic 
viscosity  for  a  number  of  liquids  and  gases.  Adapted  from  Physical  Acoustics  -  Vol.  14,  p.  490  (1979)  including 
adaptations  from  Streeter  (1961)  and  Moody  (1944)  cited  therein.  Dashed  box  indicates  scope  of  present  contract. 
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Fig.  6.  Graphs  showing  (a)  sound  speed  vs  temperature,  (b)  density  vs  tern 
perature,  (c)  kinematic  viscosity  vs  temperature  and  (d)  kinematic  viscosity 
vs  sound  speed,  for  JP-4  and  JP-5.  50 
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Fig.  9.  Comparison  of  various  types  of  ultrasonic  flow  cells.  Features 
or  characteristics  of  these  different  types  of  cells,  plus  a  few  others,  and 
identification  of  types  IB,  2A  and  8A  as  the  most  commonly  used  types  at 
this  time,  are  discussed  in  ref.  10. 
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SCHEMATIC  OF  PANAMETRICS  FUEL  MASS  FLOWMETER 
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Fig.  18.  M  system  block  diagram. 
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Fig.  19.  Photograph  of  prototype  flow  velocimeter  electronics  similar  to  the  instruments  which  were  used 
to  test  the  reference  (at  rest)  flow  velocimeter  cell  in  the  vibration  tests  at  ATL  and  the  axial  path  offset  flow 
cell  in  the  final  fuel  substitute  tests  at  GE. 
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AXIAL  PROJECTION  OF  PATH  IN  FLOW,  L  ,  METERS 


FLOW  VELOCITY  V,  m/s 


It  is  virtually  impossible  to  cha racteri/.e  the  performance  of  an  ultrasonic  flowmeter 
with  just  one  number.  Even  if  one  agrees  on  terminology  for  accuracy,  repeatability 
response  time,  etc.  ,  the  performance  still  may  depend  on  factors  beyond  the  control 
of  the  manufacturer.  Conditions  upstream  or  downstream,  type  of  fluid(s),  and  the 
temperature  and  pressure  range,  ami  other  application  factors  such  as  pipe  size,  wall 
material,  thickness  and  roughness  differ  from  case  to  case  to  case. 

In  most  cases  the  ultimate  limit  on  performance  is  imposed  by  the  repeatability  error. 
In  the  standard  Model  6000  electronics,  this  limit  is  typically  1  ns.  But  to  be  mean¬ 
ingful  this  must  be  translated  into  a  flow  error.  The  above  graph  was  constructed  so 
the  user  can  quickly  and  easily  determine  the  flow  repeatability  error  to  be  expected 
from  a  Model  6000  flowmeter,  for  a  given  set  of  conditions. 

Enter  the  left-hand  graph  at  c  and  E.  On  left  ordinate  read  error  in  V,  per  nanosecond 
error  in  At. 

Right-hand  graph  yields  the  percentage  error,  such  as  %  of  reading,  %  or  span  or  % 
full  scale,  according  to  the  V  (abscissa)  value  selected.  Example  1:  water,  c=  1500  m/s, 
L=  0.15m(6"),  V  »  3m/s  (=:10ft/s).  Assuming  error  in  At  is  Ins,  read  error  in  V 
as  7. 5  mm  / s  *  0. 25%  of  reading.  If  full  scale  were  lOm/s  read  error  in  V  as  0.075%  full 
scale.  Example  2:  JP  - 5 ,  c  =  1000  m/s,  L  =  10cm  (4" ),  V=7.62m/s(25ft/s).  On  left- 
hand  graph  read  error  in  V  as  5  mm/s  per  ns.  On  right-hand  graph  estimate  error  as 
0.07%  of  reading,  per  ns  error  in  At.  If  error  in  At  =15  ns,  then  error  in  V  would  be 
about  1%  of  reading.  Note  use  of  graduations  to  assist  logarithmic  interpolation. 
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Fig  22.  Test  oscillograms  of  received  signals  transmitted  through  JP-5 
at  GE  at  various  flow  rates  are  shown  in  top  two  recordings.  Bottom  two 
recordings  were  obtained  at  the  beginning  of  the  water  calibration  tests  at 
the  Foxboro  Company,  23  September  1980,  at  a  velocity  of  38  ft/s  (M  = 

60,  000  lb/hr).  Observe  that  At  =  1.1/isforL  =  100  mm  and  V  =  38ft/s. 
Jitter  appears  negligible  in  the  "double  exposure"  of  the  upstream  and  down¬ 
stream  received  signals.  Compare  with  Fig.  K4,  top  oscillogram,  bottom 
trace,  obtained  using  the  fuel  substitute  7024BII  at  31,450  lb/hr(14,309  g/  r) 


Fig.  23.  Flow  simulators  comprised  of  differential  path  cells,  (a)  Pulse  echo. 

(b)  Through  transmission,  (c)  Folded  transmission,  or  "pitch-and-catch, "  arrange 
ment.  (d)  Pulse  echo  with  reflector  motion  at  45°  to  interrogation  axis. 
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Fi>j.  28.  AH  flow  loop  with  offset  flow  cell  used  in  10-day  repeatability 
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Fig.  29.  Data  for  10-day  repeatability  test,  using  water  in  offset  flow  cell, 
July  1980.  ID  =  7.9  mm.  L  =  1 32  mm.  P  =  168  mm,  t  =  1  (jls.  At 
the  maximum  V,  a  =  0.2%  of  reading;  at  minimum  V,  <T  -  2%  reading 


11.  6  m /s 


ACTUAL  VELOCITY,  FT/S 

Fig.  30.  Water  calibration  test  data  obtained  by  the  Foxboro  Company,  Foxboro, 
MA,  on  23  and  24  September  1980- 

Flow  cell  type:  same  as  previous  figure.  Above  data  obtained  after  eub- 
jecting  cell  and  transducers  to  Mil  Spec  vibration  test  5  to  2000  H/.  and 
up  to  20  g /s. 
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APPENDIX  A.  OBJECTIVES /EXTRACTS  FROM  NADC  STATEMENT 
OF  WORK 


3.  1  Fuels:  The  transmitter  must  be  capable  of  the  mass  flow 
rates  of  JP-4  and  J P-5  aviation  fuels  and  all  mixtures 
thereof  without  manual  calibration. 

3.  1.  1  Fuel  Temperature  Range:  The  flowing  fuel  temperature 
range  may  vary  from  -38°C  to  +150°C. 

3.  1.  2  Fuel  Substitution:  An  approved  commercial  fluid  of  higher 
flash  point  but  same  specific  gravity  may  be  used  for  test 
and  calibration. 

3.  1.  3  Contaminated  Fuel:  The  transmitter  shall  be  operated  as 
specified  in  sec.  4.  10.  No  degradation  in  subsequent 
accuracy  tests  shall  be  evidenced  or  damage  observed 
upon  inspection. 

3.  2  Mass  Flow  Rate  Range:  The  transmitter  shall  be  capable  of 
accepting  a  flow  rate  range  and  transmitting  signals  ana¬ 
logous  to  flow  rates  of  100  PPH  to  25,  000  PPH. 

3.  3  Accuracy:  Over  the  flow  range  stated  in  item  3.  2,  the 

transmitter  accuracy  shall  be  +  5%  of  point  at  low  end  of 
flow  range  and  linearly  improving  to  +_  0.  5%  of  point  at 
3,  000  PPH.  Then  it  shall  remain  at  this  value  over  the 
rest  of  the  flow  range  to  maximum  flow  (25,  000  PPH). 

3.  4  Hysteresis:  Upscale  and  downscale  readings  shall  agree 
to  within  +_  0.  3%  of  each  other. 

3.  5  Repeatability:  Repeatability  is  here  defined  as  the  capa¬ 
bility  of  reproducing,  at  successive  points  in  time,  the 
same  signal  output  for  the  same  mass  flow  input. 

The  variations  shall  not  exceed: 

+_  2  PPH  at  low  end  of  flow  range 
+  3  PPH  at  test  point  of  3,  000  PPH 
+_  15  PPH  at  test  point  of  15,  000  PPH 
+_  25  PPH  at  test  point  of  25,  000  PPH 
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3.  6  Response  Time:  The  transmitter  shall  be  capable  of 

being  rapidly  slewed  from  500  to  25,  000  PPH  and  having 
its  output  come  to  within  +  1%  of  final  steady- state  value 
in  two  seconds. 

3.7  Signal  Output  Format:  The  transmitter  signal  output  shall 
be:  a  0  to  5  volt  square  wave  with  a  frequency  directly 
proportional  to  mass  rate  of  flow  (M)  and  further  shall  be 
compatible  with  the  requirements  of  MIL-STD-  1553  for 
multiplex  operation. 

3.  8  Vibration:  The  transmitter  and  any  associated  signal 
processing  circuitry  shall  be  capable  of  meeting  the 
accuracy  requirements  of  3.  3  while  under  vibration  at 
room  temperature  as  in  sec.  4.  9. 

3.  9  Power:  The  transmitter  shall  be  designed  and  built  to 
accept  28  volts  D.  C.  as  defined  by  MIL-STD-704. 

3.  9.  1  Optional  Power  Requirement:  In  addition,  the  design, 

but  not  the  present  fabricated  product,  shall  incorporate 
the  option  of  being  able  to  operate  on  400  Hz,  115  volt 
power  as  defined  by  MIL-STD-704. 

3.  10  Desired  Design  Attributes 

3.  10.  1  Wiring:  To  consist  of  not  more  than  4  conductors,  2  for 
power,  ?  for  signal.  Signal  and  power  common  shall  not 
be  carried  on  the  same  conductor.  Design  should  in¬ 
corporate  the  option  for  optical  signal  transmission. 

3.  10.  2  Flow  Connectors:  All  fluid  transmitting  connectors 

shall  be  of  the  quick  disconnect  type,  sized  to  accom¬ 
modate  those  flow  rates  cited  herein  and  stressed  to 
withstand  pressures  to  1600  PSI. 

4.  Acceptance  Tests  and  Evaluation 

4.  1  The  contractor  shall  prepare  and  submit  to  the  Navy  a 

comprehensive  program  of  test  procedures  designed  to 
verify  the  performance  requirements  of  section  3. 
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2  The  contractor  shall  then  conduct  tests  to  verify  the 
performance  requirements  of  section  3,  In  addition 
to  contractor-performed  tests  the  equipment  will  be 
examined  and  tested  by  the  Navy. 

3  Test  Instrumentation:  For  the  purpose  of  verifying 
the  output  performance  requirements  a  commercially  - 
available  frequency-type  or  counter-type  meter  may  be 
used. 

4  Fuel  Temperatures  and  Densities:  Fuel  temperature 
shall  vary,  as  required,  over  the  range  specified  in 

3.  1.  1.  Fuel  densities  shall  vary  over  the  range  taken 
by  the  minimum  density  of  JP-4  at  high  temperature 
and  the  maximum  density  taken  by  JP-5  at  low 
temperature. 

5  Accuracy:  Accuracy  tests  shall  first  be  performed  at 
room  temperature.  These  tests  shall  be  repeated  at  the 
specified  fuel  temperature  extremes  with  an  allowable 
variation  from  the  established  room  temperature  error 
of  not  more  than  jt  10%  of  that  error  for  any  position  in 
the  flow  range. 

6  Hysteresis:  The  transmitter  shall  be  subjected  to  a 
series  of  runs  to  verify  the  hysteresis  requirements. 
Room  temperature  accuracy  tests  shall  be  run  at 
increasing  flow  rates  from  minimum  to  maximum. 
Master  transmitter  vs  test  transmitter  readings  shall 
be  taken  at  intermediate  points  moving  up  through  the 
flow  range.  The  flow  shall  be  decreased  and  master 
vs  test  transmitter  readings  again  taken  at  the  same 
flow  po  as  indicated  by  the  master.  Upscale  and 
downsca,  readings  of  the  test  transmitter  shall  be 
within  jt  0.  3%  of  each  other  for  any  flow  rate. 

7  Repeatability:  Not  less  than  ten  room  temperature 
accuracy  tests  shall  be  run,  each  separated  from  the 
previous  by  one  day  in  time.  The  day  to  day  variations 
shall  not  exceed: 

1  PPH  at  extreme  low  end  of  flow  range 
+  3  PPH  at  test  point  of  3,  000  PPH 
15  PPH  at  test  point  of  15,  000  PPH 
+  25  PPH  at  test  point  of  25,  000  PPH 
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Several  intermediate  test  points  in  addition  to  the  above 
shall  also  be  recorded. 

4.  8  Response  Time:  The  response  time  test  shall  be  the 

measurement  of  the  capability  of  the  transmitter  to  track 
a  step  change  in  flow  rate  due  to  the  activation  of  a  quick 
opening  valve. 

4.  9  Vibration:  The  transmitter  shall  ue  tested  for  accuracy 
while  under  vibration.  This  shall  be  done  for  a  selected 
mass  flow  rate  in  the  7,  000  to  15,  000  PPH  range.  The 
vibration  test  shall  be  performed  according  to  MIL-STD-810 
method  514.  1  as  required  by  MIL-E-5007D  para  4.  6.  2.  3.  9 
except  that  the  maximum  temperature  test  shall  be  deleted. 

4.  10  Contaminated  Fuel  Test:  A  room-temperature  accuracy 

test  shall  be  run  contiuously  for  300  hours  at  maximum  flow 
rate  with  fuel  having  contaminants  of  the  size,  nature  and 
quantity  shown  in  Table  X  of  MIL-E-5007D. 
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APPENDIX  B.  LIST  OF  MATHEMATICAL  SYMBOLS 

a  constant 

A  area 

b  constant 

c  sound  speed 

C  polarizability 

f  frequency 

K  meter  factor 

L  axial  projection  of  path  in  flowing  liquid 

M  mass  flow  rate 

N  multiplier 

P  path  length  in  liquid 

Q  volumetric  flow  rate 

Re  Reynolds  number 

t,  T  transit  times 

V  flow  velocity 

V  area  average  flow  velocity 

w  window  (subscript) 

Z  acoustic  impedance 

a  attenuation  coefficient;  thermal  expansion  coefficient 

c  dielectric  constant;  error 

a  standard  deviation 

A  difference 

p  density 

v  kinematic  viscosi  ty 
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APPENDIX  C.  DERIVATION  OF  V  EQUATIONS 


Let  the  zigzag  and  straight  velocimeter  geometries  be  represented 
as  follows: 


Let  P  =  path  length  in  fluid. 

L  =  axial  projection  of  effective  path  in  flowing  fluid. 

t  =  total  time  delay  in  window,  shim,  wedge,  buffer,  impedance 

w 

matcher,  etc.  +  reciprocal  electrical  delays  including  cables 
and  circuit  components,  i.e.,  all  non-fluid  delays. 

t  =  time  difference,  upstream  minus  downstream,  measured  at  zero  flow; 
or  the  difference  between  tj  and  t^  measured  at  steady  (non  zero) 
average  flow  velocity  V  and  constant  sound  speed  c,  for  transposed 
cables,  that  is. 


t 

z 


)  -  (t 

CA  ‘ 


t,)  PS  (*,  -  *?)  n 

Z  1  Z  no  flow 


-  A 


2 


Aj=  electrical  +  nonreciprocal  path  delays  in  configuration  1 
A^-  electrical  +  nonreciprocal  path  delays  in  configuration  2 
tj  =  transit  time  measured  when  A  transmits 
1 2  =  transit  time  measured  when  B  transmits 


(A  transmits). 
(B  transmits). 


If  the  system  were  reciprocal,  and  if  there  were  no  electrical  delays,  then 
at  a  flow  velocity  V, 


where 


2  l.V 
C2 


c  =  2P/(t  +t,  -  2t  ) 
1  2  w 


(1) 

(2) 
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Therefore 


V  = 


2 

c  At 

2L 


t,  -  t. 


/2P\  2  1 

=U  l(« .+  t  -  zrp 


If  the  system  is  not  reciprocal. 


V  = 


2pZ  1»;-V-1.,-AI 


which  if  desired,  may  be  rearranged  as 


V  = 


2P2  ‘2-‘it&rA2 

L 


2t  ) 
w 


If  the  system  is  reciprocal,  A^  =  A^  -  A  and 


V  = 


2P2  t2  '  t\ _ 

L  ftj  +t2-  2A-  2tw)2 


(3) 


(4) 


(6) 


The  transmitter  is  programmed  to  transmit  an  equal  number  of 
times  in  each  configurivtion.  The  times  for  these  transmissions  to  occur 
are  Nt^  and  Nt2-  By  counting  cycles  from  a  10  MHz  clock  during  each 
phase-locked  configuration,  counts  proportional  to  tj  and  t £  are  obtained. 

The  pP  is  programmed  to  manipulate  these  counts  according  to  Eq.  (5).  * 

Provided  V  and  c  are  constant,  resolution  is  improved  by  making  N  large. 

The  upper  limit  on  N  is  imposed  by  whichever  of  the  following  is  smallest: 

(a)  Max  interval  in  which  V  or  c  is  constant; 

(b)  Half  the  required  response  time  (e.  g.  ,  1/2  of  2  sec  in  present  contract); 

(c)  Max  no.  of  counts  that  can  be  accumulated,  y  clock  frequency. 


■f 

If  the  system  turns  out  to  be  reciprocal,  then  Eq.  (6)  would  be  used  auto¬ 
matically.  Note  that  if  cable  lengths  were  100  ft,  2  A  a:  250  ns.  Assuming 
t  +  t-  s  250  ps,  neglecting  the  cable  delays  would  introduce  an  error  in 
V  of  0.2%.  However,  for  cable  lengths  less  than  10  ft,  and  P  >  6  inches, 
A  could  be  neglected,  because  its  contribution  to  errors  in  V  would  be 
<  2A  /tj  =  25  ns/lOOps  =  0.025%. 
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APPENDIX  D.  FLOW  VELOCIMETER  ELECTRONICS 


This  appendix  relates  in  general  to  ultrasonic  flowmeters.  More 
particularly,  it  relates  to  "Model  6000"  circuitry  for  precise  determination 
of  the  differences  in  transit  time  of  ultrasonic  wave  propagation  upstream 
and  downstream  within  a  fluid  flowing  in  a  conduit.*  Jet  propulsion  fuels  such 
as  JP-4,  JP-5  or  their  mixtures,  flowing  in  a  flow  cell  or  spoolpiece,  are 
of  prime  concern. 

One  method  of  determining  the  velocity  of  flow  of  fluids  within  a 
conduit  utilizing  ultrasonic  wave  transducers  has  been  to  place  a  pair  of 
transducers  apart  in  the  direction  of  flow  and  to  determine  the  difference 
in  propagation  time  of  ultrasonic  waves  between  these  transducers  in  the 
upstream  direction  and  the  downstream  direction.  The  flow'  velocity  of  the 
fluid, V,  adds  to  the  sound  velocity  c  in  the  downstream  direction  and  sub¬ 
tracts  from  it  in  the  upstream  direction.  One  effective  technique  for  ultra- 
aonically  measuring  sound  speeds,  and  hence  flow,  is  known  as  the  sing- 
around  method.  In  this  method  an  ultrasonic  wave  is  initiated  at  one  of  the 
transducers  and  is  transmitted  to  the  other  transducer.  Upon  receiving 
this  transmitted  ultrasonic  wave,  the  second  transducer  generates  an  elec¬ 
trical  signal,  which  is  used  to  again  trigger  the  propagation  of  an  ultrasonic 
wave  from  the  first  transducer.  This  process  is  repeated  and  the  repetition 
frequency  of  the  output  signals  from  the  first  transducer  then  corresponds 
to  the  transit  time  of  ultrasonic  waves  passing  from  the  first  transducer 
to  the  second  transducer.  When  the  same  technique  is  used  with  respect 
to  ultrasonic  waves  propagated  in  the  opposite  direction,  a  pair  of  frequencies 
result.  The  difference  between  the  frequencies  Af  is  proportional  to  the 
flow  velocity  of  the  fluid. 

One  problem  witli  the  singaround  technique  is  that  of  rapidly  and 
accurately  measuring  a  small  difference  in  frequency.  Other  problems 
arise  from  reverberation  effects  in  the  conduit  or  fluid  or  ultrasonic  trans¬ 
ducers  themselves,  rendering  it  still  more  difficult  to  obtain  a  correct 
measure  of  Af.  Additionally,  the  presence  of  nonfluid  path  delays  such  as 
delays  due  to  buffer  rods  between  the  transducers  and  the  fluid  leads  to  v 
inaccuracy  in  interpreting  Afin  terms  of  V. 

The  delivered  equipment  avoids  the  above  problems  as  follows. 

Broadly  speaking,  in  the  standard  Panametrics  Model  6000  Ultra¬ 
sonic  Flowmeter,  a  special  phase  locked  loop  (PLL)  techfiique  accurately 
determines  the  flow  velocity.  In  this  approach,  as  in  the  singaround 
approach,  a  pair  of  transducers  are  placed  apart  in  the  direction  of  flow. 

An  ultrasonic  wave  is  initiated  at  one  transducer.  Its  reception  at  the 
other  transducer  locks  the  half-period  of  a  voltage  controlled  oscillator 

:;:Pedersen,  N.  E.  ,  Bradshaw,  J.  E.  ,  Matson,  J.  E.  and  Lynnworth,  L.  C.  , 
Ultrasonic  Flowmeter,  U.  S.  Patent  (1981);  ser.  no.  086402  (Oct.  1979). 
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to  the  transit  time.  The  transit  time  for  waves  transmitted  in  the  opposite 
(upstream)  direction  is  similarly  locked  onto.  The  difference  in  half-periods 
then  provides  for  a  highly  accurate  measure  of  flow  velocity,  capable  of 
compensating  for  nonfluid  delays,  as  well  as  for  reverberations  and  ringing 
within  the  spoolpiece  (transducers,  or  fluid). 


For  a  typical  ultrasonic  flowmeter  configuration  employing  reciprocal 
transducers,  the  flow  velocity  V  along  the  interrogated  fluid  path  P  may  be 
computed  from  the  times  of  flight  tj  and  measured  between  the  two  trans¬ 
ducers: 


V  = 


cZAt 

2L 


(t  +  t  -  2t  )2 
1  2  w 


(1) 


where  L  =  axial  projection  of  path  in  the  flowing  fluid,  and  t  =  sum  of  all 
nonfluid  path  delays  (e.  g.  ,  windows  in  front  of  transducers,  cable  delays, 
electronic  delays). 

Equation  (1)  given  for  V  above  is  an  approximation  that  is  applicable 
to  the  extent  that  <<  cZ. 

The  flow  velocity  V  may  be  converted  to  V,  the  area  averaged  flow 
velocity,  by  utilizing  area  averaging  cells  or  other  means. 

To  understand  how  the  present  PLL  method  differs  from  a  conventional 
singaround  method,  note  that  in  the  Model  6000,  the  transit  time  between 
transducers  and  the  half  period  of  one  cycle  of  the  output  of  a  voltage  controlled 
oscillator  (VCO)  is  compared.  The  resultant  difference  signalis  used  to  control 
the  VCO  until,  after  repeated  ultrasonic  transmissions,  the  frequency  of  the 
VCO  is  such  that  its  period  is  equal  to  twice  the  transit  time  of  the  ultrasonic 
waves. 


In  this  arrangement  it  is  the  output  from  the  VCO  which  can  initiate 
the  repetition  pulses  of  ultrasonic  waves  from  the  transducer.  In  fact,  the 
period  of  the  VCO  need  not  be  equal  to  twice  the  transit  time  between  trans¬ 
ducers,  hut  may  instead  be  made  proportional  to  it  so  that  the  VCO  is  operating 
at  either  a  higher  frequency  (so  that  its  cycle  period  is  a  precise  rational 
fraction  of  the  transit  time  between  transducers)  or  may  be  operated  at  a 
lower  frequency  (so  that  its  cycle  period  is  an  exact  multiple  of  this  transit 
time).  In  this  technique,  however,  it  is  not  the  frequency  which  is  measured, 
but  rather  the  elapsed  time  of  a  number  of  cycles  of  the  VCO.  This  is  accom¬ 
plished  by  generating  at  a  relatively  high  repetition  rate  a  train  of  clock  timing 
pulses  (e.  g.  10  MHz),  at  least  one  order  of  magnitude  higher  than  the  inter¬ 
rogation  frequency,  and  measuring  the  total  accumulated  number  of  these 
timing  pulses  over  the  integer  N  periods  of  the  VCO  output.  This  technique 
substantially  reduce*  the  error  in  determination  of  the  flow  of  velocity,  for  a 
given  clock  frequency  (say  10  MHz). 
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If  response  time  is  unimportant,  the  number  of  periods,  N,  can 
be  selected  to  be  a  convenient,  large  fixed  integer,  say  1024.  This  value 
is  appropriate  for  the  measurement  of  the  flow  of  many  liquids  having  sound 
speeds  ranging  between  300  and  3000  m/s,  or  gases  having  sound  speeds 
from,  say,  100  to  1000  m/s.  For  faster  response,  a  small  N  is  selected, 
e.  g.  ,  32,  64,  or  100.  In  contrast  to  certain  earlier  methods  wherein  N  was 
adjusted  to  compensate  for  variations  in  sound  speed,  the  present  invention 
compensates  for  sound  speed  through  synchronism  (pulse  repetition  period 
proportional  to  time  of  flight)  and  instead  uses  a  generally  fixed  N  as  a  means 
for  optimizing  response  time  and  resolution. 

DESCRIPTION  OF  PREFERRED  EMBODIMENTS 

Figure  Dl  shows  one  form  of  a  flow  cell  with  transducers  4a  and  4b 
coupled  to  buffer  rods  3a  and  3b  which  are  installed  in  threaded  nipples  2a  and 
2b  welded  to  the  pipe  1.  The  transducers  are  cable- connected  to  flowmeter 
electronics  14.  The  transducers  are  cable- connected  to  flowmeter  electronics 
14.  The  pipe  10  has  an  inside  diameter  D,  an  area  A  =  nD^/4,  and  provides 
an  axially-projected  path  length  L  between  the  wetted  ends  5a  and  5b  of  the 
buffer  rods,  and  a  fluid  path  P  along  tilted  diameters.  The  fluid  at  rest  has 
a  sound  speed  c,  and  flows  at  a  velocity  V.  The  pipe  10  may  also  be  pene¬ 
trated  by  thin  rod  waveguides  connected  to  torsional  mode  sensors  6a  and  6b. 

Sensor  6b  is  encased  in  a  hermetically  sealed  sheath  7.  Pipe  wall  penetrations 
are  sealed  by  compression  fittings  8  which  are  of  conventional  design  and  so 
are  not  detailed.  The  purpose  of  the  optional  torsional  sensors  is  to  sense 
fluid  density  and/or  temperature.  However,  they  are  not  discussed  any  further 
in  this  appendix. 

In  Figs.  2a  and  2b,  there  are  shown  several  timing  diagrams.  Those 
in  Fig.  2a  pertain  to  upstream  ultrasonic  interrogation  waveforms,  while  those 
in  Fig.  2b  pertain  to  the  downstream  interrogation.  The  timing  diagrams 
illustrate  the  interrogation  repetition  rate  at  which  ultrasonic  waves  are 
emitted  from  the  transducers  4a  and  4b  as  well  as  timing  pulses  employed  to 
measure  the  period  of  these  interrogation  signals  over  several  cycles.  Turning 
to  Fig.  2a  the  square  wave  of  audio  frequency  f  j  represents  the  repetition 
rate  for  upstream  interrogation,  whose  period  tj  is  to  be  measured.  To 
achieve  high  accuracy,  the  system  employs  a  clock  frequency  f  which  is  in 
the  ultrasonic  range,  i.  e.  ,  f  >  >  f  .  The  number  of  clock  pulses  counted 
during  t^  is  t  f  +  1.  To  reduce  the  fractional  error  due  to  the  +  1  uncertainty, 
f  can  be  increased  or  the  count  can  be  extended  for  a  longer  time  than  t^, 
for  example,  for  Nt  ,  where  N  is  an  integer  greater  than  1,  and  usually  much 
greater  than  1.  In  this  illustration,  by  way  of  example,  N  =  2,  so  that  the 
counting  interval  T^  =  Zty  If  the  same  counting  multiplier  N  is  used  in  down¬ 
stream  interrogation,  then  tj  -  can  be  determined  by  counting  at  the  fc  rate 
for  intervals  T^  and  T^. 

^Figures  Dl-5  cited  in  this  appendix  are  ail  contained  at  the  end  of  this 
appendix,  and  are  not  to  be  confused  with  the  figures  from  the  main  text. 
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By  choosing  a  large  N,  the  +  1  uncertainty  in  the  measured  count 
is  rendered  insignificant.  The  time  required  to  make  such  measurements, 
again,  is  N(t.  +  t^)  or  some  multiple  of  this  time  if  more  than  one  pair  of 
interrogations  upstream  and  downstream  are  averaged.  If  a  given  time  is 
allotted  to  making  a  V  determination,  as  may  be  imposed  by  dynamic  flow 
considerations  and  a  corresponding  short  response  time,  there  are  several 
ways  that  the  +_  1  uncertainty  associated  with  the  single  clock  tick  of  frequency 
f  might  be  reduced.  These  ways  include:  use  a  clock  of  higher  f  ;  use  2 
clocks  of  slightly  different  frequencies  f  and  f  ^  in  a  vernier  mo%e;  or  use 
but  one  clock,  but  supplement  it  with  an  analog  interpolation  circuit. 

In  Fig.  D3  there  is  illustrated  in  block  diagrammatic  form  suitable 
circuitry  for  the  flowmeter  electronics  illustrated  in  Fig.  Dl.  In  the  circuit 
of  Fig.  D3  VCO  20  provides  its  output  to  a  divide  circuit  22,  with  the  divided 
output  from  the  VCO  being  provided  simultaneously  to  a  configuration  switch 
24  and  to  a  one-shot  multivibrator  26,  as  well  as  to  one  input  of  a  phase  detec¬ 
tor  28.  The  output  from  the  one-shot  26  is  provided  through  multiplexer  30 
to  either  of  two  amplifiers  34  and  36.  The  output  of  amplifier  34  is  coupled 
to  the  upstream  transducer  4a  and  the  output  of  amplifier  36  is  coupled  to 
the  downstream  transducer  4b.  The  output  of  the  one-shot  multivibrator  26 
is  also  connected  to  delay  circuit  38  whose  output  feeds  automatic  gain  control 
amplifier  40  in  receiver  41  and  a  comparison  aero  crossing  detector  44.  The 
inputs  to  the  automatic  gain  control  amplifier  40  are  received  from  either 
amplifier  46  coupled  to  the  output  from  transducer  4a,  or  from  amplifier  48, 
coupled  to  the  output  from  transducer  4b.  The  output  from  the  zero  crossing 
detector  becomes  an  input  to  the  phase  detector  28.  The  output  from  the 
one-shot  multivibrator  26  is  provided  as  a  reset  signal  to  the  phase  detector 
28. 


The  output  from  the  phase  detector  28  is  coupled  through  differential 
integrator  50  to  a  multiplexer  52  at  the  control  input  of  the  VCO  20.  Multi¬ 
plexer  52  provides  a  switching  function  to  couple  either  one  of  two  capacitors 
56  or  58  to  the  control  input  of  VCO  20. 

Typically  a  10  MHz  clock  60  is  coupled  through  logic  circuit  62  to 
either  upstream  24  bit  counter  66  or  downstream  24  bit  counter  68.  The 
outputs  from  the  counters  66  and  68  are  coupled  through  a  buffer  stage  70 
to  a  computer  (not  shown). 

The  operation  of  the  circuit  of  Fig.  D3  is  best  understood  in  conjunction 
with  the  timing  diagrams  of  Figs.  D4  and  5.  In  operation  the  VCO  20  provides 
an  output  signal  illustrated  in  Fig.  4  in  which  the  down-going  edge  of  the  wave¬ 
form  after  passing  through  divider  22  (which  for  purposes  of  this  initial  dis¬ 
cussion  will  be  considered  to  have  a  dividing  factor  of  one)  triggers  one  shot 
multivibrator  26  which  produces  the  transmit  pulse.  The  same  waveform 

edge  provides  a  signal  to  a  flip-flop  circuit,  FF  and  FF  ,  resetting  them 
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in  the  phase  detector  28.  The  transmit  pulse  is  passed  to  multiplexer  30, 
which  depending  upon  which  control  signal  it  is  receiving  from  the  con¬ 
figuration  switch  24,  may  pass  that  signal  from  the  one  shot  multivibrator  26 
to  amplifier  34  to  initiate  an  ultrasonic  pulse  from  transducer  4a.  Or,  in  the 
other  mode,  it  passes  that  same  pulse  through  amplifier  36  to  transducer 
4b  to  initiate  an  ultrasonic  pulse  in  the  opposite  direction.  Configuration 
switch  24  operates  to  change  the  state  of  the  entire  circuit  to  operate  in  one 
mode  to  determine  the  time  period  t^  for  upstream  interrogations.  Or,  in 
the  other  mode  it  determines  the  time  period  for  downstream  interrogations. 
The  quantity  by  which  the  divide  circuit  22  divides 'the  output  from  the  VCO 
20  defines  the  number  of  cycles  N  for  which  the  transit  time  is  measured. 

Thus  the  configuration  switch  24  switches  back  and  forth  between  upstream 
and  downstream,  counting  every  N  cycles  from  the  VCO  20. 

The  operation  of  the  overall  circuit  of  Fig.  D3  is  to  control  the  VCO 
such  that  it  is  operating  at  a  frequency  which  is  equal  to  1/2  the  system 
frequency,  where  the  system  frequency  is  defined  as  the  frequency  whose 
period  equals  the  time  of  flight  of  an  ultrasonic  pulse  from  one  transducer  to 
the  other.  Basically,  this  is  accomplished  by  detecting  the  zero  crossings 
of  received  pulses  in  the  zero  crossing  detector  44  and  supplying  them  to  a 
flip-flop  circuit,  FF^  in  the  phase  detector  28.  This  phase  detector  deter¬ 
mines  which  came  first,  the  rising  edge  of  the  FF^  output  or  the  rising  edge 
of  the  received  frequency  FF-,  output.  This  phase  detector  28  then  provides 
an  output  pulse,  whose  width  is  the  difference  in  the  arrival  time  of  the  two 
edges,  to  either  the  plus  or  minus  input  of  the  differential  integrator  50.  It 
is  the  output  of  this  differential  integrator  50  applied  through  multiplexer  52 
which  is  the  frequency- controlling  input  to  the  VCO  20. 

As  the  width  of  the  phase  detector  output  pulse  approaches  zero,  the 
VCO  frequency  will  approach  half  of  the  system  frequency  and  the  circuit  be¬ 
comes  locked  with  the  VCO  output  frequency  tracking  the  system  frequency. 
Since  the  phase  detector  is  always  reset  at  the  same  time  that  the  falling 
edges  of  the  transmit  and  receive  frequencies  are  generated,  it  will  always 
phase  detect  the  succeeding  two  edges  which  are  the  50%  duty  cycle  point  of 
the  transmit  frequency  and  the  zero  crossing  detected  output  of  the  received 
pulse.  This  ensures  that  the  entire  circuit  will  not  lock  on  harmonics  of 
either  the  receive  or  transmit  frequencies  and  that  every  interrogation  of 
the  flow  will  yield  correction  information.  This  feature  provides  for  very 
fast  tracking  of  the  system  frequency  (which  it  will  be  understood  changes 
with  changes  in  the  flow  velocity)  and  enables  non  50%  duty  cycle  frequencies 
to  be  phase  detected. 

The  circuit  operating  in  this  fashion  usually  disregards  transducer 
ringing  signals  because  the  50%  duty  cycle  halves  the  pulse  repetition  fre¬ 
quency,  compared  to  a  traditional  singaround.  This  means  that  the  time 
between  tranmits  is  exactly  twice  the  time  of  flight,  allowing  twice  the  time 
for  the  ringing  to  decay. 
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In  order  to  cancel  the  effects  of  triple  or  multiple  transit  reflections 
of  the  ultrasonic  waves,  or  unusually  long  ringdown  times  (e.  g.  ,  in  buffer 
rods)  the  circuit  may  be  configured  to  skip  a  number  M  of  transmitted  ultr- 
sonic  pulses  to  thereby  allow  enough  time  for  multiple  transits  to  die  out 
while  maintaining  the  VCO  at  exactly  1/2  the  system  frequency.  This  is 
accomplished  by  utilizing  the  divider  22  for  a  relatively  high  number  M  of 
cycles.  Thus,  after  the  50%  duty  cycle  edge  of  the  transmit  frequency  and 
the  received  edge  of  the  received  frequency  have  been  compared,  another 
transmit  pulse  to  the  transducer  is  not  initiated  and  therefore  the  flip-flops 
FFj  and  FF^  are  not  reset  until  M  cycles  have  passed.  This  is  allowable 
because  the  phase  detector  will  only  operate  on  the  two  succeeding  edges 
following  the  reset  pulse.  See  Fig.  D5.  Dividing  by  M  >  2  is  called  a  skip  mode. 

For  longer  ultrasonic  path  lengths  the  time  between  the  detected 
pulses  increases  and  storage  of  the  corrected  VCO  input  voltages  is  mandated. 
This  is  accomplished  by  the  differential  integrator  50  which  uses  high  im¬ 
pedance  field  effect  transistor  switches  to  keep  capacitor  leakage  low.  The 
storage  capacitors  56  and  58  are  alternatively  switched  into  the  circuit, 
depending  upon  whether  the  measurement  mode  is  upstream  or  downstream, 
as  will  be  explained  below. 

With  increased  conduit  diameters,  the  effects  of  reflections  and 
ringing  become  small  enough  due  to  beam  spreading  attenuation  or  sound 
absorption  so  that  operation  in  the  skip  mode  is  not  required.  In  this  case, 
upon  receiving  the  ultrasonic  signal  (zero  crossing),  a  new  transmit  can  be 
immediately  initiated.  At  the  same  time  a  counter  is  incremented  so  that 
a  total  of  N  increments  is  used  to  acumulate  N  samples  of  the  time  of  flight, 
whereas  in  the  PLL  incrementing  a  counter  on  each  receive  for  N  increments 
will  yield  2N  samples  (in  twice  the  time  for  equal  times  of  flight)  of  the  time 
of  flight.  The  number  N  is  implemented  by  selecting  a  fixed  divided  output 
of  divider  22.  By  dividing  the  VCO  frequency  directly,  the  regular  mode  and 
the  skip  mode  will  both  yield  the  same  number  of  times  of  flight  for  a  given 
divisor. 


As  previously  indicated,  the  output  of  the  divide  circuit  22  is  used  to 
toggle  the  conf  iguration  switch  24,  changing  the  overall  circuit  from  an  up¬ 
stream  counting  mode  to  a  downstream  counting  mode.  The  output  of  the 
configuration  switch  24  is  indicated  in  one  mode  as  a  straight  arrow  and  in 
the  second  as  an  arrow  with  a  small  circle.  Throughout  the  circuitry  this 
convention  is  applied  so  that  the  configuration  of  the  multiplexers  30  and  52, 
as  well  as  that  of  the  logic  circuit  6  2  and  the  gating  of  amplifiers  46  and  48 
are  all  shown  as  controlled  by  the  toggle  output  from  the  configuration  switch 
24.  The  output  of  the  configuration  switch  24  applied  to  the  logic  circuit  62 
operates  to  gate  the  oscillator  clock  60  sway  from  the  upstream  24  bit  counter 
66  to  the  downstream  24  bit  counter  68.  The  numbers  in  these  counters  can 
then  be  used  by  a  computer  to  compute  the  velocity  of  flow  V.  By  switching 
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the  multiplexer  52,  the  configuration  switch  24  switches  the  capacitor  used 
on  the  output  from  the  differential  integrator  50  and  therefore  allows  one 
set  of  capacitors  (each  of  the  capacitors  56  and  58  is  actually  a  pair  of 
capacitors)  to  hold  the  control  voltage  for  the  upstream  operation,  and  the 
other  set  of  capacitors  58  is  used  to  hold  the  control  voltage  for  the  down¬ 
stream  operation. 

Any  errors  associated  with  the  switching  in  and  out  of  these  capacitors 
and  the  subsequent  need  of  a  number  of  cycles  to  recharge  the  capacitors  to 
the  correct  value  are  eliminated  as  follows.  The  Model  6000  only  gates  on 
the  oscillator  clock  60  to  the  counters  for  a  particular  number  of  counts  at 
the  end  of  a  number  of  cycles  of  the  VCO.  In  fact,  any  number  of  cycles  of 
transmits  can  be  ignored  by  these  counters.  The  computer  will  divide  by 
only  the  number  of  cycles  actually  used  to  gate  on  the  counters.  Such  a 
method  allows  the  circuit  to  well  establish  its  operating  frequency  before  the 
counters  are  employed. 

The  S  P  electronics  produces  a  voltage  which  is  linearly  related  to  p 
in  the  range  of  1.  5  to  6.  5  volts.  The  output  of  the  S  P  electronics  is  con¬ 
nected  to  a  Teledyne  analog  to  digital  converter.  The  A/D  converter  upon 
command  of  microprocessor  performs  a  conversion  of  input  voltage  to  digital 
format  and  then  signals  the  computer  that  the  conversion  is  complete.  The 
computer  commands  and  receives  a  conversion  for  every  flow  sample.  Using 
the  formula  for  the  S  P  output,  the  computer  converts  the  digital  code  of  the 
■A/  D  first  to  a  voltage  and  then  to  a  value  of  density  in  lb/  gal  (this  is  part  of 
its  program).  The  computer  also  has  calculated  the  flow  velocity  from  the 
latest  upstream  and  downstream  transit  times  (tj,  t^) -  It  then  multiplies  the 
velocity  times  the  cross  sectional  area  of  the  cell  (a  number  that  is  stored  in 
its  permanent  memory),  and  finally  multiplies  that  result  by  the  density 
obtained  from  the  conversion.  The  data  point  is  then  output  to  the  output  devices. 

The  logic  circuit  referred  to  on  p.  D-4  consists  of  three  flip  flops  and 
a  few  logic  gates.  Their  purpose  is  to  steer  the  10  MHz  clock  to  one  counter 
or  the  other  depending  upon  the  state  of  one  of  the  flip  flops,  whose  state  is 
changed  every  time  the  configuration  switch  is  changed.  The  other  two  flip 
flops  arc  used  to  signal  the  computer  that  a  count  has  been  completed,  so  that 
the  computer  can:  1)  read  which  counter  is  ready  (the  state  of  the  first  flip 
flop);  2)  read  the  counter;  and  3)  reset  the  counter  for  the  next  interval. 
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APPENDIX  F:  EXTRACT  FROM  NASA  RFP 


EXHIBIT  "B* 

DESIGN  GUIDELINES  AND  SPECIFICATIONS 

1.  Measurand  and  method  of  iMiurentnt  -  Fuel  mass 
flow  (both  mass  flow  and  total  mass  flow  accumulated) 
shall  be  measured  either  directly,  or  by  using  a  composite 
system  of  separately  measuring  volumetric  flow  and  f uel 
density,  or  by  measuring  a  combination  of  quantities  from 
which  mass  flow  can  be  calculated.  For  a  flowmeter  system 
that  measures  volumetric  flow,  the  preferred  method  of 
measuring  density  is  to  use  a  densitometer  rather  than  a 
correlation  between  density  and  temperature.  The  reasc- 
is  that  fuels  that  are  within  the  MIL  specifications  car 
still  show  variations  between  batches  that  can  produce  un¬ 
certainties  up  to  ± 1/2  percent  in  density  when  determinec 
by  a  temperature  measurement.  This  uncertainty  exceeds 
the  contract  goal. 

2.  Types  of  fuels  -  Fuels  of  interest  are  JP-3,  JP-4, 
JP-5,  JP-8  (Type  A-l)  and  Type  A. 

3.  Flow  range  -  Typical  ranges  of  fuel  flow  between 
engine  full  power  and  idle  are  between  50:1  to  100:1  de¬ 
pending  on  whether  or  not  the  engine  has  an  afterburnoi  . 
The  absolute  value  of  full  scale  flowrate  varies  with 
engine  thrust?  but  for  the  purpose  of  this  contract,  a 
flowmeter  with  a  nominal  full  scale  of  6.3  1/s  (40,noo 
lb/hr)  and  a  50:1  operating  range  is  of  primary  interest. 

A  full  scale  range  of  0.5  1/s  (3000  lb/hr)  is  also  of 
interest  and  shall  be  considered. 

4.  Pressure  -  Flowmeters  are  subjected  to  high  pres¬ 
sures  because  tEey  are  usually  located  downstream  of  the 
fuel  pump.  Operating  pressures  up  to  7000  kPa  (1000  psi) 
shall  be  considered.  Flowmeter  bodies  shall  be  hydro¬ 
static  pressure  tested  to  1.5  times  the  maximum  operating 
pressure. 

5.  Pressure  drop  -  At  maximum  fuel  flow  the  maximum 
pressure  drop  across  the  flowmeter  shall  be  68  kPa  (10 
psi) . 

6.  Fuel  temperature  -  The  flowmeter  shall  be  capable 
of  operating  over  a  fuel  temperature  span  from  -55°C  to 
130°C . 

7.  Ambient  temperature  -  The  flowmeter  shall  be  capa¬ 
ble  of  operating  over  an  ambient  temperature  span  from 
-55°C  to  130°C. 

8.  Accuracy  -  The  total  error  band  for  mass  flow  mea- 
surement  shall  be  no  greater  than  0.25  percent  of  read¬ 
ing.  The  project  goal  of  the  error  band  is  0.1  percent  of 
reading. 
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9.  Resolution  -  The  maximum  value  of  resolution  of 
fuel  flow  measurement  at  1/SO  of  full  scale  shall  be  n.2^ 
percent. 

10.  Ambient  pressure  -  The  flowmeter  and  signal  condi¬ 
tioning  electronics  are  usually  located  in  unpressurized 
regions  of  the  aircraft.  The  units  shall  be  capable  of 
satisfactory  operation  in  an  external  pressure  environment 
between  100  kPa  and  7  Pa  (15  to  0.001  psi). 

11.  Vibration  Characteristics  -  The  flowmeter  shall  be 
capable  of  satisfactory  operation  in  the  following  vibra¬ 
tion  -  frequency  envelopes  ±1.2  mm  (5  to  14  Hz), 

±1  g  (14  to  23  Hz),  ±0.45  mm  (23  to  90  Hz),  and  ±15  g  (90 
to  2  kHz)  . 

12.  Size  and  Weight  -  Because  of  space  limitations 

associated  with  flight  applications,  the  flowmeter  shall 
fit  within  the  cylindrical  envelope  outlined  by  the  AN 
fitting  nuts  (MS-33656)  associated  with  the  nominal  fuel 
line  tube  size  for  the  flow  range  specified  in  paragraph  3 
of  Exhibit  "B".  A  protrusion  from  the  side  of  the 
flowmeter  is  acceptable  but  shall  be  no  larger  than  a 
volume  of  the  following  dimensions:  3  fuel  line  tube 

diameters  long,  2  fuel  line  tube  diameters  high,  and  3/4 
tube  diameters  wide.  The  maximum  length  of  the  flowmeter 
including  end  connections  (see  paragraph  19  of  Exhibit 
"B")  shall  be  8  fuel  line  tube  diameters.  The  maximum 
size  of  the  signal  conditioning,  which  may  be  located 
remotely  from  the  flowmeter,  shall  be  1000  cm3.  The 
maximum  weight  of  the  flowmeter  assembly  including  any 
required  valves  and  manifolds  (but  not  including  signal 
conditioning)  shall  be  5  kg. 

13.  Material  -  Wetted  parts  of  the  flowmeter  shall  be 
compatible  with  the  fuels  listed  above  and  at  the  pres¬ 
sures  and  temperatures  listed  above  without  suffering  cor¬ 
rosion,  brittleness,  seal  leakage,  or  other  degrading  pro¬ 
perties. 

14.  Response  -  Time  constant  of  the  flowmeter  shall 
not  exceed  25  ms. 

15.  Failure  mode  -  Because  the  safety  of  the  aircraft 
is  of  paramount  importance,  any  failure  of  the  flowmeter 
shall  not  cut  off  fuel  supply  or  otherwise  interfere  with 
proper  engine  performance. 

16.  Power  -  The  flowmeter  and  signal  conditioning  (if 
requiredl  shall  operate  on  28  V  dc. 

17.  Output  -  The  output  (or  outputs)  of  the  flowmeter 
including  signal  conditioning  (if  required)  shall  be  a 
voltage  or  frequency  which  is  a  single  valued  function  of 
fuel  mass  flow  or  of  quantities  from  which  flow  can  be 
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computed.  The  output  (or  outputs)  shall  be  compatible 
with  digital  processing  techniques.  The  data  processing 
of  the  output  signal  (or  signals)  is  not  considered  part 
of  this  contract. 

18.  pressure  pulsations  -  The  flowmeter  performance 
shall  be  unaffected  by  maximum  fuel  line  pressure  fluctua¬ 
tions  of  ±2  percent  of  the  fuel  pressure  for  frequencies 
of  fluctuations  above  10  Hz. 

19.  Mounting  and  position  sensitivity  -  The  flowmeter 
shall  have  AN  Series  s7°  male  flared  tube  (MS-33656)  end 
connections.  Flowmeter  performance  shall  be  unaffected  by 
changes  in  operating  attitude. 

20.  Overrange  capability  -  The  flowmeter  performance 
shall  be  unaffected  after  being  subjected  to  a  fuel  flow 
of  125  percent  of  full  scale. 

21.  Calibration  -'it  is  likely  that  insitu  calibration 
of  the  flowmeter  in  the  aircraft  will  not  be  done  because 
of  the  complexity  of  such  a  procedure.  However,  the 
flowmeter,  including  identical  flight  fuel  system  upstream 
and  downstream  tubing  sections,  shall  be  calibrated  as  an 
assembly  on  a  flow  stand  of  sufficient  precision  to 
determine  flowmeter  accuracy. 
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APPENDIX  G.  CLAMP-ON  AND  OFFSET  ALTERNATIVES  FOR  THE 
FLOW  VELOCIMETER 

The  V  portion  of  the  electronics  may  be  easily  adapted  to  flow 
cell  configurations  other  than  the  delivered  square-holed  sleeve.  As 
mentioned  in  the  body  of  this  report  in  connection  with  Figs.  9,  10  and  12, 
two  such  alternatives  are  the  clamp-on  and  the  offset.  The  main  features 
of  these  alternatives  are  compared  with  the  square-holed  sleeve  in 
Table  Gl. 


Table  Gl. 

Comparison  of  Three  Flow  Cell  Designs. 

Tyre 

Main  Advantages 

Main  Disadvantages 

Clamp-On 

Ease  of  installation  on  Inaccuracy  due  to  profile 

existing  pipe;  ease  of  uncertainty  and  wall- 

maintenance;  no  cost  borne  noise;  L  usually 

for  cell  (except  for  clamp),  small  when  D  is  small. 

Offset 

100%  of  flow  is  area- 
averaged  for6<D^ 

30  mm;  resolution 
improved  in  proportion 
to  L;  L  a  P  -  D;  low 
cost. 

APi  end  effects  due  to 
offset  geometry;  in¬ 
stallation  must  avoid  gas 
entrapment. 

Square-Holed 

Sleeve 

100%  of  flow  is  area- 
averaged  without  need 
for  offset. 

Expensive;  P  >  \[2.  L; 
installation  must  avoid 
gas  entrapment. 

Clamp-On  (Externally  Mounted)  Transducers 

Externally- mounted  transducers  do  not  penetrate  the  pipe  wall. 
These  wedge-type  angle  beam  transducers.  Fig.  Gl,  typically  provide  low- 
velocity  (c.)  oblique  incidence  waves  (e.  g.  ,  vertically  polarized  shear  at 
an  angle  oF  incidence  9 ^  =  60°)  which  are  coupled  to  the  pipe's  exterior. 

If  the  pipe  wall  thickness  W  is  large  compared  to  the  wavelength  A  therein 
the  preferred  wave  in  the  wall  is  also  shear,  and  of  velocity  C2*  Ii  the 
wedge  and  pipe  are  made  of  the  same  material,  c^=  c^. 

If  the  wedge  and  wall  are  made  of  elastically  dissimilar  materials, 
and  in  particular  if  the  wedge  is  a  relatively  low-speed  material  compared 
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to  the  pipe,  and  further,  if  the  pipe  wall  is  not  too  thick,  (e,  g.  ,  W 

comparable  to  ^21  a  a  c*uas Lamb  wave  may  be  launched  in  the 

wall  at  the  refractecTangle  0^  =  90°.  In  other  words,  c./sin0.  =  C2/sin90o 

The  analysis  of  this  problem  and  the  implementation  of  its  solution  are 

simplest  if  c.  is  not  a  function  of  temperature.  Hence,  an  "isopaustic" 

wedge  is  preferred  in  which  c,  <  c„,  ,  for  the  pipe  and  low  -speed 

T  ,  ,  ,  .  ,  .1  2 Lamb 

Lamb  wave  mode  of  interest. 

In  the  delivered  system,  the  V  portion  of  the  electronics  may  be 
modified  for  use  with  clamp-on  transducers  of  the  type  described  above. 
One  simple  form  of  this  solution  is  adequate  only  to  the  extent  that  the 
profile  is  developed,  L  =  constant  and  \>  -  constant.  This  v  constraint 
means  Re  is  a  function  of  V  only,  and  so  the  meter  factor  K  can  be 
programmed  as  follows:  K,aminar  =  0.  7500,  Ktransltlonai  =  0.84, 

^turbulent  1  '»•  01 1  >“8  R'»' 

Derivation  of  clamp-on  equations  corresponding  to  Lamb  or  SV 
waves  in  the  wall  follow. 


Table  G2. 

Typical  values  for  clamp-on 

parameters  when  obliquely 

incident  wave  is  a  vertically  polarized  shear  wave,  for 
Lamb  and  SV  waves  in  the  wall  and  notation. 

Parameter 

Lamb  Wave  in  Wall 

SV  Wave  in  Wall 

0,.  deg 

0,.  deg 

60 

60 

90 

60 

Qy  deg 

45  to  60 

20  to  30 

c  ,  m/s 
cz,  m/s 

1500  to  2000 

3000 

~  2000 

~  2500 

c^,  m / s 

1000  to  2000 

1000  to  2000 

tj.HS 

t2’ 

t3,ps 

20 

20 

(L-S)/c 

(L-S)/c 

P/C3Z 

P/C3 

P 

D/cos0^ 

D/cos0^ 

L 

D  tan0^ 

D  tan0^ 

W 

Comparable  to  \ 

Large  compared  to  \ 

'These  derivations  are  based  on  models,  Fxg.Gl,  in  which  a  multi¬ 
plicity  of  equal-delay  paths  exist  between  the  transducers,  such  that 
constructive  reinforcement  occurs  for  a  range  of  transducer  spacings. 
The  equal-delay  concept  is  based  in  part  on  observations  reported  by 
L.  C.  Lynnworth  in  the  1979  Ultrasonics  Symp.  Proc.  ,  pp.  376-379, 
IEEE  Cat.  No.  79CH1482-9SU  (1979). 
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Table  G2.  (cont'd) 

Special  Notation  for  Clamp-On  Derivations 

t 

average  transit  time,  or  one-way  transit  time  at  no¬ 
flow 

c 

z 

axial  projection  of  sound  speed  in  pipe  wall 

delay  outside  pipe 

tz 

delay  along  pipe  wall 

*3 

delay  in  fluid 

S 

spacing  between  intersections  of  transducer  center- 
lines  and  pipe  exterior 

V9i 

incident  sound  speed  and  angle  of  incidence 

c3.#3 

refracted  sound  speed  and  angle  of  refraction  in  fluid 

D,  L,  P,  W, 

z  dimensions  as  shown  in  illustration.  Fig.  Gl. 

Lamb  waves  in  wall 

Snell's  Law  states: 


Cj/sinflj  =  c2/sin#2  =  c3/8in<?3  (Gl) 

The  refracted  wave  in  the  wall  propagates  at  the  velocity  c  =  c  /sin#  . 

L  1 

Assuming  t^  is  measured  by  butting  the  probes  (angle  beam  trans¬ 
ducers)  with  the  transducer  centerlines  aligned,  and  with  F  measured  as 
the  average  transit  time  f(tUp  +  tcjown) /zl  ,  one  may  write 


Now 


and 


l  = 


V*2 


S  -  D  tan  6 
c^/sintfj 
D  sin#^ 

Cj  sin cos  #3 


(G2) 

(G3) 

(G4) 


Therefore 


t-t 


1 


(sin#^)  (S-D  tan#3) 
c 


D  sin#j 
c  sin 8  cos  0 


(G5) 
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Thus, 


c ,  )  g  j 

-  Dsi^  •  5  =  .rn53c^e3  -  U“e3  =  CO,<?3 


0  =  COt 
3 


'  fe)(^  -  •) 


where 


V  =  c3  At/2L 


c  =  (c  /sin0  )(sin0  ),  L  =  Dtan0  ,  At  =  t  -  t 

3  113  3  upstream  downstream 


an<^  V  =  KV  where  K  =  meter  factor.  The  volumetric  flow  rate  is  Q  =  VA  = 
ttD  v /4.  The  mass  flow  rate  is  M  =  pQ,  where  p  =  density. 

SV  Waves  in  Wall 

This  derivation  is  similar  to  the  Lamb  wave  case,  the  principal 
difference  being  the  "refracted"  wave  in  the  relatively  thick  wall  zigzags 
at  an  effective  velocity  in  the  z  direction  given  by  c  =  c^  sin0^,  assuming 
for  simplicity  that  wedge  and  wall  have  identical  sfiear  wave  velocities. 


As  before,  one  may  write  for  Fig.  Gl  (c). 


t  =  tl  +  t2  +  t3 
D  sin  0 

t  = - - - 

3  Cj  sin 6 ^  cos  0 ^ 


(G2-) 


(C4>) 


but  now 


_  S-L-x 
t2  “ c  sin 0 


2  W  _  S-L-x  2W 


2SinW2  c2sin97  CjSinS^  c^  sin0j 


whe  re 


x  =  S-  L-  2W  tan  9 ^  and  L  =  D  tanS^ 


Therefore  S-Dtan0  -x  Dsin0  x  sin0 

£  —  £  —  "  - ■ —  "  ’  11  -f 

1  c  j  sin0j  c  ^  sin  0^  cos  9^  c^ 


(G51  ) 


c  sin0  w 

- D - “  Vd  sin20J  =  sin20i  cot03+  —  sin  20 } 


(G6-) 
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or 


D  sin 9 


(t  -  t, )  -  —  =  cot  0,  + 


1 


2  W 
D 


cot  8. 


(G6") 


After  solving  this  equation  for  9^,  one  proceeds  to  calculate  c^,L,  V, 
V,  Q  and  M  as  before. 

Offset  Flow  Cell 


The  delivered  transducers,  or  others  of  similar  design,  may  be 
installed  in  tees  connected  at  the  ends  of  a  section  of  pipe  or  tubing  (Fig,  12). 
The  effective  interrogation  length  L  is  very  nearly  proportional  to  the  length 
of  the  sectionbetween the  tees.  Hence,  this  cell  design  offers  an  economical 
way  to  achieve  high  resolution.  For  example,  in  the  delivered  system, 

L  10  cm.  An  order-of-magnitude  increase  in  L  is  achieved  with  a  1-m 
long  pipe  section.  The  repeatability  in  V  corresponding  to  L  =  1  m  may  be 
determined  from  the  nomogram,  Fig.  21,  or  from  the  equation 

AV  =  c2At/2L 

where  we  take  the  electronic  repeatability  to  be  At  =  1  ns  and,  for  example, 
c  =  1500  m/s.  Hence 

AV  =  (2.  25-  106)  (10  ^)/2  as  1  mm/s. 


If  this  offset  design  were  to  be  utilized  with  the  delivered  system, 
the  mainelectronic  change  would  be  in  the  constants  prom,  in  which  the 
new  L,  P  and  A  (area)  would  have  to  be  entered.  The  upper  limit  on  L 
is  imposed  either  by  space  constraints  or  by  attenuation. 

For  the  offset  axial  flow  cell,  or  for  any  other  conduit  of  circular 
cross  section,  the  volumetric  flow  rate  Q  may  be  graphed  as  a  function  of 
duct  inside  diameter  D  and  area  averaged  flow  velocity  V  as  in  Fig.  G2. 

For  small  engine  applications  or  other  situations  where  there  is 
not  enough  axial  length  between  existing  ports  to  install  a  standard  axial 
offset  flow  cell,  a  folded  path  flow  cell  may  be  considered.  By  routing 
the  flow  oblique  or  orthogonal  to  the  initial  flow  direction,  opportunities 
may  be  created  to  achieve  a  large  enough  L  to  obtain  the  required  At,  and 
thereby  obtain  the  required  sensitivity  to  flow.  Two  examples  of  folded 
path  offset  flow  cells  are  shown  in  Fig.  G3. 
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Fig.  Gl.  Models  for  derivations,  (a)  Obliquely  incident  broadband 
SV  pulse  launches  quasi-Lamb  wave  in  wall.  Lamb  wave  leaks  across 
fluid,  generating  a  series  of  paths  having  equal  delays,  (b)  Equivalent 
(simplified)  schematic  for  Lamb  waves  in  wall,  (c)  Diagram  analo¬ 
gous  to  (b)  except  W  is  larger  than  pulse  width,  preventing  generation 
of  Lamb  waves.  Instead,  SV  pulse  zigzags  along  wall.  One  of  the 
equal-delay  paths  is  shown. 
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Fig.  G2.  Volumetric  flow  rate  Q  in  gpm  and  liters /s  as  a  function  of  D  and  V. 

It  may  be  of  interest  to  note  that  V  (ft/s)  ^  0.4Q/D^  for  D  in  inches 
and  Q  in  gpm.  If  Q  is  in  liters /second  and  D  is  in  cm,  V  =  0.4Q/irD^ 
m/s  =  1.27Q/D^  m/s. 
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Fig.  G3.  Folded  path  offset  flowmeter  configurations  enable  flow  to  be  inter¬ 
rogated  ultrasonically  in  at  least  two  noncolinear  conduit  segments,  (a)  Sym¬ 
metrical  design:  axes  of  all  conduit  segments  lie  in  one  plane;  adjacent  segments 
are  orthogonal;  two  main  parts  of  folded  path  are  parallel,  spaced  apart  a  dis¬ 
tance  X,  of  equal  length  (Yj  =  Y^  >>  X)  and  orthogonal  to  aligned  inlet  and  outlet 
axes,  (b)  As  in  (a)  except:  axes  not  all  in  same  plane;  inlet  and  outlet  ports  are 
"reversed"  and  misaligned  slightly,  enabling  the  distance  between  ports  Xp  to  be 
as  small  as  the  conduit  inside  diameter  D.  (c)  Unsymmetrical,  nonplanar  design 
uses  skewed  segments  of  unequal  length  (X  +  Yj  #  Y^)  to  adapt  to  arbitrary 
configuration  of  inlet  and  outlet  ports.  In  the  folded  path  configurations  shown, 
the  effective  interrogation  path  length  L  exceeds  the  distance  between  inlet  and 
outlet  ports.  n  a 
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rig.  Ob.  Chart  recording  showing  jitter  as  a  function  of  flow  velocity  V  and  as  a  function  of  the  RC  time 
constant  T  of  the  analog  output  (I  =  4  to  20  mA)  of  a  Model  6000  ultrasonic  flowmeter.  The  cell  used  in  this 
test  was  the  axial  offset  type,  using  the  AH  flow  loop  of  Fig.  28,  p.  73. 
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Table  G3.  Jitter  observed  in  offset  flow  cell  of  Fig. 28,  measured  by 
Model  6000  ultrasonic  flow  velocimeter.  Date  of  test: 

16  April  1981. 

VELOCITY 

OF  FLOW, 
FT/S 

5.  89 

4.  28 

1.  88 

0 

TIME 

CONSTANT, 

SECONDS 

1.  97 

0.47 

0.47 

0.47 

0.47 

ANALOG  OUTPUT  CURRENT, 

I,  milliamperes 

_ 

7.  15 

7.  16 

7.  15 

7.  15 

7.  15 

7.  14 

7.  13 

7.  12 

7.  J  3 

7.  12 

7.  18 

7.  17 

7.  18 

7.  16 

7.  17 

7.  14 

7.  13 

7.  13 

7.  12 

7.  11 

6.  28 

6.29 

6.  30 

6.29 

6.28 

6.29 

6.  28 

6.27 

6.  30 

6.28 

5.01 

5.  01 

5.  01 

5.  01 

5.  01 

5.  01 

5.01 

5.01 

5.  01 

5.  01 

4.  06 

4.  06 

4.  06 

4.  06 

4.  06 

4. 06  j 
4.06  j 

4.  06 

4.  06 

4.  06 

average 

OF  10 
READINGS 

7.  14 

7.  15 

6.  286 

5.  01 

4.  06 

STANDARD 
DEVIATION  (a) 
OF  ABOVE 

10  READINGS 

.014 

.  026 

. 00966 

0 

0 

~  1  % 

0.  20 

0.  36 

0.  15 

0 

0 

APPENDIX  H.  HYDROSTATIC  PRESSURE  TEST  REPORT 


INSPECTION  REPORT 


SUMMARY: 


H-l 


NEXT  OPERATION: 


NADC-80254-60 


APPENDIX  L  TEMPERATURE  TEST  DATA 


Date 

Transducer 
Temp.  ,  C 

Echo  Amplitude 
mV 

Oscillogram 
I  - 

26  Aug.  1980 

23 

1580 

1 

-40 

1470 

2 

+150 

1480 

3 

+175 

1500 

4 

2  Sept.  1980 

25 

1500 

5 

-196 

1400 

6 

-196 

700 

7 

-196 

750 

8 

3  Sept.  1980 

-100 

850 

9 

-83 

1100 

10 

-55 

1200 

11 

Notes:  -196°C  achieved  in  liquid  nitrogen 

-100°C  achieved  after  liquid  nitrogen  evaporated 

High  temperature  achieved  in  environmental  chamber 
o  o  .  o 

-40  C,  -55  C,  -83  C  in  environmental  chamber 


Remarks 

Start  test 
-0.  6  dB 


Resume  test 

0.  5  hr  dwell 
2.  5  hr  dwell 


>  10  hr  dwell 
at  T  <  -100°C 

-3  dB 
-2  dB 
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Report  of  Test  on 

ZIG-ZAG  FLOW  CELL 
SINE  VIBRATION  TESTING 
for 

PANAMETRICS 


Associated  Testing  Laboratories,  Inc. 

Wayne,  New  Jersey  07470 
Burlington,  Massachusetts  01803 

Date _ September  18,  1980 


Prepared 

Checked 

Approved 

By 

R.  Montvitt 

E.  E.  Kulcsar 

Signed 

Date 
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Administrative  Data 


1.0 

Purpose  of  Test: 

To  subject  the  Flow  Cell  to  sine  vibration  testing  in  accordance 
with  MIL-STD-8 10C  and  the  procedures  of  this  report. 

2.0 

Manufacturer:  Panametrics 

221  Crescent  Street 

Waltham,  MA  02154 

3.0 

Manufacturer's  Type  of  Model  No.: 

Zig-Zag  Flow  Cell:  P/N 

N6  2  26  9- 7  8-C-0069 . 

4.0 

Drawing,  Specification  or  Exhibit: 

MIL-STD-810C ,  Procedure  I, 

Part  1,  Cat.  B  Equipment, 

Curve  L. 

5.0 

Quantity  of  Items  Tested: 

One  (1)  (S/N  001) 

6.0 

Security  Classifcation  of  Items: 

Unclassified 

7.0 

Date  Test  Completed: 

September  18,  1980 

8.0 

Test  Conducted  By.  Associated  Testing  Laboratories,  Inc. 

9.0 

Disposition  of  Specimens:  Returned 

to  Panametrics. 

10.0 

Abstract: 

The  Flow  Cell  was  subjected  to  sine  vibration  testing  as 
outlined  in  this  report.  There  was  no  evidence  of 
physical  damage  as  a  result  of  the  testing. 


Report  No - _X~2248_~11 


Page 
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Burlington,  Massachusetts  01803 
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VIBRATION  TEST 


TEST  PROCEDURE 


The  Zig-Zag  Flow  Cell  was  subjected  to  a  vibration  test  in  accor¬ 
dance  with  MIL-STD-810C,  Method  514.2,  Procedure  I,  Part  1,  Cat. 

B  Equipment,  Curve  L,  as  follows. 

The  unit  was  securely  attached  to  a  vibration  test  fixture  which, 
in  turn,  was  securely  attached  to  the  table  of  a  vibration  exci¬ 
ter.  A  control  accelerometer  was  mounted  to  the  fixture  for  con¬ 
trolling  the  input  vibration  amplitude  of  the  unit.  A  monitor 
accelerometer  was  mounted  to  the  unit  for  monitoring  the  vibration 
amplitude  on  the  unit.  The  unit  was  then  subjected  to  a  resonance 
search  over  the  frequency  range  of  5  to  2000  Hz  at  the  levels  in 
Table  I  below. 

The  unit  was  then  subjected  to  30  minute  resonance  dwells  at  the 
four  most  severe  resonant  frequencies  which  were  determined  using 
the  above  procedure.  The  resonance  dwells  were  performed  at  the 
levels  in  Table  I. 

At  the  conclusion  of  the  resonance  dwells  in  each  axis,  the  unit 
was  subjected  to  vibration  cycling  over  the  frequency  range  of 
5  to  2000  Hz  at  the  levels  in  Table  I.  The  frequency  range  of  5 
to  2000  Hz  and  return  to  5  Hz  was  traversed  during  a  time  inter¬ 
val  of  20  minutes.  The  unit  was  subjected  to  a  total  of  three 
hours  of  vibration  per  axis,  including  resonance  dwell  and  cycl¬ 
ing  time. 


Table  I 


Frequency  (Hz) 

5-14 
14  -  23 
23  -  104 
104  -  2000 


Amplitude 

0 . 1  inch  d . a . 
l.Og's 

0.036  inch  d.a. 
20.0g's 


The  above  procedure  was  performed  in  each  of  three  mutually  per¬ 
pendicular  axes. 
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VIBRATION  TEST 


TEST  RESULTS 


There  was  no  evidence  of  physical  damage  noted  as  a  result  of 
the  above  testing.  During  the  resonance  search,  the  following 
resonant  frequencies  were  determined:  X  Axis  1950;  Y  Axis 
1425  Hz;  Z  Axis  1530  Hz  and  1900  HZ.  A  resonance  dwell  was 
performed  at  the  above  mentioned  frequencies  for  a  period  of  30 
minutes.  The  Flow  Cell  was  returned  to  Panametrics  for  further 
post  test  evaluation.  X-Y  plots  generated  from  the  above  test¬ 
ing  may  be  referred  to  in  the  appendix  of  this  report. 

Explanatory  note  added  April  1981  by  Panametrics  concerning 
strip  chart  recordings  on  pages  J-l6,  J -  1 7 : 

The  deflections  recorded  on  pages  J-16  and  J-17  represent  ana¬ 
log  signals  proportional  to  flow  velocity  V,  as  explained  on  p.  34. 
The  excursions  for  the  cell  at  rest  are  typically  +_  1  minor  division 
For  the  cell  tinder  vibration,  the  same  magnitude  is  observed 
except  for  the  group  of  resonances  observed  near  the  graduation 
marked  2000  Hz,  during  each  logarithmic  sweep  of  frequency. 

At  the  resonances,  the  high  noise  interfered  with  the  signal, 
causing  erroneous  V' s  to  be  indicated.  Although  not  available  at 
the  time  of  these  vibration  tests  in  September  1980,  later  models 
of  the  flow  velocimi  ter  include  better  filtering  and  bad  data  re¬ 
jection  algorithms,  which  arc  designed  to  provide  be  tter  noise' 
immunity.  The  resonances  themselves,  however,  need  to  be 
addressed  in  terms  of  mechanical  design,  to  minimize  their 
occurrence  or  strength,  rather  than  rely  on  electronic  or  soft¬ 
ware  remedies. 
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PHOTOGRAPH  OF  AXIS  DESIGNATION 


Z  AXIS 
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FURTHER  COMMENTS  ON  VIBRATION  TEST  PROCEDURES 

AND  RESULTS 


1.  Maximum  Flow  Rate 

The  maximum  water  flow  rate  through  the  cell  was  limited  to  3  ft/s 
(~1  m/s)  since,  at  ATL,  the  source  was  a  city  water  pressure  tap  some 
50  ft  (~15  m)  from  the  Ling  vibration  equipment.  Garden  hoses  were  used 
to  bring  water  to  the  delivered  cell  to  connect  between  this  cell  and  the 
second  cell  at  rest  (distance  about  10  ft  or  3  m),  and  from  the  second  cell 
to  an  outdoor  drain.  Thus,  the  total  length  of  garden  hose  was  110  ft  (—30  m). 

2.  Maximum  g  Load 

On  4  and  5  September  1980  several  preliminary  vibration  tests  were 
conducted  at  ATL  using  three  split  collars  as  simple,  low-cost  mounts. 

These  collars  were  bolted  to  the  driver  plate  using  three  1/2  in.  -  13  steel 
bolts.  During  y-axis  vibration,  however,  the  bolts  failed  several  times. 
Therefore,  the  mount  was  redesigned  as  shown  in  the  ATL  report  photographs. 
Prior  to  bolt  failure,  the  maximum  vibration  level  sensed  by  an  acceler¬ 
ometer  on  the  cell  was  140  g's. 

It  is  noteworthy  that  the  flowmeter  transducers  and  electronics  op¬ 
erated  satisfactorily  up  to  this  140  g  level.  Also,  after  the  preliminary 
testing  was  terminated,  the  flow  cell  was  disassembled  and  examined  for 
internal  failures.  None  were  found.  The  cell  was  reassembled,  remounted 
in  the  new  fixture  and  subjected  to  the  final  vibration  tests  documented  in  the 
AT L  report. 

3.  Vibration- Induced  Noise 

At  a  few  discrete  frequencies,  not  necessarily  the  major  structural 
resonances,  the  noise  exceeds  the  signal  at  the  receiver  transducer.  This 
audio  frequency  noise  "gets  through"  despite  ultrasonic  bandpass  filtering, 
apparently  by  means  of  nonlinear  mixing  in  the  transducer  element.  The 
approach  to  eliminating  this  problem  (as  evidenced  in  the  recordings  on  pp. 
J-l6,  J-17)  would  be  to  determine  the  source  and  modify  its  dimensions 
and/or  its  supports  and  damping. 
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APPENDIX  K.  FINAL  SYSTEM  TEST  IN  FUEL  SUBSTITUTE  7024  BII 

On  31  March  1981  and  1  April  1981,  at  GE,  Lynn,  Massachusetts, 
AETD  Bldg,  40,  Test  Stand  76,  the  M  system  was  tested  in  three  runs 
covering  flow  ranges  from  95  to  31,450  lb/hr. 

The  minimum  flow  rate  was  determined  by  the  smallest  available 
turbine  flowmeter,  whose  calibration  extended  down  to  114  lb/hr.  The 
maximum  was  determined  by  the  piping  and  sensor  configuration,  i,  e.  , 
31,450  lb/hr  was  the  maximum  M  achievable.  In  runs  1  and  2,  the  M 
limit  of  ~  27,  000  lb/hr  was  imposed  by  the  ANC  16  turbine  sensor,  which 
probably  would  have  been  damaged  if  higher  M  rates  would  have  been 
attempted.  The  high  pressure  boost  pump  used  in  these  tests  was  manu¬ 
factured  by  Roth,  model  No.  31  TEH  9479  BF  (92  gpm,  250  psi,  760  TDH). 

The  fluid  was  7024  BII,  a  fuel  substitute  commonly  used  in  cali¬ 
bration  tests.  The  Re  ranged  from  ^100  to  ~  ]0^,  thereby  covering 
laminar,  transitional  and  turbulent  regimes.  Temperature  ranged  from 
62.  4  to  79.  8°F.  Pressure  ranged  from  20  to  220  psig.  No  vibration 
isolators,  static  mixers  or  flow  straightener  s  were  used.  See  Fig.  K1  , 
which  shows  several  views  of  the  test  arrangement  and  test  equipment. 

Table  K1  lists  the  reference  turbine  flowmeters  used  for  'he 
three  runs,  their  ranges,  test  dates  and  cell  parameters.  Table 
K2  lists  all  relevant  data.  Dual  or  multiple  entries  represent  the  extremes 
of  turbine  and  ultrasonic  variability  at  each  flow  rate,  or  difference  in 
readings  of  two  different  turbine  meters  at  a  given  flow  rate.  The  apparent 
larger  spread  of  ultrasonic  values  is  attributed  to  the  sound  beam's  inertia - 
less  response  to  turbulence,  compared  to  the  response  of  a  turbine  flow¬ 
meter  which  is  smoothed  by  virtue  of  the  turbine's  mass  and  readout  means. 

The  spread  of  ultrasonic  values,  or  jitter,  has  been  reported  pre¬ 
viously  by  Pedersen  etal.  (1^77)’"  and  independently  by  Schmidt  (1980).  ^ 
Schmidt' s  chart  recordings  demonstrate  that  the  jitter  decreases  markedly 
as  the  integration  or  averaging  time  increases  from,  say,  1  s  to  1 0  s  to  100  s, 
in  crude  oil  clamp-on  tests  in  a  24  inch  pipeline.  Pedersen  etal.  ,  quantified 
the  decrease  in  standard  deviation  as  a  function  of  integration  time  interval, 
based  on  ultrasonic  measurements  along  a  45°  path  in  a  24  inch  natural  gas 
pipeline.  Integration  times  were  0.5,  5,  50  and  500  minutes. 


:;:N.  E.  Pedersen,  J.  E.  Bradshaw,  L.  C.  Lynnworth  and  P.  R.  Morel,  293- 
318,  especially  E'ig.  7,  in:  L.  K.  Irwin  (ed.),  NBSSpec.  Publ.  484,  Proc. 
Symp.  on  Flow  in  Open  Channels  and  Closed  Conduits  (1977). 

.’;;:;T.R.  Schmidt,  Clamp-On  Ultrasonic  Flowmeters  -  Application  Considera¬ 
tions  and  Field  Test  Results,  Proc.  ISA  1980  Conference,  pp.  111-126, 
ISA,  ISBN  87664-491-4,  C.  I.  80-510  (Oct.  1980). 
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Schmidt  has  suggested  (private  communication,  April  1981)  that 
ultrasonic  measurements  in  a  more  nearly  axial  direction  than  his  clamp-on 
geometry  allowed,  ought  to  be  less  responsive  to  the  radial  components  of 
turbulence,  and  hence  steadier.  This  is  probably  true.  However,  in  tests 
with  the  axial  offset  flow  cells  used  in  the  present  programs,  e.  g.  ,  the  A  H 
tests,  jitter  increased  as  V  increased.  This  jitter  may  have  been  mainly 
due  to  inlet  and  outlet  turbulence  near  the  tees,  which  had  no  flow  straight- 
eners  or  other  flow  control  surfaces  in  the  acoustic  path.  Perhaps  the  use 
of  vanes  or  tubes,  especially  in  the  offset  path  near  the  inlet  (as  suggested 
by  Fig.  11)  will  prove  helpful  in  reducing  tui  bulence  and  the  associated 
jitter.  (See  pp.  G-10,  G- 1 1  for  axial  path  jitter  test  data.  ) 

The  reference  M  is  computed  by  the  Cox  system  based  on  a  selected 
turbine  sensor  and  the  measured  fluid  temperature.  Uncertainty  is  stated 
to  be  1/2%  of  reading,  presumably  three  standard  deviations.  Because 
different  turbine  sensors  were  used  in  each  run  above  600  lb/hr,  the  difference 
between  the  ultrasonic  M  and  the  reference  M  is  not  the  same  in  each  run. 

From  the  raw  M  data  it  is  clear  that  the  ultrasonic  M  readout  is 
proportional  to  the  actual  M  as  measured  by  the  various  turbine  sensors. 

The  constant  of  proportionality  was  determined  by  comparing  the  averages 
of  all  the  valid  ultrasonic  data  shown  in  column  4  of  Table  K2.  The  M 
values  of  column  4  are  plotted  versus  the  corresponding  average  turbine  M 
in  Fig.  K2.  In  this  plot,  representing  three  runs  covering  M  from  95  to 
3  1,450  lb/hr,  the  results  are  evidently  linear  with  respect  to  full  scale 
flow.  Deviations  from  linearity,  expressed  as  percent  of  reading,  are 
largest  in  the  bottom  tenth  of  the  meter' s  range,  i.  e.  ,  below  3000  lb/hr 
(Table  K3). 

Table  K2  also  lists  the  flow  velocity  V,  computed  from  the  Cox  M, 
fluid  temperature  T,  pipe  area  A,  and  specific  gravity  SG  =  0.767  -  0.000435 
(T  -  70fF).  V  =  Q/A;  Q  =  M/ p  =  M/  62.  4  SG  cu  ft/  hr  =  1 .  998  x  10"  3  M/ SG 
gpm.  V  was  also  measured  ultr asonically  in  a  vertical  axial  offset  flow  cell 
built  from  a  1  inch  schedule  40  pipe  x  12  inches  long  and  two  tees,  with 
transducers  at  each  end.  These  transducers  had  buffers  (similar  to  those  in 
the  zigzag  cell)  to  prevent  gas  trapped  in  the  top  tee  from  blocking  the 
acoustic  wave.  The  corrected  ultrasonic  V  is  plotted  versus  V  calculated 
from  the  turbine  data  in  Fig.  K3.  Remarks  on  linearity  are  similar  to  remarks 
for  the  zigzag  cell.  For  both  cells,  no  hysteresis,  temperature,  pressure 
or  drift  effects  were  observed. 


Besides  the  numerical  data  of  Table  K2,  oscillographic  data  showing 
waveforms  for  the  zigzag  and  axial  offset  cells  are  of  interest  (Fig.  K4). 
These  waveforms  illustrate  several  aspects  of  the  ultrasonic  measurements, 
such  as  basic  time  intervals  tj ,  t ^  and  At  =  -  tj;  waveform  stability  in 

the  zigzag  cell  at  the  maximum  flow  rate  achievable  in  the  test  stand  for 
the  given  piping  configuration  (31,450  lb/hr);  waveform  jitter  in  the  offset 
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cell  (without  flow  control  surfaces)  as  observed  between  18,000  and  25,000  lb/ hi 
The  similarity  and  stability  of  the  upstream  and  downstream  waveforms 
and  high  signal  to  noise  ratio  for  the  zigzag  cell  at  all  flow  rates  in  these 
M  tests  (maximum  flow  velocity  Vmax  =  26.  2  ft/s  or  8  m/s  in  the  1  inch 
square  nole)  indicates  that  the  present  design  is  usable  to  even  higher  flow 
rates.  (The  Foxboro  water  flow  tests  previously  demonstrated  that  flows 
up  to  at  least  38  ft/s  (11.6  m/s)  could  be  covered  with  the  present  1  inch 
square  hole  zigzag  design.  ) 

By  comparing  the  raw  ultrasonic  M  and  turbine  M  data  lor  7024  1311, 
one  can  generate  correction  factors  which  comprise  a  look-up  table  lineari¬ 
zation  algorithm.  This  was  programmed  into  the  delivered  equipment.  This 
algorithm  serves  a  function  equivalent  to  the  dynamic  profile  compensation 
program  that  the  model  6 000  velocimeter  uses  in  conventional  tilted  diameter 
interrogations.  For  example,  when  the  beam  diameter  d  is  much  less  than 
the  duct  diameter  D  and  if  the  turbulent  profile  follows  the  universal  velocity 
distribution  law: 


Re 

$ 

2000 

K  - 

0.  750 

2000  <  Re 

< 

4000 

K  - 

0.  840 

Re 

4000 

K  = 

1/(1.119  -  0.  01 1  log  Re) 
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Table  Kl. 

Turbine  flowmeter 

and  ultrasonic  flow  c 

ell  parameter  data. 

Cox  turbine 

Date 

am  or 

Run 

Range, 

flowmeter 

used 

pm 

no. 

-PPh 

ANC  8-6, 

31  March 

am 

1  & 

114  to 

S/N23030 

1981 

pm 

2 

2280 

ANC  16, 

31  March 

pm 

2 

600  to 

S/N  10382 

1981 

24,000 

ANC  24, 

31  March 

am 

1 

990  to 

S/N23803 

1981 

59, 000 

ANC  24, 

1  April 

am 

3 

990  to 

S/N23801 

1981 

59,000 

Ultrasonic  cell 

L 

mm 

P 

mm 

A 

mm^ 

lw 

±1 

X 

mm 

Flow  axis 

orientation 

Square  hole 
zigzag 

102 

274.  3 

651.  6 

10 

30 

Horizontal 

Axial  offset 

353.  3 

404.  8 

557.  6 

6.  5 

20 

V  ertical 
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Table  K2,  continued.  Run 
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Table  K3. 


Rms  deviation  from  linearity  expressed  as 
percent  of  reading. 


Range  of 


M,  lb/hr 

Run  1 

Run  2 

Run  3 

Average 

0-2999 

11.  77 

6.  22 

- 

7.  91 

3000-5999 

0.  68 

1.  69 

1.  3 

1.  29 

6000-8999 

- 

1.  63 

1.  7 

1.  66 

9000-11999 

0.  32 

0.  43 

0.  71 

0.  51 

12000-14999 

- 

0 

0.  6 

0.  6 

15000-17999 

1.  5 

0.  2 

0.  2 

0.  77 

18000-20999 

- 

0.  75 

0 

0.  65 

21000-23999 

- 

1.  4 

0.  16 

00 

d 

24000  - 2699q 

- 

1.  27 

0.  08 

0.  90 

27000-29999 

- 

0.  7 

0.  41 

0.  53 

30000-32999 

L  ... 

- 

- 

0.  4 

0.  4 

Different  turbine  sensors,  each  normally  accurate  to  +  0.  5%  of  reading, 
were  used  above  600  lb/hr  in  each  run.  See  Table  Kl. 
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Fig.  K-l.  Mass  flowmeter  under  final  fuel  substitute  (7024  B  II)  test  at 
GE-I.ynn.  (a)  Genera!  view,  (b)  Close-up  of  flow  cell,  other  sensors 
and  piping,  (c)  Operator  reading  mass  flow  rate,  (d)  Close-up  of  mass 
flow  rate  instrument  (left)  and  volumetric  flow  rate  instrument  associated 
with  vertical  axial  offset  cell.  Dates  of  test:  31  March  and  1  April  1981. 
Photos  courtesy  of  General  Electric,  Lynn,  MA. 
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Fig.  K2.  Corrected  Panamctrics  mass  flow  rate  (M)  versus  Cox  turbine 
flowmeter  mass  flow  rate,  runs  1,  2  and  3.  Dates  of  tests: 

31  March  and  1  April  1981,  Site;  GE-Lynn. 
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V  CALCULATED  FROM  TURSINE  M,  FT/S 


Fig.  K3.  Corrected  ultrasonic  V  versus  value  calculated  from  turbine 
flowmeters.  Dates  of  tests:  31  March  and  1  April  1981. 
Site:  GE-Lynn, 
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Cell:  Zigzag  square  hole 

M  =  3 1, 480  pph 

(14,  309  kg/hr) 

t  =  220(j.s;  At  =  1.  2ps 

T  =  62.  4°F 


Cell:  Axial  offset 

M  =  18,  000  pph 
(8182  kg/hr) 

t  =  300  ps;  =  2.  3 1 


Fig.  K4.  Oscillograms  of  received  ultrasonic  signals  in  zigzag  and 
offset  cells,  using  7024  BI1  fluid  at  GE-Lynn.  IDate  of  test:  31  March  1981. 


NADC- 80254- 60 


NAVY  CALIBRATION  CORRECTION 


The  following  procedure  has  been  implemented  in  the  software  of  the  Navy 
flowmeter  to  correct  the  flow  measurement  using  calibration  data:* 


V  =  V  x  K. 

corrected  x 


where 


2  /t  -  T  1 
2  x  P  up  down 

L  X  (T  +  T 

up  down 


w 


V.  <  V  <  V.  , 
1  1  +  1 


V. 

K. 

i_ 

_ l 

_ i 

1 

0 

1.  000 

2 

.2 

.  970 

3 

.4 

.  940 

4 

.  6 

.  920 

5 

.  8 

.  900 

6 

1.0 

.  890 

7 

1.2 

.  885 

8 

1.4 

.  880 

9 

1.  6 

.  875 

10 

1.  8 

.  868 

1 1 

2.0 

.  865 

12 

3.  0 

.  850 

13 

4.0 

.  847 

14 

5.  0 

.  847 

15 

6.  0 

.  850 

16 

7.  0 

.  855 

17 

0 

.  860 

18 

7.  0 

.  862 

19 

10.0 

.  865 

20 

11.0 

.  867 

21 

12.  0 

.  869 

22 

13.  0 

.  870 

^Calibration  data 

obtained  using 

7024  BII  fuel  substitute 

at  GE,  I 

on  3 1  March  and 

1  April  1981. 

Data  in  NADC-80254-60, 

MA 
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APPENDIX  L.  NAYY  MASS  FLOWMETER  INITIAL  SET-UP. 

1)  The  mass  flowmeter  consists  of  two  electronic  circuits.  The  first 

is  the  Simmonds  Precision  electronics  associated  with  the  SP  densitometer 
cell.  The  second  is  the  Panametrics'  ultrasonic  flowmeter  associated  with 
the  zigzag  flow  cell. 

2)  SP  connections: 

a)  The  SP  density  flow  cell  has  two  connections  on  it.  One  is  marked 
H  and  the  other  is  unmarked.  There  are  two  cables  that  connect  the 
cell  to  the  SP  electronics.  The  connectors  are  keyed  at  their  con¬ 
nection  to  the  electronics.  The  cable  that  connects  to  the  noninsulated 
BNC  connector  connects  to  the  connector  labeled  H  at  the  cell.  The 
cable  connection  at  the  cell  for  H  also  has  a  length  of  white  shrink 
tubing  over  it  which  the  other  cable  does  not.  Both  cables  that  inter¬ 
connect  the  cell  and  electronics  are  white  and  both  must  be  connected 
for  correct  operation. 

b)  The  connection  of  the  SP  electronics  to  the  flowmeter  electronics  is 
do,.e  with  three  single  conductor  cables.  On  the  SP  unit  there  are 
three  binding  posts  labeled:  +15  (RED),  GND  (BLACK)  and  OUTPUT 
(YELLOW).  On  the  flowmeter  electronics  on  the  right  side  there  is 

a  12  position  double  row  terminal  block  (see  Fig.  19).  The  bottom 
three  terminals  are  labeled:  +15V,  ANALOG  INPUT,  and  CIRCUIT 
GND.  There  are  connections  already  on  one  row  of  the  three  terminals 
that  connect  into  the  flowmeter.  The  SP  terminal  labeled  +15  connects 
to  the  flowmeter  terminal  labeled  +15V,  the  SP  terminal  labeled  GND 
connects  to  the  flowmeter  terminals  labeled  CIRCUIT  GND  and  the 
SP  terminal  labeled  OUTPUT  connects  to  the  flowmeter  terminal  labeled 
ANALOG  INPUT. 

3)  The  flowmeter  to  transducer  conections  are  made  with  two  coaxial  cables. 

At  the  flowmeter  the  cables  attach  to  the  BNC  connectors  labeled  XMTA/RCVB 
and  XMTB/RCVA.  At  the  flow  cell  the  cables  attach  to  the  transducers  with 
Microdot  connectors.  There  is  no  preference  for  which  cable  attaches  to  which 
transducer  unless  the  digital  display  is  used.  In  this  case  the  transducer 
located  upstream  with  respect  to  flow  through  the  flow  cell  is  connected  to  the 
BNC  terminal  at  the  flowmeter  labeled  XMTB/RCVA,  and  the  downstream 
transducer  to  the  BNC  connector  labeled  XMTA/RCVB. 

4)  The  line  power  connections  for  the  flowmeter  are  made  on  a  10  position, 
single  row  terminal  block  on  the  right  side  of  the  flowmeter.  The  bottom 
three  terminals  are  used  and  are  labeled  N,  L,  and  GND  .L  The  120  VAC 

rs*  ? 

connections  are  made  to  N  and  L  while  the  earth  ground  connection  is  made 
to  the  GND  terminal. 
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